[petsc-users] sqrt: DOMAIN error during run

TAY wee-beng zonexo at gmail.com
Mon Nov 21 13:10:50 CST 2011


Hi,

I have extended my 2D CFD code to 3D.

When I tried to do a step by step debug using Compaq visual Fortran 
(CVF) in windows, there is no error.

However, if I execute the code directly, the following error happens:

- sqrt: DOMAIN error
Image              PC        Routine            Line        Source
ibm2d_high_Re_wo_  00B95C69  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00B95AC7  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00B95C21  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00B9CEA8  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00BBF61D  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00BB7E30  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00BB1B98  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00B79319  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  0085C29A  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  0063C656  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  008EC67D  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  005C95DC  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00549D2C  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00500BEE  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00428FE3  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00463073  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00BC5309  Unknown               Unknown  Unknown
ibm2d_high_Re_wo_  00BB1035  Unknown               Unknown  Unknown
kernel32.dll       7C817067  Unknown               Unknown  Unknown

Incrementally linked image--PC correlation disabled.

I repeated running the code in Linux and calling "KSPGetConvergedReason" 
gives:

P Diverged

Does anyone know why this error occurs? It's strange that the step by 
step debug works and gives a reasonable solution but executing the code 
gives error.

Also, when changing from 2D to 3D, is there any additional stuff I need 
to add in besides adding another dimension?

-- 
Yours sincerely,

TAY wee-beng



More information about the petsc-users mailing list