To make it work I needed to find a better initial guess than the default (0.5,0.5). One working so far implementation of this idea is in the branch_codim_two branch. Unfortunately it needs some documentation and cleaning. It's in my TODO list but spare time is not on my side at the moment.
I assume your implementation is in fe/mapping_q1.cc in function MappingQ1::transform_real_to_unit_cell_internal? I.e., the one here:
http://dealii.org/svn/dealii/branches/branch_codim_two/deal.II/source/fe/mapping_q1.cc ? Best W. ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
