Ana,
how do you calculate the value of linear basis function and and order at the
given quadrature using if statement?
//Define basis functions
template <int dim>
double FEM<dim>::basis_function(unsigned int node, double xi){
/*"basisFunctionOrder" defines the polynomial order of the basis function,
"node" specifies which node the basis function corresponds to,
"xi" is the point (in the bi-unit domain) where the function is being
evaluated.
double value = 1.; //Store the value of the basis function in this variable
return value;
}
Are you asking how this is done in deal.II, or how to do it by hand?
In deal.II, we build 1d Lagrange polynomials of the correct order, and then
evaluate them as tensor products at the given point xi/eta.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
