getmatrix.h

Go to the documentation of this file.

#ifndef GETMATRIX_FILE
#define GETMATRIX_FILE
 
#include "geom.h"
#include <vector>
 
void CalcElem(int const ial[3], double const xc[], double ske[3][3], double fe[3]);
 
void AddElem(int const ial[3], double const ske[3][3], double const fe[3],
             int const id[], int const ik[], double sk[], double f[]);
 
 
// #####################################################################
class CRS_Matrix
{
    public:
       explicit CRS_Matrix(Mesh const & mesh);
 
       ~CRS_Matrix()
       {}
 
       void Derive_Matrix_Pattern();
 
       void CalculateLaplace(std::vector<double> &f);
 
       void ApplyDirichletBC(std::vector<double> const &u, std::vector<double> &f);
 
       void GetDiag(std::vector<double> &d) const;
 
       void Mult(std::vector<double> &w, std::vector<double> const &u) const;
 
       void Defect(std::vector<double> &w,
                   std::vector<double> const &f, std::vector<double> const &u) const;
 
       int Nrows() const
          {return _nrows;}
 
       void Debug() const;
 
       int fetch(int row, int col) const;
 
       void AddElem_3(int const ial[3], double const ske[3][3], double const fe[3], std::vector<double> &f);
 
       bool Compare2Old(int nnode, int const id[], int const ik[], double const sk[]) const;
 
    private:
       Mesh const & _mesh;      
       int _nrows;              
       int _nnz;                
       std::vector<int> _id;    
       std::vector<int> _ik;    
       std::vector<double> _sk; 
 
};
 
 
#endif