cuikdim.c

Go to the documentation of this file.

#include "filename.h"
#include "error.h"
#include "parameters.h"
#include "atlas.h"
#include "samples.h"


#include <stdio.h>
#include <math.h>

int main(int argc, char **arg)
{
  Tparameters parameters;
  TAtlasBase world;
  Tfilename fparam;
  double *configuration;
  unsigned int k;

  if (argc>1)
    { 
      /* Init parameters */
      CreateFileName(NULL,arg[1],NULL,PARAM_EXT,&fparam);
      fprintf(stderr,"Reading parameters                : %s\n",GetFileFullName(&fparam));
      InitParametersFromFile(GetFileFullName(&fparam),&parameters);
      
      /* Read the world from file */
      CS_WD_INIT(&parameters,arg[1],&world);

      /* Read start sample */
      ReadOneSample(&parameters,arg[1],CS_WD_GET_NUM_SYSTEM_VARS(&world),&configuration);

      /* Estimate the dimesion of the solution manifold*/
      k=CS_WD_MANIFOLD_DIMENSION(&parameters,configuration,&world);
      fprintf(stderr,"N_DOF = %u\n",k);

      /* Print the N_DOF determined in the atlas initialization */
      DeleteFileName(&fparam);
      DeleteParameters(&parameters);
      CS_WD_DELETE(&world);
    }
  else
    {
      fprintf(stderr,"Use:\n");
      fprintf(stderr,"  cuikdim <problem_name>  \n");
      fprintf(stderr,"    <problem_name> the input problem.\n");
    }

  return(EXIT_SUCCESS);
}