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