Hi Barry, In that case, will it be very difficult to use paramesh together with PETSc?
I may thinking that as long as I can create a matrix from paramesh and store it as a PETSc format matrix, everything should be fine when I try to solve Ax=b. Is it the same for the parallel case too? Or is it better if I direct the questions to the paramesh developers? Thanks! Yours sincerely, TAY wee-beng On 17-11-2010 18:28, Barry Smith wrote: > libmesh has adaptive refinement and all kinds of cool stuff and lives on > top of PETSc. > > Barry > > But it is C++. > > > On Nov 17, 2010, at 8:04 AM, TAY wee-beng wrote: > >> Hi, >> >> I am thinking of adding adaptive mesh refinement to my current fortran code. >> I am using an immersed boundary code with staggered cartesian grid. I am >> using PETSc to solve my momentum equations and PETSc/hypre to solve my >> poisson equation. The code now runs in parallel. >> >> Since I'm programming in fortran, some other c++ packages may not be >> feasible. I am considering packages like paramesh and AGRIF. I wonder if it >> is possible to integrate them together with PETSc. Is there anything which I >> need to consider too? >> >> Thanks! >> >> -- >> Yours sincerely, >> >> TAY wee-beng >>
