- 新着としてマーク
- ブックマーク
- 購読
- ミュート
- RSS フィードを購読する
- ハイライト
- 印刷
- 不適切なコンテンツを報告
Hello,
The attached simple program "ti.f90" gives
1 =[1.00000000000000000,1.00000000000000000] 0.00E+00 -> GOOD
3 =[3.00000000000000000,3.00000000000000000] 0.00E+00 -> GOOD
1/3 =[0.33333333333333331,0.33333333333333337] 0.56E-16 -> GOOD
4/3 =[1.33333333333333326,1.33333333333333348] 0.22E-15 -> GOOD
1+1/3=[1.33333333333333348,1.33333333333333348] 0.00E+00 -> BAD !
The compilation was done using
ifort -D_LINUX $R/include/mkl_interval.f90 ti.f90
-L$R/lib/32 -lmkl_lapack -lmkl_ia32
Is it a mkl/8.0.1/lib/32/libmkl_lapack.a library bug or
I'm wrong (I have not find any detailed documentation)
Thanks
Best regards
--
Fabrice
The attached simple program "ti.f90" gives
1 =[1.00000000000000000,1.00000000000000000] 0.00E+00 -> GOOD
3 =[3.00000000000000000,3.00000000000000000] 0.00E+00 -> GOOD
1/3 =[0.33333333333333331,0.33333333333333337] 0.56E-16 -> GOOD
4/3 =[1.33333333333333326,1.33333333333333348] 0.22E-15 -> GOOD
1+1/3=[1.33333333333333348,1.33333333333333348] 0.00E+00 -> BAD !
The compilation was done using
ifort -D_LINUX $R/include/mkl_interval.f90 ti.f90
-L$R/lib/32 -lmkl_lapack -lmkl_ia32
Is it a mkl/8.0.1/lib/32/libmkl_lapack.a library bug or
I'm wrong (I have not find any detailed documentation)
Thanks
Best regards
--
Fabrice
コピーされたリンク
1 返信
- 新着としてマーク
- ブックマーク
- 購読
- ミュート
- RSS フィードを購読する
- ハイライト
- 印刷
- 不適切なコンテンツを報告
Yes, this is a bug. Thank you for pointing it out.
We are currently working on this, and, hopefully,the stuation will be fixed in MKL 8.1.