
<html>

  <head>


    <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">

  </head>

  <body bgcolor="#FFFFFF" text="#000000">

    Hi all,<br>

    I am implementing a time marching loop (My problem is not amenable 

    to usage of the TS module) where the coefficient matrix A wont

    change with time. I was wondering, in the sequence of commands<br>

    <br>

    <div align="justify"><i>KSPCreate(PETSC_COMM_SELF, &ksp); </i><i><br>

      </i><i>KSPSetOperators(ksp, A,  A, SAME_PRECONDITIONER); </i><i><br>

      </i><i>KSPGetPC(ksp, &pc);</i><i><br>

      </i><i>PCSetType(pc, PCLU); </i><i><br>

      </i><i>PCFactorSetMatSolverPackage(pc, MATSOLVERMUMPS);</i><i><br>

      </i><i>KSPSetFromOptions(ksp);</i><i><br>

      </i><i>KSPSolve(ksp, b, x);</i><i><br>

      </i><i>KSPDestroy(&ksp);</i><br>

      <br>

      would it be  efficient/possible to move some of them out of the

      loop.  <br>

      If so, which ones can I move out.<br>

      <br>

      Thanks,<br>

      Shriram<br>

      <br>

    </div>

  </body>

</html>