I have a non-symmetric matrix (attached). When I use in-core, I got wrong answer but occ produces correct answer. Can someone give me hints what's going on?
To test, I modified the pardiso_unsym_c.c (attached) to read my matrix in.
Exact solution should be ones. pardiso_unsym_c_sol_one_ic.res contains the solution using ic and pardiso_unsym_ooc.res using ooc.
I ran the example on linux using command:
make sointel64 interface=lp64 function=pardiso_unsym_c
PARDISO version is : 103000115
To test, I modified the pardiso_unsym_c.c (attached) to read my matrix in.
Exact solution should be ones. pardiso_unsym_c_sol_one_ic.res contains the solution using ic and pardiso_unsym_ooc.res using ooc.
I ran the example on linux using command:
make sointel64 interface=lp64 function=pardiso_unsym_c
PARDISO version is : 103000115
链接已复制
5 回复数
Hi there. I know it doesn't explain why you got different results on IC versus OOC, but I noticed that while you mention editing pardiso_unsym_c.c you actually attached pardiso_sym_c.c and it doesn't see to be edited to read your matrix. Could you reload your example program for IC that provides the wrong results? It would be helpful in finding any errant settings.Thanks!
Todd
yes, this the real bug.This issue has been submitted to our internal development tracking
database.
Here is a bug tracking number for your reference:DPD200169895. I will inform you once a new update becomes
available.--Gennady