[BACK]Return to phc6.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / phc

Diff for /OpenXM/src/phc/phc6.c between version 1.1 and 1.2

version 1.1, 1999/10/08 02:12:14 version 1.2, 2000/10/31 12:21:28
Line 58  char *phc_which(char *s);   /* search a path for the f
Line 58  char *phc_which(char *s);   /* search a path for the f
 struct phc_object phc_complexTo(long double r, long double i);  struct phc_object phc_complexTo(long double r, long double i);
   
   
 int phc_scan_for_string(FILE *fp, char str[], int lenstr);  int phc_scan_for_string(FILE *fp, char str[]);
 struct phc_object phc_scan_solutions(FILE *fp, int npaths, int dim );  struct phc_object phc_scan_solutions(FILE *fp, int npaths, int dim );
 struct phc_object phc_scan_output_of_phc(char *fname);  struct phc_object phc_scan_output_of_phc(char *fname);
 struct phc_object phc_call_phc(char *sys);  struct phc_object phc_call_phc(char *sys);
Line 77  main(int argc, char *argv[]) {
Line 77  main(int argc, char *argv[]) {
   int message = 0;    int message = 0;
   for (i=1; i<argc; i++) {    for (i=1; i<argc; i++) {
     if (strcmp(argv[i],"-v") == 0) {      if (strcmp(argv[i],"-v") == 0) {
       phc_verbose = 1;        phc_verbose = 1; message=1;
     }else if (strcmp(argv[i],"-g") == 0) {      }else if (strcmp(argv[i],"-g") == 0) {
       phc_overwrite = 0;        phc_overwrite = 0;
     }else if (strcmp(argv[i],"-i") == 0) {      }else if (strcmp(argv[i],"-i") == 0) {
Line 98  main(int argc, char *argv[]) {
Line 98  main(int argc, char *argv[]) {
   }    }
   while (1) {    while (1) {
     if (message) printf("dim= ");      if (message) printf("dim= ");
     if (scanf("%d",&dim)<0) break;          if (fgets(input,INPUTSIZE,stdin) <= 0) break;
       sscanf(input,"%d",&dim);
     sprintf(input,"%d\n",dim);      sprintf(input,"%d\n",dim);
     if (message) printf("Input %d equations please.\n",dim);      if (message) printf("Input %d equations please.\n",dim);
     for (i=0; i<dim; i++) {      for (i=0; i<dim; i++) {
Line 128  main(int argc, char *argv[]) {
Line 129  main(int argc, char *argv[]) {
   }    }
 }  }
   
 int phc_scan_for_string(FILE *fp, char str[], int lenstr)  int phc_scan_for_string(FILE *fp, char str[])
      /*       /*
   **  Scans the file fp for a certain string str of length lenstr+1.    **  Scans the file fp for a certain string str of length lenstr+1.
   **  Reading stops when the string has been found, then the variable    **  Reading stops when the string has been found, then the variable
   **  on return equals 1, otherwise 0 is returned.    **  on return equals 1, otherwise 0 is returned.
   */    */
 {  {
   char buf[lenstr+1];  #define BUF_SIZE 1024
     char buf[BUF_SIZE];
   char ch;    char ch;
   int index,i,compare,npaths,dim,found;    int index,i,compare,npaths,dim,found;
     int lenstr;
     lenstr = strlen(str);
     if (lenstr >= BUF_SIZE-1) {
           fprintf(stderr,"Too long string in phc_scan_for_string\n");
           exit(-1);
     }
   index = -1;    index = -1;
   found = 0;    found = 0;
   while ((fscanf(fp,"%c",&ch)!=EOF) && found == 0)    while (((ch = fgetc(fp))!=EOF) && found == 0)
     {  
       if (index == -1 && ch == str[0])  
         {          {
           index = 0;            if (index == -1 && ch == str[0])
           buf[index] = ch;  
         }  
       else  
         {  
           if (index == lenstr)  
             {  
               compare = 0;  
               for (i=0; i<lenstr+1; i++)  
                 {                  {
                   if (buf[i]!=str[i])                    index = 0;
                     {                    buf[index] = ch;
                       compare = compare+1;  
                     }  
                 }                  }
               if (compare == 0)            else
                 {                  {
                   found = 1;                    if (index == lenstr)
                           {
                             compare = 0;
                             for (i=0; str[i] != '\0'; i++)
                                   {
                                     if (buf[i]!=str[i])
                                           {
                                             compare = compare+1;
                                           }
                                   }
                             if (compare == 0)
                                   {
                                     found = 1;
                                   }
                             index = -1;
                           }
                     else
                           if (index > -1 && index < lenstr)
                             {
                                   index = index+1;
                                   buf[index] = ch;
                             }
                 }                  }
               index = -1;            if (found == 1) break;
             }  
           else  
             if (index > -1 && index < lenstr)  
               {  
                 index = index+1;  
                 buf[index] = ch;  
               }  
         }          }
       if (found == 1) break;  
     }  
   return found;    return found;
 }  }
 struct phc_object phc_scan_solutions(FILE *fp, int npaths, int dim )  struct phc_object phc_scan_solutions(FILE *fp, int npaths, int dim )
