<DIV>Hi</DIV>  <DIV> </DIV>  <DIV>I used DA object to write our code for our 3D Reacting flow solver and it was quite useful.</DIV>  <DIV>because of the reaction inside the domain,it takes different amount of time to solve stiff equations.In part of the domain,which you don't have reaction, solving equations will be fast,but in the part that you have reaction, it takes much longer time to solve the equation.</DIV>  <DIV>therefore,there is inherently load imbalance inside the domain.</DIV>  <DIV> </DIV>  <DIV>when you use the DA object to decompose the mesh,I could not not find a way to enforce this issue through the DA functions.</DIV>  <DIV> </DIV>  <DIV>it would be great if you can have a function like</DIV>  <DIV> </DIV>  <DIV>DACreate1d(MPI_Comm comm,DAPeriodicType wrap,PetscInt M,PetscInt dof,PetscInt s,<STRONG>PetscInt weight</STRONG>,PetscInt *lc,DA *inra)</DIV>  <DIV> </DIV>  <DIV>through using <STRONG>weight,
 </STRONG>or something like that, you can dynamically balance your domain as you march in time.(for example after each 50 iterations,you can time your code for each sub-domain,and then you can define weight=time_subdomain/time_max_subdomains.)</DIV>  <DIV> </DIV>  <DIV>Is it possible to have something like this?or do you have any suggestion?</DIV>  <DIV> </DIV>  <DIV>thanks</DIV>  <DIV> </DIV>  <DIV>mehdi</DIV>  <DIV> </DIV>  <DIV> </DIV>  <DIV> </DIV>  <DIV><BR> </DIV><p>

      <hr size=1>Never miss a thing.  <a href="http://us.rd.yahoo.com/evt=51438/*http://www.yahoo.com/r/hs"> Make Yahoo your homepage.</a>