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>
