115 c_vector<double, ELEMENT_DIM>& rPhi,
118 c_vector<double, ELEMENT_DIM*PROBLEM_DIM> ret;
119 c_vector<double, PROBLEM_DIM> neumann_bc_values;
121 for (
unsigned i=0; i<ELEMENT_DIM; i++)
123 for (
unsigned problem_dim = 0; problem_dim<PROBLEM_DIM; problem_dim++)
125 double neumann_bc_value = this->mpBoundaryConditions->GetNeumannBCValue(&rSurfaceElement, rX, problem_dim);
127 ret(PROBLEM_DIM*i + problem_dim) = mScaleFactor * rPhi(i) * neumann_bc_value;