Hi,
I am attempting to calculate the curl of my computed solution, but have
been rather unsuccessful so far. I am attempting to use the
get_function_curls method in the FEValuesViews class to do so, but it
will not work for me. It compiles, but when I run it, I get an error
message:
408: undefined reference to `void dealii::FEValuesViews::Vector<2,
2>::get_function_curls<dealii::Vector<double> >(dealii::Vector<double>
const&, std::vector<dealii::Tensor<1, 1, double>,
std::allocator<dealii::Tensor<1, 1, double> > >&) const'
The segments of code that generate this error are:
std::vector<Tensor<1,1> > fe_function_curls(n_q_points, Tensor<1,1>());
...
[408] fe_views.get_function_curls(solution, fe_function_curls);
solution is the Vector<double> containing the coefficients of the
computed solution. I am using dim = 2, so Tensor<1,1> should be the
correct curl type.
What I really don't understand is that other methods in the
FEValuesViews class work just fine--for instance, I can call
get_function_divergences and get_function_gradients (with a
Tensor<2,dim> vector as the input, so this i the FEValuesViews method,
not the FEValues one). Can anyone provide some insight on this problem?
I have tried to compute the curl from the get_function_gradients method,
but because I'm not using first-order elements, I get a bogus result. I
can't figure out how to associate the entries from solution
(coefficients of global elements) to the local shape functions. If
someone could tell me that, I could compute the curls from the shape
functions (as I do when I assemble the system matrix).
Thank you!
Best,
-J
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
