<html>

  <head>

    <meta content="text/html; charset=windows-1252"

      http-equiv="Content-Type">

  </head>

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

    <br>

    <pre wrap=""> <i>Understood. But if you can define/describe the exact operator split method (i.e. what defines A and B) then we see how it could possibly be handled within TS. 

   Our goal is, when possible, to make TS flexible enough to support some operator split methods, we can only do this by specific examples.

   Thanks

    Barry</i></pre>

    <br>

    Au = (k grad(u)_x,  is a discretization of the flux in the x

    direction, while B is for flux in the y direction. These matrices

    are assembled by solving local elliptic problems in each cell

    following a  multiscale finite volume method.<br>

     To begin, the operator split is done as follows:<br>

    <br>

    u_t + Au = f ( considering only x dir fluxes) as (u* - u_n)/tau   +

    A u* = f and solve for u*<br>

    <br>

    then u_t + B u = f (considering only y dir fluxes) as (u - u*)/tau +

    B u = f and solve for u at current level.<br>

    <br>

    Since TS module cant be used, I shall have to do the time marching

    by writing my code. Perhaps I should be reusing the KSP solvers,

    since the sparsity pattern etc  will be unchanged at each time step.<br>

    <br>

    Thanks,<br>

    Shriram<br>

    <br>

  </body>

</html>