MUMPS interface function

zmumps - call to MUMPS

Calling Sequence

[id]=zmumps (id [,mat])

Input Parameters

Output Parameters

Internal Parameters

  • id.INST: (MUMPS reserved component) MUMPS internal parameter.
  • id.TYPE: (MUMPS reserved component) defines the arithmetic (complex or double precision).
  • Description

    The function zmumps solves systems of linear equations of the form Ax = b where A is square sparse matrix and b is a dense or sparse vector or matrix. The solver MUMPS is used and we refer the user to the MUMPS User's Guide for full details. Before the first call to zmumps, a call to initmumps must have been done:

         [id]=initmumps();
       

    Examples

    // this is a small linear system
    // whose solution is [1;2;3;4;5]
    A = sparse( [ 2  3  4  0  0;
                  3  0  -3  0  6; 
                  0 -1 1  2  0; 
                  0  0  2  0  0; 
                  0  4  0  0  1] );
    b = [20 ; 24; 9; 6; 13];
    
    // initialization of the MUMPS structure (here job=-1) 
    id=initmumps();
    [id]=zmumps(id);
    id.RHS=b;
    
    // call to MUMPS for the resolution
    id.JOB=6;
    [id]=zmumps(id,A);
    x=id.SOL
    norm(A*x-b)
    
    // Destruction of the MUMPS instance
    id.JOB=-2;
    [id]=zmumps(id);
    
       
       See also the examples provided in the directory "examples" that
       comes with the distribution of this interface.
       

    See Also

    initmumps,  dmumps,  

    References

    http://graal.ens-lyon.fr/MUMPS/

    http://www.enseeiht.fr/apo/MUMPS/