Пример #1
0
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;
}
Пример #2
0
int cjacobia_(Complex*a,int*n,double*d,Complex*u,Complex*v)
     { return  cJacobiA(a,*n,d,u,v);}