this gives me the following cycle:
93 94: 0.119366
93 95: 0.301687
95 93: -0.242756
95 94: 0.493702
93,94,95 are the dof indices of one of the corners on the inside of
your mesh (the three components x,y,z).
In other words, we get:
x=ay+bz
z=cx+dy
I am not sure how to go from here. @Wolfgang, do you have an idea?
On Mon, Jun 25, 2012 at 1:38 PM, Timo Heister <[email protected]> wrote:
>> I can't print out the constraint matrix before constraints.close() is
>> called, correct? I tried and got an empty file.
>
> You need to make sure the file is flushed before you hit the
> exception. Just put the code in a block like this then it should work:
> {
> std::ostringstream constraint_filename;
> constraint_filename << "constraint_matrix.txt";
> std::ofstream constraint_output (constraint_filename.str().c_str());
> constraints.print(constraint_output);
> }
>
> --
> Timo Heister
> http://www.math.tamu.edu/~heister/
--
Timo Heister
http://www.math.tamu.edu/~heister/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii