int cDiagonalH(int nDim,...) { va_list ap; COMPLEX*MassM=malloc(sizeof(COMPLEX)*nDim*nDim); int i,j,id; va_start(ap,nDim); for(i=0;i<nDim;i++)for(j=i;j<nDim;j++) MassM[sMT(i,j)]=va_arg(ap, COMPLEX); va_end(ap); if(idCur>idLIM) idCur=0; id=idCur++; extendData(id,nDim,3); FError=FError|cJacobiH(MassM, nDim, ev[id], cV[id]); free(MassM); return id; }
int cjacobih_(Complex*a,int*n,double*d,Complex*v){return cJacobiH(a,*n,d,v);}