int cDiagonalA(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=0;j<nDim;j++) MassM[i*nDim+j]=va_arg(ap, COMPLEX);
va_end(ap);
if(idCur>idLIM) idCur=0;
id=idCur++;
extendData(id,nDim,4);
FError=FError|cJacobiA(MassM, nDim, ev[id], cU[id], cV[id]);
free(MassM);
return id;
}
int cjacobia_(Complex*a,int*n,double*d,Complex*u,Complex*v)
{ return cJacobiA(a,*n,d,u,v);}