Interface to calling the Jacobi function

[Nathan Shin]

Has anyone made a calling routine to call this Jacobi function?

My questions are :

1. How did you set up s and tau in the ROTATE definition;
2. How did you use **a to input the matrix


#define ROTATE(a,i,j,k,l) g=a[i][j];h=a[k][l];a[i][j]=g-s*(h+g*tau);a[k][l]=h+s*(g-h*tau);

void jacobi(float **a, int n, float d[], float **v, int *nrot)


Thank you in advance.

Nathan Shin