getmatrix.h

Go to the documentation of this file.

#ifndef GETMATRIX_FILE
#define GETMATRIX_FILE

void Get_Matrix_Pattern(const int nelem, const int ndof_e, int const *const ia,
                        int &nnz, int* &id, int* &ik, double* &sk);



void CalcElem(const int ial[3], const double xc[], double ske[3][3], double fe[3]);

void AddElem(const int ial[3], const double ske[3][3], const double fe[3],
        const int id[], const int ik[], double sk[], double f[]);



void DebugMatrix(const int nnode, const int id[], const int ik[], const double sk[]);

void DebugVector(const int nnode, const double v[]);

void GetMatrix (const int nelem, const int ndof_e, const int ia[], const int nnode, const double xc[],
                const int nnz, const int id[], const int ik[], double sk[], double f[]);



void ApplyDirichletBC(const int nx, const int ny, const int neigh[],
       const double u[], const int id[], const int ik[] , double sk[], double f[]);



int fetch(const int row, const int col, const int id[], const int ik[]);


void Defect(double w[], const double f[], const double u[],
            const int nnode, const int id[], const int ik[], const double sk[]);


void CrsMult(double w[], const double u[], const int nnode, const int id[], const int ik[], const double sk[]);



void GetDiag(const int nnode, const int id[], const int ik[], const double sk[], double d[]);



#endif