3D implementation about multiphase-stokes HOT 6 OPEN

Sound good. Can we talk about this on Friday at our meeting so I can try to implement it?

Sure. We can discuss it on Friday. My suggestion is to finish up convergence testing so we can push out a paper before you try this.

from multiphase-stokes.

@bindi-nagda it would be very helpful for this if you wrote down the details of the smoother discretization somewhere (maybe your dissertation?)

from multiphase-stokes.

Yeah, I've included the details of the smoother discretization in my dissertation. I can send you a copy of it.

from multiphase-stokes.

Following up on this, initial tests indicate that we spend almost 90% of our time solving the linear system in the smoother. This suggests that this preconditioner definitely won't scale to 3D, as it will only get more expensive (solving a 13 by 13 system on each grid cell as opposed to a 9 by 9).

from multiphase-stokes.

Because we spend almost 90% of our time in the smoother, we may be able to significantly decrease run time by implementing our own LU or SVD solver in the smoother routine instead of using the lapack implementation. The compiler would then be free to optimize the matrix solve because the system size will be known at compile time.

from multiphase-stokes.

Because we spend almost 90% of our time in the smoother, we may be able to significantly decrease run time by implementing our own LU or SVD solver in the smoother routine instead of using the lapack implementation. The compiler would then be free to optimize the matrix solve because the system size will be known at compile time.

Sound good. Can we talk about this on Friday at our meeting so I can try to implement it?

from multiphase-stokes.