C++ (Cpp) rJacobiA示例

示例#1

文件： chDiag.c 项目： leloulight/micro_omega_4

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;
}

示例#2

文件： fortran.c 项目： restrepo/micromegas_old

int rjacobia_(double*a,int*n,double*d,double*u,double*v)
     {return rJacobiA(a,*n,d,u,v);}