dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

I've started trying to remove the PETSc ifdefs from the ODE solvers,
but ran into some problems. The ODE solver uses the calls

    Vector x;
    real* xx = x.array();
    // operate on xx
    x.restore(xx);

This is supported only by PETSc vectors.

There are a couple of different solutions:

1. Don't use Vector::array().

2. Add Vector::array() to uBlas vectors.

3. Do 2 and add it to GenericVector as well.

Any thoughts? I have a feeling 1 is the best solution, but it's also
the one that requires the most work (for me to rewrite the ODE
solvers).

/Anders