1//---------------------------------------------------------------------
2//---------------------------------------------------------------------
3#include "header.h"
4
5void exact_solution(double xi,double eta,double zeta,double dtemp[]) {
6
7//---------------------------------------------------------------------
8//---------------------------------------------------------------------
9
10//---------------------------------------------------------------------
11//     this function returns the exact solution at point xi, eta, zeta
12//---------------------------------------------------------------------
13
14      int m;
15#define dtemp(m) dtemp[m-1]
16
17      for (m = 1; m <= 5; m++) {
18         dtemp(m) =  ce(m,1) +
19           xi*(ce(m,2) + xi*(ce(m,5) + xi*(ce(m,8) + xi*ce(m,11)))) +
20           eta*(ce(m,3) + eta*(ce(m,6) + eta*(ce(m,9) + eta*ce(m,12))))+
21           zeta*(ce(m,4) + zeta*(ce(m,7) + zeta*(ce(m,10) +
22           zeta*ce(m,13))));
23      }
24
25      return;
26}
27
28
29