int rDiagonalA(int nDim,...) { va_list ap; REAL* MassM=malloc(sizeof(REAL)*nDim*nDim); int i,j,id; va_start(ap,nDim); for(i=0;i<nDim;i++) for(j=0;j<nDim;j++) MassM[i*nDim+j]=va_arg(ap,REAL); va_end(ap); if(idCur>idLIM) idCur=0; id=idCur++; extendData(id,nDim,2); FError=FError|rJacobiA(MassM, nDim, ev[id], rU[id],rV[id]); free(MassM); return id; }
int rjacobia_(double*a,int*n,double*d,double*u,double*v) {return rJacobiA(a,*n,d,u,v);}