91 setupArrays(
const Teuchos::RCP<const PointRule> & pr)
97 int num_nodes = point_rule->topology->getNodeCount();
98 int num_cells = point_rule->workset_size;
99 int num_space_dim = point_rule->spatial_dimension;
101 if (point_rule->isSide()) {
102 TEUCHOS_ASSERT(
false);
105 int num_points = point_rule->num_points;
107 coords_ref = af.template buildStaticArray<Scalar,IP,Dim>(
"coords_ref",num_points, num_space_dim);
109 node_coordinates = af.template buildStaticArray<Scalar,Cell,NODE,Dim>(
"node_coordinates",num_cells, num_nodes, num_space_dim);
111 jac = af.template buildStaticArray<Scalar,Cell,IP,Dim,Dim>(
"jac",num_cells, num_points, num_space_dim,num_space_dim);
112 jac_inv = af.template buildStaticArray<Scalar,Cell,IP,Dim,Dim>(
"jac_inv",num_cells, num_points, num_space_dim,num_space_dim);
113 jac_det = af.template buildStaticArray<Scalar,Cell,IP>(
"jac_det",num_cells, num_points);
115 point_coords = af.template buildStaticArray<Scalar,Cell,IP,Dim>(
"point_coords",num_cells, num_points, num_space_dim);