[Dune] constrained linear algebraic system Ax=b

[Carsten Gräser]

Hi,
you can only use those solvers and preconditioners
with square matrices. Please consult your favorite
book on numerical linear algebra for more details.

Carsten

Am 28.08.2015 um 19:22 schrieb 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 --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 473 bytes
Desc: OpenPGP digital signature
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20150901/cbdbe0af/attachment.sig>