[petsc-users] Adding adaptive mesh refinement to current fortran code

[TAY wee-beng]

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
>>