Example of the CellTools class. More...
|void||vField (double &v1, double &v2, double &v3, const double &x, const double &y, const double &z)|
|Evaluation of a 3D vector field in physical coordinate frame. |
|int||main (int argc, char *argv)|
Example of the CellTools class.
Definition in file example_04.cpp.
Given a vector field u(x,y,z) and a face workset we want to compute the flux of u on every face in this workset. A face workset is a set of faces that are images of the same reference face. It is defined by the following items: 1. cell topology of a parent cell 2. a set of nodes in physical frame defining the parenct cells in the workset 3. subcell dimension and ordinal, relative to the reference cell in 1)
Given a face workset, the key steps to accomplish the task, , are as follows: 1. Obtain cubature points on workset faces, i.e., in physical frame; 2. Obtain (non-normalized) face normals at cubature points on workset faces 3. Evaluate the vector field u(x,y,z) at cubature points on workset faces 4. Compute dot product of u(x,y,z) and the face normals, times the cubature weights
|void vField||(||double &||v1,|
|const double &||x,|
|const double &||y,|
|const double &||z|