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

[Barry Smith]

   I don't know. Never used paramesh.

   Barry


On Nov 18, 2010, at 4:04 AM, TAY wee-beng wrote:

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