- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Link Copied
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
It's part ofincompressible Navier-stokes solverfor multi-phase flow and is inconvenient to provide.
I have checked out some other error and now I can useMKL for my problemcorrectly.
Thank you very much!
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hello.
I'm having a similar problem when i'm trying to run one of the examples which i found in MKL example folder. The difference is, that when I use Neumman boundary conditions (NNNNNN) at the compilation i get the following error:
user1@user1:~/Documents/IntelFortran/mkl_poisson_solver$ make soia32 function=d_Helmholtz_3D_f
rm -f *.mod *.o
make mkl_dfti.mod mkl_trig_transforms.mod mkl_poisson.mod d_Helmholtz_3D_f.res _IA=ia32 EXT=so RES_EXT=so
make[1]: Entering directory `/home/sat/Documents/IntelFortran/mkl_poisson_solver'
mkdir -p ./_results/intel_parallel_ia32_so
ifort -mia32 -w -c /opt/intel/composer_xe_2013.1.117/mkl/include/mkl_dfti.f90
mkdir -p ./_results/intel_parallel_ia32_so
ifort -mia32 -w -c /opt/intel/composer_xe_2013.1.117/mkl/include/mkl_trig_transforms.f90
mkdir -p ./_results/intel_parallel_ia32_so
ifort -mia32 -w -c /opt/intel/composer_xe_2013.1.117/mkl/include/mkl_poisson.f90
mkdir -p ./_results/intel_parallel_ia32_so
ifort -mia32 -w source/d_Helmholtz_3D_f.f90 -L"/opt/intel/composer_xe_2013.1.117/mkl/lib/ia32" -Wl,--start-group "/opt/intel/composer_xe_2013.1.117/mkl/lib/ia32"/libmkl_intel.so "/opt/intel/composer_xe_2013.1.117/mkl/lib/ia32"/libmkl_intel_thread.so "/opt/intel/composer_xe_2013.1.117/mkl/lib/ia32"/libmkl_core.so -Wl,--end-group -L"/opt/intel/composer_xe_2013.1.117/mkl/../compiler/lib/ia32" -liomp5 -lpthread -ldl -o _results/intel_parallel_ia32_so/d_Helmholtz_3D_f.out
export LD_LIBRARY_PATH="/opt/intel/composer_xe_2013.1.117/mkl/lib/ia32":/opt/intel/composer_xe_2013.1.117/compiler/lib/ia32:/opt/intel/mic/coi/host-linux-release/lib:/opt/intel/mic/myo/lib:/opt/intel/composer_xe_2013.1.117/mpirt/lib/ia32:/opt/intel/composer_xe_2013.1.117/compiler/lib/ia32:/opt/intel/composer_xe_2013.1.117/mkl/lib/ia32:"/opt/intel/composer_xe_2013.1.117/mkl/../compiler/lib/ia32"; _results/intel_parallel_ia32_so/d_Helmholtz_3D_f.out >_results/intel_parallel_ia32_so/d_Helmholtz_3D_f.res
1
make[1]: *** [d_Helmholtz_3D_f.res] Error 1
make[1]: Leaving directory `/home/user1/Documents/IntelFortran/mkl_poisson_solver'
make: *** [soia32] Error 2
-----------------------------------------------------------------------------------------------------------------
and if i try to execute, i get the following message:
MKL POISSON LIBRARY WARNING:
The problem is degenerate due to rounding errors. The approximate solution
that provides the minimal Euclidean norm of the solution will be computed.
Double precision 3D Helmholtz example FAILED to compute the solution...
I have set the boundaries to zero:
bd_ax(iy,iz) = 0.0D0
bd_bx(iy,iz) = 0.0D0
bd_ay(ix,iz) = 0.0D0
bd_by(ix,iz) = 0.0D0
bd_az(ix,iy) = 0.0D0
bd_bz(ix,iy) = 0.0D0
and the q coefficient to zero, in order to create a poisson equation.
thank you in advance
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi,
The example check that error have to be 0 after each Helmholtz call. In your case error became equal to 1 that corresponded to warning about degenerate problem. To make your Helmholtz example passed correct change line
if (stat!=0) goto FAILURE;
on
if (stat<0) goto FAILURE;
Thank you,
Alexander Kalinkin
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page