<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">All,<br>
<br>
I have been working with SAMRAI and Petsc to develop an implicit hyperbolic solver. Both packages seem very well done -- I am new to both, however. The basic physics problem is the solution to the linear wave equation (or rather something very close to it).
I am using SAMRAI to track the mesh and Petsc to solve the resulting implicit time stepping problem. I am giving the MPIcomm object from SAMRAI to Petsc but have noticed that SAMRAI distributes the processors across the domain seemingly randomly. Petsc on
the other hand distributes them sequentially for a matrix. In other words, SAMRAI might give processors with adjacent rank numbers nonadjacent patches while Petsc always assigns the FIRST n rows to the rank 0 process (where n is the number of rows requested
by the rank 0 process), the NEXT m rows to the rank 1 process, etc. This means that a nice, diagonally banded matrix is no longer banded (because we have reordered the rows according to the SAMRAI domain decomposition). Naturally, the solution to the linear
problem is the same. However, I am concerned that I have lost the apparent diagonal structure of the matrix. Is this a problem for the linear solvers? I know that Petsc is used for implicit problems together with various meshing codes -- this is a common
thing. I would like to use the library in its intended fashion -- is the proper thing to do to simply fill the rows you own and forget the fact that the structure is no longer banded? Am I missing something?<br>
<br>
Thanks in advance for any help you might offer,<br>
<br>
Todd<br>
<br>
</div>
</body>
</html>