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.
May I know if it is possible to integrate paramesh together with PETSc?
Are there people who manage to do that? Is there anything which I need
to consider too?
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?
Thank you very much.
--
Yours sincerely,
TAY wee-beng
|
