72
72
double value(int n, double *wt, Func<double> *u_ext[], Func<double> *u, Func<double> *v, GeomVol<double> *e, Func<double> **ext) const{
73
73
double result = int_grad_u_grad_v<double, double>(n, wt, u, v) / Reynolds;
75
result += int_u_v<double, double>(n, wt, u, v) / time_step;
79
77
Ord ord(int n, double *wt, Func<Ord> *u_ext[], Func<Ord> *u, Func<Ord> *v, GeomVol<Ord> *e, Func<Ord> **ext) const
81
79
Ord result = int_grad_u_grad_v<Ord, Ord>(n, wt, u, v) / Reynolds;
83
result += int_u_v<Ord, Ord>(n, wt, u, v) / time_step;