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

Reply via email to