dolfin team mailing list archive

["Garth N. Wells" <g.n.wells@xxxxxxxxxx>]

On Thu, 2006-06-01 at 10:39 +0200, DOLFIN wrote:
> One or more new changesets pushed to the primary DOLFIN repository.
> A short summary of the last three changesets is included below.
> 
> changeset:   1981:5dd79b41a3d85cde459ff80e85edd968aeb90f78
> tag:         tip
> user:        "Garth N. Wells <g.n.wells@xxxxxxxxxx>"
> date:        Thu Jun  1 09:33:27 2006 +0200
> files:       src/kernel/la/uBlasKrylovSolver.cpp
> description:
> Add uBlasKrylovSolver.cpp
> 
> 
> changeset:   1980:f2770086f902eff6b19f85f6367d45054a998f09
> user:        "Garth N. Wells <g.n.wells@xxxxxxxxxx>"
> date:        Thu Jun  1 09:31:16 2006 +0200
> files:       src/kernel/la/dolfin/uBlasKrylovSolver.h src/kernel/la/dolfin/uBlasLUSolver.h
> description:
> Add uBlasKrylovSolver and uBlasLUSolver.
> 
> 
> changeset:   1979:6fef8c669a6a21635a3ff05a341d0b4fa9d2ae11
> user:        "Garth N. Wells <g.n.wells@xxxxxxxxxx>"
> date:        Thu Jun  1 09:15:15 2006 +0200
> files:       src/kernel/la/DenseVector.cpp src/kernel/la/GMRES.cpp src/kernel/la/KrylovSolver.cpp src/kernel/la/LU.cpp src/kernel/la/Makefile.am src/kernel/la/Makefile.in src/kernel/la/PETScKrylovSolver.cpp src/kernel/la/PETScLU.cpp src/kernel/la/dolfin/DenseVector.h src/kernel/la/dolfin/GMRES.h src/kernel/la/dolfin/KrylovSolver.h src/kernel/la/dolfin/LU.h src/kernel/la/dolfin/Makefile.am src/kernel/la/dolfin/Makefile.in src/kernel/la/dolfin/PETScKrylovSolver.h src/kernel/la/dolfin/PETScLU.h src/kernel/la/dolfin/Preconditioner.h src/kernel/la/dolfin/dolfin_la.h
> description:
> Add prefix "PETSc" to PETSc-based solvers and data structures. "KrylovSolver" and "LU" are typedefs to "PETScKrylovSolver" and "PETScLU" if PETSc is enabled, otherwise they are aliases to uBlasKrylovSolver and uBlasLUSolver.
> 
> This should make things simple throughout the code. If a PETSc-specific feature is required, use PETScFoo and add an "#ifdef", otherwise just use Foo.
> 

I'm having some trouble with the base class LinearSolver which has the
virtual function

 virtual uint solve(const SomeMatrix& A, SomeVector& x, const
SomeVector& b) = 0;
  
The problem is that I'd like to derive different solvers from it which
accept different matrix and vector types. For example

  PETScLU::solve(const PETScMatrix& A, PETScVector& x, const
PETScVector& b)
  uBlasLUSolver::solve(const uBlasMatrix& A, uBlasVector& x, const
uBlasVector& b)
 
Any ideas how to do this? I think that the Barton-Nackman trick (where
the base class is a template) which we tested for matrices might do the
job, at the expense of the complication introduced by the templates. 

Garth

> 
> -------------------------------------------------------
> For more details, visit http://www.fenics.org/hg/dolfin
> 
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/cgi-bin/mailman/listinfo/dolfin-dev