Hi
Could you give some more information as to why it does not work? What are V and
S? Normally you would have a solution U = [V , S] and then:
fe[vc].get_function_values (U, val_v); // for the velocity part
fe[sc].get_function_vaiues (U, val_s); // for the concentration part
Also you should use an FEValues object to extract the data:
FEValues<dim> fe_values_ref (fe, qf_cell, update_values | update_gradients
| update_JxW_values);
const FEValuesExtractors::Vector displacements(0);
fe_values_ref[displacements].get_function_values (solution_n,
local_solution_values_face);
Andrew
On 29.03.2011, at 18:10, Heitor Pina wrote:
> Hi all
>
> I define a FESystem consisting of dim (velocity, say) + 1 (concentration,
> say) components as
>
> fe (FE_Q<dim>(2), dim,
> FE_Q(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
