Satyendra, the simple answer might be: use FEValues::shape_value_component() instead of FEValues::shape_value(). The long answer is in the document that Wolfgang linked to.
Best, Guido Satyendra Tomar wrote: > Hi, > > I'm puzzled over how to compute the RHS term with the test functions form > the Raviart-Thomas (RT) space. Such a problem arises when working with the > operator I - grad div (which results in the bilinear form being an H(div) > inner product), and the RHS function a vector in H(div) (at the discrete > level). > > I've seen the examples 8, 20, and 22. In 20, the RHS is a scalar function > (and the inner product is computed with scaler test function). In other two > examples, it is a vector but is composed of primitive FE, whereas what is > required here a non-primitive FE. Using the example 20, the left hand side > is rather straight forward. However, how should I proceed for the RHS? > > Bests, > Satyendra > > > > ------------------------------------------------------------------------ > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
