cuikerror.c

Go to the documentation of this file.

#include "box.h"
#include "random.h"
#include "defines.h"
#include "filename.h"
#include "chart.h"

#include <stdlib.h>

int main(int argc, char **arg)
{ 
  if (argc>1)
    {
      Tparameters parameters;
      TAtlasBase wcs;

      Tfilename fparam;
      Tfilename fsamples;
      FILE *fs;
      char *solName;

      unsigned int i,nv;
      double *s;

      boolean end;
      int token;

      unsigned int rep;
      unsigned int ns;
      double *fullS;
      
      /*Init parameters*/
      CreateFileName(NULL,arg[1],NULL,PARAM_EXT,&fparam);
      fprintf(stderr,"Reading parameters                : %s\n",GetFileFullName(&fparam));
      InitParametersFromFile(GetFileFullName(&fparam),&parameters);
      
      /* Read the problem file (cuik or world) */
      CS_WD_INIT(&parameters,arg[1],&wcs);
      nv=CS_WD_GET_NUM_SYSTEM_VARS(&wcs);
      NEW(s,nv,double);
      
      /* Open the the sample file */
      if (argc>2)
        solName=arg[2];
      else
        solName=arg[1];
      rep=(unsigned int)(GetParameter(CT_REPRESENTATION,&parameters));
      if (rep==REP_JOINTS)
        CreateFileName(NULL,solName,NULL,JOINTS_EXT,&fsamples);
      else
        CreateFileName(NULL,solName,NULL,LINKS_EXT,&fsamples);
      fprintf(stderr,"Reading samples from              : %s\n",GetFileFullName(&fsamples));
      fs=fopen(GetFileFullName(&fsamples),"r");
      if (!fs)
        Error("Could not open the file with the samples");

      end=FALSE;
      ns=1;
      do {
        token=fscanf(fs,"%lf",&(s[0]));
        if ((token==EOF)||(token==0))
          end=TRUE;
        else
          {
            for(i=1;i<nv;i++)
              {
                if (fscanf(fs,"%lf",&(s)[i])==EOF)
                  Error("Not enough values in the sample");
              }
            CS_WD_REGENERATE_SOLUTION_POINT(&parameters,s,&fullS,&wcs);
            fprintf(stderr,"  Error %u: %.16f\n",ns,CS_WD_ERROR_IN_EQUATIONS(fullS,&wcs));
            free(fullS);
            ns++;
          }
      } while (!end);

      free(s);
      fclose(fs);
      DeleteFileName(&fsamples);

      CS_WD_DELETE(&wcs);
      DeleteParameters(&parameters);

      DeleteFileName(&fparam);
    }
  else
    {
      fprintf(stderr,"Use:\n");
      fprintf(stderr,"  cuikerror <problem_name> [<sol_file>]\n");
      fprintf(stderr,"    <problem_name> the name from which to define the  .world,\n");
      fprintf(stderr,"    <sol_name> [optional] file with the the solution points\n");
    }
 
  return(EXIT_SUCCESS);
}