Markus,

Thank you! I was able to update using svn and now everything works perfectly.

Best,
-J

On 04/15/2012 03:17 AM, Markus Bürg wrote:
Hello Justin,

I have checked in a fix for this bug in the current svn version.

Best Regards,
Markus



On 15.04.2012 00:17, Justin Droba wrote:
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

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


--
-----------------------------------------------------------------------------
Justin Droba
Office: A111 Wells Hall
Ph.D. Candidate in Applied Mathematics
Michigan State University
Department of Mathematics
East Lansing, MI 48824
-----------------------------------------------------------------------------
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to