void SurfaceScalarGradientOperator<EvalT, Traits>:: evaluateFields(typename Traits::EvalData workset) { Intrepid2::Vector<MeshScalarT> Parent_Grad_plus(3); Intrepid2::Vector<MeshScalarT> Parent_Grad_minor(3); for (int cell=0; cell < workset.numCells; ++cell) { for (int pt=0; pt < numQPs; ++pt) { Intrepid2::Tensor<MeshScalarT> gBasis(3, refDualBasis,cell, pt,0,0); Intrepid2::Vector<MeshScalarT> N(3, refNormal,cell, pt,0); gBasis = Intrepid2::transpose(gBasis); // in-plane (parallel) contribution for (int node(0); node < numPlaneNodes; ++node) { int topNode = node + numPlaneNodes; // the parallel-to-the-plane term for (int i(0); i < numPlaneDims; ++i ){ Parent_Grad_plus(i) = 0.5*refGrads(node, pt, i); Parent_Grad_minor(i) = 0.5*refGrads(node, pt, i); } // the orthogonal-to-the-plane term MeshScalarT invh = 1./ thickness; Parent_Grad_plus(numPlaneDims) = invh * refValues(node,pt); Parent_Grad_minor(numPlaneDims) = -invh * refValues(node,pt); // Mapping from parent to the physical domain Intrepid2::Vector<MeshScalarT> Transformed_Grad_plus(Intrepid2::dot(gBasis, Parent_Grad_plus)); Intrepid2::Vector<MeshScalarT> Transformed_Grad_minor(Intrepid2::dot(gBasis,Parent_Grad_minor)); // assign components to MDfield ScalarGrad for (int j(0); j < numDims; ++j ){ surface_Grad_BF(cell, topNode, pt, j) = Transformed_Grad_plus(j); surface_Grad_BF(cell, node, pt, j) = Transformed_Grad_minor(j); } } } } for (int cell=0; cell < workset.numCells; ++cell) { for (int pt=0; pt < numQPs; ++pt) { for (int k(0); k< numDims; ++k){ grad_val_qp(cell, pt, k) = 0; for (int node(0); node < numNodes; ++node) { grad_val_qp(cell, pt, k) += surface_Grad_BF(cell, node, pt, k)* val_node(cell,node); } } } } }
void DOFGradInterpolation_noDeriv<EvalT, Traits>:: evaluateFields(typename Traits::EvalData workset) { //Intrepid Version: // for (int i=0; i < grad_val_qp.size() ; i++) grad_val_qp[i] = 0.0; // Intrepid::FunctionSpaceTools:: evaluate<ScalarT>(grad_val_qp, val_node, GradBF); for (std::size_t i=0; i < grad_val_qp.size(); ++i) grad_val_qp(i)=0.0; for (int cell=0; cell < workset.numCells; ++cell) for (int qp=0; qp < numQPs; ++qp) for (int dim=0; dim<numDims; dim++) for (int node=0 ; node < numNodes; ++node) grad_val_qp(cell,qp,dim) += val_node(cell, node) * GradBF(cell, node, qp, dim); }
void DOFVecGradInterpolation<EvalT, Traits>:: evaluateFields(typename Traits::EvalData workset) { // This is needed, since evaluate currently sums into //for (int i=0; i < grad_val_qp.size() ; i++) grad_val_qp[i] = 0.0; for (std::size_t cell=0; cell < workset.numCells; ++cell) { for (std::size_t qp=0; qp < numQPs; ++qp) { for (std::size_t i=0; i<vecDim; i++) { for (std::size_t dim=0; dim<numDims; dim++) { // For node==0, overwrite. Then += for 1 to numNodes. ScalarT& gvqp = grad_val_qp(cell,qp,i,dim); gvqp = val_node(cell, 0, i) * GradBF(cell, 0, qp, dim); for (std::size_t node= 1 ; node < numNodes; ++node) { gvqp += val_node(cell, node, i) * GradBF(cell, node, qp, dim); //grad_val_qp(cell,qp,i,dim) += val_node(cell, node, i) * GradBF(cell, node, qp, dim); } } } } } // Intrepid::FunctionSpaceTools::evaluate<ScalarT>(grad_val_qp, val_node, GradBF); }
void DOFVecGradInterpolation<PHAL::AlbanyTraits::Jacobian, Traits>:: evaluateFields(typename Traits::EvalData workset) { int num_dof = val_node(0,0,0).size(); int neq = num_dof / numNodes; for (std::size_t cell=0; cell < workset.numCells; ++cell) { for (std::size_t qp=0; qp < numQPs; ++qp) { for (std::size_t i=0; i<vecDim; i++) { for (std::size_t dim=0; dim<numDims; dim++) { // For node==0, overwrite. Then += for 1 to numNodes. ScalarT& gvqp = grad_val_qp(cell,qp,i,dim); gvqp = FadType(num_dof, val_node(cell, 0, i).val() * GradBF(cell, 0, qp, dim)); gvqp.fastAccessDx(offset+i) = val_node(cell, 0, i).fastAccessDx(offset+i) * GradBF(cell, 0, qp, dim); for (std::size_t node= 1 ; node < numNodes; ++node) { gvqp.val() += val_node(cell, node, i).val() * GradBF(cell, node, qp, dim); gvqp.fastAccessDx(neq*node+offset+i) += val_node(cell, node, i).fastAccessDx(neq*node+offset+i) * GradBF(cell, node, qp, dim); } } } } } // Intrepid::FunctionSpaceTools::evaluate<ScalarT>(grad_val_qp, val_node, GradBF); }