Hello
I have coupled two problems 2D Navier-Stokes and 1D evolution equation on 
boundary of 2D. Assembling and solving of 2D problem is OK, problem is when i 
try to call function for solving of 1D problem. I have searched in matrices of 
1D problem in time and they are all nontrivial, but it gives me vector of 
trivial solutions. Assembling of 1D problem is therefore OK, but in solving, 
there is some problem.
In method solve for 1D i call system_matrix2.n_nonzero_elements() and it gives 
me nonzero number(13).
But how is possible, that it gives me a trivial solution of 1D problem.

code of solve2
{
SparseDirectUMFPACK A_direct;
   A_direct.initialize(system_matrix2);
   A_direct.vmult (solution2, system_rhs2);
}

Thank You
Marek Capek

P.S.: Right hand side is initialized from previous time iteration, therefore it 
is zero. 
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to