triangle.morley.1.bas_fun.c

主要特点 * Flexible finite element space construction based template element library 基于模板单元库灵活的构造有限

/*************************************************************************
* 	tirangle.1.bas_fun.c : by R.Lie
*/

void lambda_1(const double * p, const double ** v, void * value)
{
	double * val = (double *)value;
	double area = (v[1][0] - v[0][0])*(v[2][1] - v[0][1])
		- (v[1][1] - v[0][1])*(v[2][0] - v[0][0]);
	val[0] = (v[1][0] - p[0])*(v[2][1] - p[1])
		- (v[1][1] - p[1])*(v[2][0] - p[0]);
	val[0] /= area;
	val[0] = 1.0 - 2*val[0];
};

void lambda_2(const double * p, const double ** v, void * value)
{
	double * val = (double *)value;
	double area = (v[1][0] - v[0][0])*(v[2][1] - v[0][1])
		- (v[1][1] - v[0][1])*(v[2][0] - v[0][0]);
	val[0] = (v[2][0] - p[0])*(v[0][1] - p[1])
		- (v[2][1] - p[1])*(v[0][0] - p[0]);
	val[0] /= area;
	val[0] = 1.0 - 2*val[0];
};

void lambda_3(const double * p, const double ** v, void * value)
{
	double * val = (double *)value;
	double area = (v[1][0] - v[0][0])*(v[2][1] - v[0][1])
		- (v[1][1] - v[0][1])*(v[2][0] - v[0][0]);
	val[0] = (v[0][0] - p[0])*(v[1][1] - p[1])
		- (v[0][1] - p[1])*(v[1][0] - p[0]);
	val[0] /= area;
	val[0] = 1.0 - 2*val[0];
};

void gradient_lambda_1(const double * p, const double ** v, void * value)
{
	double * val = (double *)value;
	double area = (v[1][0] - v[0][0])*(v[2][1] - v[0][1])
		- (v[1][1] - v[0][1])*(v[2][0] - v[0][0]);
	val[0] = -2.0*(v[1][1] - v[2][1])/area;
	val[1] = -2.0*(v[2][0] - v[1][0])/area;
};

void gradient_lambda_2(const double * p, const double ** v, void * value)
{
	double * val = (double *)value;
	double area = (v[1][0] - v[0][0])*(v[2][1] - v[0][1])
		- (v[1][1] - v[0][1])*(v[2][0] - v[0][0]);
	val[0] = -2.0*(v[2][1] - v[0][1])/area;
	val[1] = -2.0*(v[0][0] - v[2][0])/area;
};

void gradient_lambda_3(const double * p, const double ** v, void * value)
{
	double * val = (double *)value;
	double area = (v[1][0] - v[0][0])*(v[2][1] - v[0][1])
		- (v[1][1] - v[0][1])*(v[2][0] - v[0][0]);
	val[0] = -2.0*(v[0][1] - v[1][1])/area;
	val[1] = -2.0*(v[1][0] - v[0][0])/area;
};

/*
*  end of file
**************************************************************************/