Is solving a Nx * Ny system equivalent to a Nx * Ny * 1 system in PETSc?

Barry Smith bsmith at
Wed Oct 14 12:41:33 CDT 2009

    Speed should be pretty much the same.

    I could interpret your question as: If one as a 3d PDE solver is  
it worth writing a totally separate 2d code or can one just run the 3d  
code with 2 dimensions? My answer to this question is, since code  
development and maintenance is very time-consuming it makes sense to  
just have one code.


On Oct 14, 2009, at 10:43 AM, Wee-Beng TAY wrote:

> Hi,
> I'm not sure if my question is trivial but I'm wondering if it makes  
> any difference in solving a Nx * Ny system and a Nx * Ny * 1 system  
> in PETSc? For e.g., 10x10 vs. 10x10x1. Will there be any problem or  
> difference in answer? What about the speed?
> -- 
> Thank you very much and have a nice day!
> Yours sincerely,
> Wee-Beng Tay

More information about the petsc-users mailing list