Line 187  struct phc_object phc_scan_solutions(FILE *fp, int npa
Line 195  struct phc_object phc_scan_solutions(FILE *fp, int npa
   struct phc_object rob,sob;    struct phc_object rob,sob;
   char ch;    char ch;
   int fnd,i,j,nsols;    int fnd,i,j,nsols;
   float res;    double res;
   long double realpart;    long double realpart;
   long double imagpart;    long double imagpart;
   long double realparts[npaths][dim];    long double realparts[npaths][dim];
   long double imagparts[npaths][dim];    long double imagparts[npaths][dim];
   nsols = 0;    nsols = 0;
   while (fscanf(fp,"%c",&ch)!=EOF)    while ((ch = fgetc(fp)) != EOF)
     {      {
       fnd = phc_scan_for_string(fp,"start residual :",15);        fnd = phc_scan_for_string(fp,"start residual :");
       if (fnd==1)        if (fnd==1)
         {          {
           fscanf(fp,"%E",&res);            fscanf(fp,"%E",&res);
           /* printf(" residual = "); printf("%E\n",res); */            /* printf(" residual = "); printf("%E\n",res);  */
           if (res < 1.0E-12) nsols = nsols+1;            if (res < 1.0E-12) {
           fnd = phc_scan_for_string(fp,"the solution for t :",19);                  nsols = nsols+1;
                   if (nsols > npaths) {
                     fprintf(stderr,"Something is wrong in phc_scan_solutions\n");
                     fprintf(stderr,"npaths=%d, nsols=%d \n",npaths,nsols);
                     exit(-1);
                   }
             }
             fnd = phc_scan_for_string(fp,"the solution for t :");
           for (i=0;i<dim;i++)            for (i=0;i<dim;i++)
             {              {
               fnd = phc_scan_for_string(fp,":",0);                fnd = phc_scan_for_string(fp,":");
               fscanf(fp,"%LE",&realpart);                fscanf(fp,"%LE",&realpart);
               fscanf(fp,"%LE",&imagpart);                fscanf(fp,"%LE",&imagpart);
               if (res < 1.0E-12)                if (res < 1.0E-12)
Line 245  struct phc_object phc_scan_output_of_phc(char *fname)
Line 260  struct phc_object phc_scan_output_of_phc(char *fname)
   int fnd,npaths,dim,i,nsols;    int fnd,npaths,dim,i,nsols;
   otp = fopen(fname,"r");    otp = fopen(fname,"r");
   if (phc_verbose) fprintf(stderr,"Scanning the %s of phc.\n",fname);    if (phc_verbose) fprintf(stderr,"Scanning the %s of phc.\n",fname);
   fnd = phc_scan_for_string(otp,"THE SOLUTIONS :",14);    fnd = phc_scan_for_string(otp,"THE SOLUTIONS :");
   fscanf(otp,"%i",&npaths);    fscanf(otp,"%i",&npaths);
   if (phc_verbose) fprintf(stderr,"  number of paths traced = %i\n",npaths);    if (phc_verbose) fprintf(stderr,"  number of paths traced = %i\n",npaths);
   fscanf(otp,"%i",&dim);    fscanf(otp,"%i",&dim);
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>