26.08.2015 Views

nonsymmetric dynamics

ML 5.0 Smoothed Aggregation User's Guide - Trilinos - Sandia ...

ML 5.0 Smoothed Aggregation User's Guide - Trilinos - Sandia ...

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Jacobi Point Jacobi. Specify damping with coarse: damping factor.<br />

Gauss-Seidel Point Gauss-Seidel. Specify damping with coarse: damping<br />

factor.<br />

symmetric Gauss-Seidel Symmetric point Gauss-Seidel. Specify damping with coarse:<br />

damping factor.<br />

Chebyshev Use degree coarse: sweeps Chebyshev polynomial as a solver.<br />

Hiptmair Two-stage Hiptmair smoother. Take coarse: node sweeps<br />

(or coarse: edge sweeps) iterations with method coarse:<br />

subsmoother type to solve the nodal (or edge) subproblem.<br />

Amesos-KLU Use Klu through Amesos. Coarse grid problem is shipped to<br />

proc 0, solved, and solution is broadcast<br />

Amesos-Superlu<br />

Use SuperLU, version 3.0, through Amesos.<br />

Amesos-UMFPACK<br />

Use Umfpack through Amesos. Coarse grid problem is shipped<br />

to proc 0, solved, and solution is broadcasted.<br />

Amesos-Superludist Use SuperLU dist through Amesos.<br />

Amesos-MUMPS<br />

Use double precision version of Mumps through Amesos.<br />

user-defined User-defined smoothing function with prototype int<br />

(*)(ML Smoother *, int inLength, double *invec,<br />

int outLength, double *outvec). See ML example<br />

ml user smoothing.cpp.<br />

SuperLU<br />

Use ML interface to SuperLU(deprecated).<br />

Table 5: ML Epetra::MultiLevelPreconditioner: Commonly used coarse matrix solvers. Unless otherwise<br />

noted, number of sweeps is specified with coarse: sweeps. To use Amesos solvers, ML must be configured<br />

with with-ml amesos, and Amesos must be properly configured. In Gauss-Seidel and symmetric<br />

Gauss-Seidel, Ifpack implementation is automatically used (if available) for matrices of type Epetra<br />

CrsMatrix.<br />

// Zoltan requires coordinates<br />

double *xcoord, *ycoord, *zcoord;<br />

// user must load the coordinates, of coarse ...<br />

MLList.set("x-coordinates", xcoord);<br />

MLList.set("y-coordinates", xcoord);<br />

The following code demonstrates repartitioning for a system arising from the eddy current<br />

approximations to Maxwell’s equations.<br />

Teuchos::ParameterList MLList;<br />

// these correspond to the main edge matrix<br />

MLList.set("repartition: enable",1);<br />

MLList.set("repartition: max min ratio",1.3);<br />

MLList.set("repartition: min per proc",1000);<br />

MLList.set("repartition: partitioner","Zoltan");<br />

// Zoltan requires coordinates<br />

double *xcoord, *ycoord, *zcoord;<br />

20

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!