getmatrix.h File Reference

#include "geom.h"
#include <cassert>
#include <vector>

Go to the source code of this file.

Classes
class	Matrix
class	CRS_Matrix
class	FEM_Matrix
class	BisectInterpolation
	‍** More...
class	BisectIntDirichlet

Functions
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[])

Function Documentation

AddElem()

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[]
	)

Adds the element stiffness matrix ske and the element load vector fe of one triangular element with linear shape functions to the appropriate positions in the symmetric stiffness matrix, stored as CSR matrix K(sk,id, ik)

Parameters

[in]	ial	node indices of the three element vertices
[in]	ske	element stiffness matrix
[in]	fe	element load vector
[out]	sk	vector non-zero entries of CSR matrix
[in]	id	index vector containing the first entry in a CSR row
[in]	ik	column index vector of CSR matrix
[out]	f	distributed local vector storing the right hand side

Warning

Algorithm requires indices in connectivity ial in ascending order. Currently deprecated.

Definition at line 496 of file getmatrix.cpp.

CalcElem()

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

Calculates the element stiffness matrix ske and the element load vector fe of one triangular element with linear shape functions.

Parameters

[in]	ial	node indices of the three element vertices
[in]	xc	vector of node coordinates with x(2*k,2*k+1) as coordinates of node k
[out]	ske	element stiffness matrix
[out]	fe	element load vector

Definition at line 445 of file getmatrix.cpp.