Invertering af en matrice vedhjælp af MKL

Jeg har ikke et kode eksempel, men jeg kan da prøve at konvertere
kaldene fra Fortran til C:

int n; // dimension
float a[n][n]; // matrix
int pivot[n]; // pivot
int status; // 0=success
float work[n]; // work storage

sgetrf(&n, &n, a, &n, pivot, &status);

sgetri(&n, a, &n, pivot, work, &n, &status);

a er både input og output.

Husk at C/C++ og Fortran har forskellig konvention omkring
to-dimensionelel arrays, så når du kalder fra C/C++ til Fortran så
kalder du errelt med den transponerede matrice.

Jeg har prøvet at indsætte din kode, men der kommer et forkert resultat ud. ??? Er det noget med c++ og fortrans forskellige håndtering af to-dimensionelel array

#include <stdio.h>
#include <mkl.h>
#include <math.h>
#include <float.h>
#include <stdlib.h>
int main()
{
const int dim=3;
int n=dim; // dimension
float a[dim][dim]; // matrix
int pivot[dim]; // pivot
int status; // 0=success
float work[dim]; // work storage
int i,j;

for(i=0;i<dim;i++)
{
    for(j=0;j<dim;j++)
    {
        a[i][j]=i*j+i+j*j+1 ;/*Bare nogle værdier */
        printf("%.2f ",a[i][j]);
    }
    printf("\n");
}


sgetrf(&n, &n, &a[0][0], &n, pivot, &status);
printf("\nsgetrf status :%i \n",status);
for(i=0;i<dim;i++)
{
    for(j=0;j<dim;j++)
    {
        printf("%.2f ",a[i][j]);
    }
    printf("\n");
}


sgetri(&n, &a[0][0], &n, pivot, work, &n, &status);
printf(" \nsgetri status :%i \n",status);
for(i=0;i<dim;i++)
{
    for(j=0;j<dim;j++)
    {
        printf("%.2f ",a[i][j]);
    }
    printf("\n");
}



return 0;

Hov vent lige, jeg tror at A bliver singulær, jeg kontrollere det lige.

jep, din kode funger. (O:
Nu virker det. Mange tak. (Husk at svar så jeg kan give dig point) (igen igen.)

svar