Quantcast
Channel: Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
Viewing all articles
Browse latest Browse all 2652

PARDISO b array

$
0
0

In most problems involving a matrix inversion, one only looks at a single b array in the Ax=b problem.

PARDISO has however an arrangement, where b is actually b[*,nhrs] in both the Geneva and the MKL versions. I checked the manuals.

But in the sample, just supplied by Intel, the b matrix in the C# example is called as a vector b[*] and passed to mkl as such. Clearly it runs and gives the correct answer, but why does the compiler not flag an error? This occurs in the other samples from other places, including the latest GENEVA sample

int mtype = 11; /* Real unsymmetric matrix */
        /* RHS and solution vectors. */
        double[] bData = new double[n];
        bData[101] = 450;
        AlignedDoubleArrayPointer b = new AlignedDoubleArrayPointer(bData);

 

Clearly not an ideal solution that one misses the second index on the array?

JMN


Viewing all articles
Browse latest Browse all 2652

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>