vector_structure.hpp

Go to the documentation of this file.

#pragma once

#include "./mpi_structure.hpp"

#include <climits>
#include <cmath>
#include <vector>
#include <iostream>
#include <iomanip>
#include <string>
#include <cstdlib>

using namespace std;

class CSysVector {
  
private:
  
    unsigned long nElm; 
    unsigned long nElmDomain; 
#ifdef HAVE_MPI
  unsigned long nElmGlobal; 
#endif
    unsigned short nVar; 
    unsigned long nBlk; 
    unsigned long nBlkDomain; 
  su2double* vec_val; 
public:
  
  CSysVector(void);
  
  CSysVector(const unsigned long & size, const su2double & val = 0.0);
  
  CSysVector(const unsigned long & numBlk, const unsigned long & numBlkDomain, const unsigned short & numVar, const su2double & val = 0.0);
  
  CSysVector(const CSysVector & u);
  
    void SetValZero(void);
  
  explicit CSysVector(const unsigned long & size, const su2double* u_array);
  
  explicit CSysVector(const unsigned long & numBlk, const unsigned long & numBlkDomain, const unsigned short & numVar,
                      const su2double* u_array);
  
  virtual ~CSysVector();
  
  void Initialize(const unsigned long & numBlk, const unsigned long & numBlkDomain, const unsigned short & numVar, const su2double & val = 0.0);
  
  unsigned long GetLocSize() const;
  
  unsigned long GetSize() const;
  
  unsigned short GetNVar() const;
  
  unsigned long GetNBlk() const;
    
  unsigned long GetNBlkDomain() const;
  
  void Equals_AX(const su2double & a, CSysVector & x);
  
  void Plus_AX(const su2double & a, CSysVector & x);
  
  void Equals_AX_Plus_BY(const su2double & a, CSysVector & x, const su2double & b, CSysVector & y);
  
  CSysVector & operator=(const CSysVector & u);
  
  CSysVector & operator=(const su2double & val);
  
  CSysVector operator+(const CSysVector & u) const;
  
  CSysVector & operator+=(const CSysVector & u);
  
  CSysVector operator-(const CSysVector & u) const;
  
  CSysVector & operator-=(const CSysVector & u);
  
  CSysVector operator*(const su2double & val) const;
  
  friend CSysVector operator*(const su2double & val, const CSysVector & u);
  
  CSysVector & operator*=(const su2double & val);
  
  CSysVector operator/(const su2double & val) const;
  
  CSysVector & operator/=(const su2double & val);
  
  su2double & operator[](const unsigned long & i);
  
  const su2double & operator[](const unsigned long & i) const;
    
  su2double norm() const;
  
  void CopyToArray(su2double* u_array);
  
  void SubtractBlock(unsigned long val_ipoint, su2double *val_residual);
  
  void AddBlock(unsigned long val_ipoint, su2double *val_residual);
  
  void SetBlock(unsigned long val_ipoint, unsigned short val_var, su2double val_residual);
  
  void SetBlock(unsigned long val_ipoint, su2double *val_residual);
  
  void SetBlock_Zero(unsigned long val_ipoint);
  
  void SetBlock_Zero(unsigned long val_ipoint, unsigned short val_var);
    
  su2double *GetBlock(unsigned long val_ipoint);
    
  su2double GetBlock(unsigned long val_ipoint, unsigned short val_var);
  
  
  friend su2double dotProd(const CSysVector & u, const CSysVector & v);
  
};

class CMatrixVectorProduct {
public:
  virtual ~CMatrixVectorProduct() = 0; 
  virtual void operator()(const CSysVector & u, CSysVector & v)
  const = 0; 
};
inline CMatrixVectorProduct::~CMatrixVectorProduct() {}

class CPreconditioner {
public:
  virtual ~CPreconditioner() = 0; 
  virtual void operator()(const CSysVector & u, CSysVector & v)
  const = 0; 
};
inline CPreconditioner::~CPreconditioner() {}

#include "vector_structure.inl"