interval.2.bas_fun.c

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

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

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

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

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

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

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

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

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