[petsc-users] sqrt: DOMAIN error during run

[TAY wee-beng]

Thanks Barry and Jed,

I found that the error is due to a divide by zero.

Some variables had zero values and were not calculated. Hence division 
by these values gave NaN. Strangely, running thru the debugger gave zero 
instead.

Yours sincerely,

TAY wee-beng


On 21/11/2011 8:33 PM, Barry Smith wrote:
>     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
>>