[petsc-users] sqrt: DOMAIN error during run

[Barry Smith]

   If it runs in one mode (like in a debugger) but crashes in another (like not in the debugger) this is often a symptom of memory corruption:  http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind

   Barry

On Nov 21, 2011, at 1:10 PM, TAY wee-beng wrote:

> 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
>