Hi,
I'm constantly getting an access violation error when I try to call Pardiso in a Visual Studio 2008 project. The basic structure of the code is:
integer value_count ivar, idum
real*8, allocatable :: myvalues(:)
integer, allocatable :: myia(:), myja(:)
INTEGER*4 pt(64)
INTEGER iparm(64)
!real*8, allocatable :: B2(:), x(:)
REAL(DP), allocatable :: B2(:), x(:)
INTEGER maxfct, mnum, mtype, phase, n, nrhs, error, msglvl
ALLOCATE(myvalues(value_count))
ALLOCATE(myja(value_count))
ALLOCATE(myia(mat_size+1))
ALLOCATE(x(mat_size))
do ivar = 1, 64
iparm(ivar) = 0
pt(ivar) = 0
end do
do ivar = 1, mat_size
!fill myvalues, myja, myia
...
x(ivar) = 0
end do
iparm(1) = 1
iparm(2) = 0
iparm(3) = 1
iparm(4) = 0
iparm(5) = 0
iparm(6) = 0
iparm(7) = 0
iparm(8) = 9
iparm(9) = 0
iparm(10) = 13
iparm(11) = 1
iparm(12) = 0
iparm(13) = 0
iparm(14) = 0
iparm(15) = 0
iparm(16) = 0
iparm(17) = 0
iparm(18) = -1
iparm(19) = -1
iparm(20) = 0
iparm(27) = 1
maxfct=1
mnum=1
nrhs=1
error = 0
msglvl = 1
mtype = -2
phase=13
!B is an input parameter of the function with type REAL(DP)
B2 = B
CALL pardiso (pt, maxfct, mnum, mtype, phase, mat_size, myvalues, myia, myja, idum, nrhs, iparm, msglvl, B2, x, error)
phase=-1
CALL pardiso (pt, maxfct, mnum, mtype, phase, mat_size, myvalues, myia, myja, idum, nrhs, iparm, msglvl, B2, x, error)
I'm working under IA32 architecture, linking and compiling works fine, only the memory error is bugging.
Does anybody have an idea?
