7 OutPut(
"Example: SinCos.h with INTERNAL_PROBLEM_IDENTITY " <<
14 void ExactU1(
double x,
double y,
double *values)
16 values[0] = sin(Pi*x);
17 values[1] = Pi*cos(Pi*x);
19 values[3] = -Pi*Pi*sin(Pi*x);
22 void ExactU2(
double x,
double y,
double *values)
24 values[0] = -Pi*y*cos(Pi*x);
25 values[1] = Pi*Pi*y*sin(Pi*x);
26 values[2] = -Pi*cos(Pi*x);
27 values[3] = Pi*Pi*Pi*y*cos(Pi*x);
30 void ExactP(
double x,
double y,
double *values)
32 values[0] = sin(Pi*x)*cos(Pi*y);
33 values[1] = Pi*cos(Pi*x)*cos(Pi*y);
34 values[2] = -Pi*sin(Pi*x)*sin(Pi*y);
35 values[3] = -Pi*Pi*sin(Pi*x)*cos(Pi*y)-Pi*Pi*sin(Pi*x)*cos(Pi*y);
38 void InitialU1(
double x,
double y,
double *values)
40 values[0] = sin(Pi*x);
43 void InitialU2(
double x,
double y,
double *values)
45 values[0] = -Pi*y*cos(Pi*x);
48 void InitialP(
double x,
double y,
double *values)
50 values[0] = sin(Pi*x)*cos(Pi*y);
56 void BoundCondition(
int i,
double t, BoundCond &cond)
62 void U1BoundValue(
int BdComp,
double Param,
double &value)
66 case 0: value=sin(Pi*Param);
68 case 1: value=sin(Pi);
70 case 2: value=sin(Pi*(1-Param));
74 default: cout <<
"wrong boundary part number" << endl;
80 void U2BoundValue(
int BdComp,
double Param,
double &value)
86 case 1: value=Pi*Param;
88 case 2: value=-Pi*cos(Pi*(1-Param));
90 case 3: value=-Pi*(1-Param);
92 default: cout <<
"wrong boundary part number" << endl;
101 void LinCoeffs(
int n_points,
double *X,
double *Y,
102 double **parameters,
double **coeffs)
107 double u1, u1x, u1y, u1lap, u2, u2x, u2y, u2lap, px, py;
109 for(i=0;i<n_points;i++)
121 u1lap = -Pi*Pi*sin(Pi*x);
122 u2 = -Pi*y*cos(Pi*x);
123 u2x = Pi*Pi*y*sin(Pi*x);
125 u2lap = Pi*Pi*Pi*y*cos(Pi*x);
126 px = Pi*cos(Pi*x)*cos(Pi*y);
127 py = -Pi*sin(Pi*x)*sin(Pi*y);
129 coeff[1] = -nu*u1lap+u1*u1x+u2*u1y+px;
130 coeff[2] = -nu*u2lap+u1*u2x+u2*u2y+py;
double RE_NR
Definition: Database.h:313
static TParamDB * ParamDB
Definition: Database.h:1134