ioset.c

二维热能方程的求解

/*** IO utility for reading setfile ***/
/*** Gererated by Dr. Kim's CSTART  ***/

#include <stdio.h>
#include <strings.h>
#include <unistd.h>
#if defined(_IBMR2) || defined (AIX) || defined(hpux)
#include <fcntl.h>
#else
#include <sys/fcntl.h>
#endif
#define MAX_LEN 80

typedef  double TYPE;
typedef  double ATYPE;
typedef  double VTYPE;
typedef  float DTYPE;

void setups(char*setfile,char*infofile,char*datafile,int*nt,int*nx,int*ny,VTYPE*at,VTYPE*bt,VTYPE*ax,VTYPE*bx,VTYPE*ay,VTYPE*by,int*itmax,VTYPE*tol,int*idELL,int*idsol,int*id_nl,int*level,VTYPE*eta,int*maxADI,int*maxGS,int*idSGS,int*maxPCG,int*idsymm,int*idpcond,int*interpol,int*id_hsrc,int*ntrue,int*ndiff,int*nonlin,int*nreact,int*nbc,int*idpsrc,VTYPE*wavet,VTYPE*wavex,VTYPE*wavey,VTYPE*tsbgn,VTYPE*tsend,int*n1,int*n2,int*ierr)
{
   char Program[MAX_LEN];
   char Author[MAX_LEN];
   char Date[MAX_LEN];
   char comment1[MAX_LEN];
   char comment2[MAX_LEN];
   char comments[MAX_LEN];
   char atype[MAX_LEN];
   char vtype[MAX_LEN];
   char array[MAX_LEN];
   char narr[MAX_LEN];
   char wksp[MAX_LEN];
   char nwksp[MAX_LEN];
   char dataext[MAX_LEN];
   char Headers[MAX_LEN];
   char Calls[MAX_LEN];
   char data[MAX_LEN];
   char ndata[MAX_LEN];

   FILE *fp000;

   if ((fp000=fopen(setfile,"r"))==NULL){
       printf("Error in SETUPS: when try to open <%s>\n",setfile);
       *ierr=1;return;}
   fscanf(fp000,"Program=%s Author=%s Date=%s nt=%d nx=%d ny=%d at=%lg bt=%lg ax=%lg bx=%lg ay=%lg by=%lg itmax=%d tol=%lg idELL=%d idsol=%d comment1=%s id_nl=%d level=%d comment2=%s eta=%lg maxADI=%d maxGS=%d idSGS=%d maxPCG=%d idsymm=%d comments=%s idpcond=%d interpol=%d id_hsrc=%d ntrue=%d ndiff=%d nonlin=%d nreact=%d nbc=%d idpsrc=%d wavet=%lg wavex=%lg wavey=%lg tsbgn=%lg tsend=%lg atype=%s vtype=%s array=%s narr=%s wksp=%s nwksp=%s dataext=%s Headers=%s Calls=%s",
          Program,Author,Date,nt,nx,ny,at,bt,ax,bx,ay,by,itmax,tol,idELL,idsol,comment1,id_nl,level,comment2,eta,maxADI,maxGS,idSGS,maxPCG,idsymm,comments,idpcond,interpol,id_hsrc,ntrue,ndiff,nonlin,nreact,nbc,idpsrc,wavet,wavex,wavey,tsbgn,tsend,atype,vtype,array,narr,wksp,nwksp,dataext,Headers,Calls);
   if(*level>=1) printf("SETUPS: from <%s>\n",setfile);
   if(*level>=3) {
      printf(" Program=%s Author=%s\n",Program,Author);
      printf(" Date=%s nt=%d\n",Date,*nt);
      printf(" nx=%d ny=%d\n",*nx,*ny);
      printf(" at=%lg bt=%lg\n",*at,*bt);
      printf(" ax=%lg bx=%lg\n",*ax,*bx);
      printf(" ay=%lg by=%lg\n",*ay,*by);
      printf(" itmax=%d tol=%lg\n",*itmax,*tol);
      printf(" idELL=%d idsol=%d\n",*idELL,*idsol);
      printf(" comment1=%s id_nl=%d\n",comment1,*id_nl);
      printf(" level=%d comment2=%s\n",*level,comment2);
      printf(" eta=%lg maxADI=%d\n",*eta,*maxADI);
      printf(" maxGS=%d idSGS=%d\n",*maxGS,*idSGS);
      printf(" maxPCG=%d idsymm=%d\n",*maxPCG,*idsymm);
      printf(" comments=%s idpcond=%d\n",comments,*idpcond);
      printf(" interpol=%d id_hsrc=%d\n",*interpol,*id_hsrc);
      printf(" ntrue=%d ndiff=%d\n",*ntrue,*ndiff);
      printf(" nonlin=%d nreact=%d\n",*nonlin,*nreact);
      printf(" nbc=%d idpsrc=%d\n",*nbc,*idpsrc);
      printf(" wavet=%lg wavex=%lg\n",*wavet,*wavex);
      printf(" wavey=%lg tsbgn=%lg\n",*wavey,*tsbgn);
      printf(" tsend=%lg atype=%s\n",*tsend,atype);
      printf(" vtype=%s array=%s\n",vtype,array);
      printf(" narr=%s wksp=%s\n",narr,wksp);
      printf(" nwksp=%s dataext=%s\n",nwksp,dataext);
      printf(" Headers=%s Calls=%s\n",Headers,Calls);
   }
   fclose(fp000);

   strcpy(infofile,"HEATSRC");
   strcpy(datafile,"HEATSRC@");
   if ((fp000=fopen(infofile,"r"))==NULL){
       printf("Error in SETUPS: when try to open <%s>\n",infofile);
       *ierr=1;return;}
   fscanf(fp000,"n1=%d n2=%d data=%s ndata=%s",n1,n2,data,ndata);
   if(*level>=1) printf("DATASETUP: from <%s>\n",infofile);
   if(*level>=3) {
      printf(" n1=%d n2=%d\n",*n1,*n2);
      printf(" data=%s ndata=%s\n",data,ndata);
      }
   fclose(fp000);

  return;
}