[petsc-users] PetscKernel_A_gets_inverse_A_

[Jed Brown]

Sorry about the slow response. Do you only want the square kernels with
static (compile-time) size? These are currently implemented as macros, but
it's not a problem to modify them to be functions. C99 provides a
well-defined mechanism to have a single version with external linkage, but
also encourage inlining.

I'm considering providing a set of kernels that would support non-square
matrices, but the naming conventions would have to change and my
applications would frequently not know the size statically (but it would
typically be ~10 or less, so calling BLAS doesn't make sense).


On Fri, Dec 21, 2012 at 3:29 AM, Aldo Bonfiglioli <
aldo.bonfiglioli at unibas.it> wrote:

> Dear all,
> would it be possible to have a unified interface (also Fortran callable)
> to the PetscKernel_A_gets_inverse_A_ routines?
> I find them very useful within my own piece
> of Fortran code to solve small dense linear system (which I have
> to do very frequently).
> I have my own interface, at present, but I need to
> change it as needed when a new PETSc version is released.
>
> Regards,
> Aldo
> --
> Dr. Aldo Bonfiglioli
> Associate professor of Fluid Flow Machinery
> Scuola di Ingegneria
> Universita' della Basilicata
> V.le dell'Ateneo lucano, 10 85100 Potenza ITALY
> tel:+39.0971.205203 fax:+39.0971.205215
>
>
> Publications list <http://publicationslist.org/aldo.bonfiglioli>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20121231/00889033/attachment.html>