← Back to team overview

dolfin team mailing list archive

Re: Status of linear algebra

 

Anders Logg wrote:
Garth, what is the status of the linear algebra? If you are happy,
I'll take another look and then we can take care of the remaining stuff.

It's pretty much done. I still need to fix the selection of the preconditioner when choosing by name,

  PETScKrylovSolver solver(PETScPreconditioner::ilu);

You'll see that there is an enum list in PETScPreconditioner and uBlasPreconditioner. I'd like to change this to just one list so it won't matter whether PETSc or uBlas is being used, and if the requested preconditioner is not available, then a default is chosen. Should then be possible to do

  KrylovSolver solver(Preconditioner::ilu);

Same story for choosing the Krylov method by name. Any suggestions?

I'm also thinking of moving the individual Krylov methods out of uBlasKrylovSolver.



Here's a suggested plan:

       1. Review/cleanup (Garth)
       2. Review/cleanup/class diagram (Anders)
       3. Fix PyDOLFIN linear algebra including typedefs (Johan J)
       3. Finish linear algebra section in manual (Johan H)
       4. Remove PETSc ifdefs from Navier-Stokes module (Johan H)

I did this yesterday. Just one PETSc ifdef left for choosing the Hypre preconditioner.

       4. Remove PETSc ifdefs from elasticity module (Johan J)

I removed them from elasticity (not elasticity-updates). The module needs some cleaning up still.

       4. Unit tests for linear algebra (Anders, Garth)
       4. Demos for linear algebra (Garth)
       4. Benchmarks for linear algebra (Garth)
       4. Review/cleanup of ODE solvers (Anders)

What more do we want to finish before the next release?


Sounds good.

Garth

Hopefully we can make a concentrated effort and make the release
quickly so we can move on to work on the mesh.

/Anders
_______________________________________________
DOLFIN-dev mailing list
DOLFIN-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/dolfin-dev



Follow ups

References