Hello Heitor,

can you give us the error message? Or specify a bit more detailled what does not work?

Thanks,
Markus



Am 29.03.11 18:10, schrieb Heitor Pina:

Hi all

I define a FESystem consisting of dim (velocity, say) + 1 (concentration, say) components as

fe (FE_Q  
<file:///C:%5CUsers%5CHitor%20Pina%5CDesktop%5CdealII700-doc%5Cdoxygen%5Cdeal.II%5CclassFE__RaviartThomas.html><dim>(2),
 dim,
     FE_Q  
<file:///C:%5CUsers%5CHitor%20Pina%5CDesktop%5CdealII700-doc%5Cdoxygen%5Cdeal.II%5CclassFE__DGQ.html>(2),
 1)

I wish to get function values at the quadrature points of a finite element velocity function V and a finite element scalar function S.

I thought that  issuing

const  FEValuesExtractors::Vector  vc(0);

const  FEValuesExtractors::Scalar   sc(dim);

fe[vc].get_function_values (V, val_v);    // for the velocity part

fe[sc].get_function_vaiues (S, val_s);    // for the concentration part

would do the trick but this does not seem to work. As I may be missing something very basic here any help would be appreciated.

--

TIA

        Heitor Pina


_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to