On 12/18/2016 01:21 AM, seyedali88 via deal.II User Group wrote:

    I don't recall, but you can always just print the individual elements of
    these
    objects in a loop.


This is exactly what I am trying now, but I receive the following as output:

INVERSE OF JACOBIAN
nan nan nan nan
nan nan nan nan

The loop I wrote is:

pcout << std::endl << "INVERSE OF JACOBIAN" << std::endl;
for (unsigned int i = 0; i < 2; ++i)
{
for (unsigned int j = 0; j < dim; ++j)
pcout << invJ[i][j] << " ";
pcout << std::endl;
}
pcout << std::endl;

Strange. Is it not filled yet or am I accessing it in a wrong manner?

Are you running in debug mode? Did you pass the appropriate update_* flag to FEValues when you asked for these objects? And what kind of cell are you on?

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
For more options, visit https://groups.google.com/d/optout.

Reply via email to