[Dune] constrained linear algebraic system Ax=b

[abdeladhim tahimi]

Hi dune developers,

I have the linear system A.x=b where:

A = [  1.0   -1.0    0.0 ; ...
      -1.0    2.0   -1.0 ; ...
       0.0   -1.0    1.0 ; ]

and,

b = [ 0.0; ...
      1.0; ...
      0.0; ]

and the condition x[1] = 0;

I added an equation to my system to get "A1.x=b1" with non-square matrix A1,

A1 = [  1.0   -1.0    0.0 ; ...
       -1.0    2.0   -1.0 ; ...
        0.0   -1.0    1.0 ; ...
        0.0    1.0    0.0 ; ]

b1 = [ 0.0; ...
       1.0; ...
       0.0; ...
       0.0; ]

I have not been able to solve this new system using BiCGSTABSolver and
SeqSSOR preconditioner, something is telling me I am taking a wrong way to
solve my original problem "A.x=b". any help?

thank you in advance for your hints.

Adim
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20150828/d87d73c4/attachment.htm>