cuiksystem.h

Go to the documentation of this file.

#ifndef CUIKSYSTEMH
#define CUIKSYSTEMH

#include "equations.h"
#include "variables.h"
#include "statistics.h"
#include "list.h"
#include "box_list.h"
#include "parameters.h"
#include "box.h"
#include "boolean.h"
#include "vector.h"
#include "csmapping.h"

#include <stdlib.h>


/************************************************************************************/
#define DEPTH_FIRST_SEARCH 1
/************************************************************************************/

/************************************************************************************/
#define BREADTH_FIRST_SEARCH 2
/************************************************************************************/

/************************************************************************************/
#define MINIMIZATION_SEARCH 3
/************************************************************************************/

/************************************************************************************/
/*#define CUT_POINT 0.5*/
#define CUT_POINT 0.45+0.1*randomDouble()
/************************************************************************************/

/************************************************************************************/
#define MPI_TREAT_BOX_TIMEOUT(cs) (cs->nvariables*SIMPLEX_TIMEOUT*10)
/************************************************************************************/

typedef struct { 
  boolean updated; 
  boolean empty;   
  Tstatistics st; 
  Tmapping *simplification; 
  /*  */
  unsigned int searchMode; 
  Tequation orig_eqMin; 
  Tequation eqMin; 
  /******************************************************************************/
  /* SIMPLIFIED AND DUMMIFIED SUB-SYTEM                                         */
  Tmapping *orig2sd; 
  Tequations equations; 
  Tvariables variables; 
  unsigned int nequations; 
  unsigned int nvariables; 
  boolean *systemVar; 
  boolean *notDummyVar; 
  unsigned int *varType; 
  /******************************************************************************/
  /* SIMPLIFIED SUB-SYTEM                                                       */
  Tmapping *orig2s; 
  Tequations simp_equations; 
  Tvariables simp_variables; 
  unsigned int simp_nequations; 
  unsigned int simp_nvariables; 
  unsigned int simp_nee; 
  Tequation **J; 
  /******************************************************************************/
  /* ORIGINAL SUB-SYTEM                                                         */
  Tequations orig_equations; 
  Tvariables orig_variables; 
  unsigned int orig_nequations;
  unsigned int orig_nvariables;
  boolean *orig_systemVar; 
  boolean *orig_notDummyVar;
  char **orig_varNames; 
} TCuikSystem;


/****************************************************************************/
/* Methods to initialize/define/re-define CuikSystems (constructors et al.) */
/****************************************************************************/

void InitCuikSystem(TCuikSystem *cs);


void InitCuikSystemFromFile(Tparameters *p,char *filename,TCuikSystem *cs);

void CopyCuikSystem(TCuikSystem *cs_dst,TCuikSystem *cs_src);

void CuikSystemMerge(Tparameters *p,TCuikSystem *cs1,TCuikSystem *cs2,TCuikSystem *cs);


boolean CmpBoxesEquation(void *b1,void *b2,void *cs);

void SetCSSearchMode(unsigned int sm,Tequation *eqMin,TCuikSystem *cs);

void AddTerm2SearchCriterion(double w,unsigned int v,double val,TCuikSystem *cs);

/****************************************************************************/
/*     Methods to add individual variables/equations to the cuik system     */
/****************************************************************************/

unsigned int AddVariable2CS(Tvariable *v,TCuikSystem *cs);

void AddEquation2CS(Tparameters *p,Tequation *eq,TCuikSystem *cs);


/****************************************************************************/
/*        Methods to access the data stored in the CuikSystem.              */
/*  These methods provide controlled access to the stored data so that the  */
/*     higher levels can not introduce inconsistancies in the CuikSystem    */
/****************************************************************************/

void GetCSVariables(Tvariables *vs,TCuikSystem *cs);

unsigned int GetCSNumVariables(TCuikSystem *cs);

unsigned int GetCSNumNonDummyVariables(TCuikSystem *cs);


void GetCSVariable(unsigned int n,Tvariable *v,TCuikSystem *cs);

void SetCSVariableRange(unsigned int n,Tinterval *r,TCuikSystem *cs);

unsigned int GetCSVariableID(char *name,TCuikSystem *cs);

unsigned int GetCSSystemVars(boolean **sv,TCuikSystem *cs);

void GetCSEquations(Tequations *eqs,TCuikSystem *cs);

unsigned int GetCSNumEquations(TCuikSystem *cs);

void GetCSEquation(unsigned int n,Tequation *eq,TCuikSystem *cs);

void GetCuikSystemJacobian(Tequation ***J,TCuikSystem *cs);

void DeleteCuikSystemJacobian(Tequation **J,TCuikSystem *cs);

/****************************************************************************/
/*                    Methods to solve the CuikSystem.                      */
/****************************************************************************/

unsigned int MaxReduction(Tparameters *p,unsigned int varMask,Tbox *b,TCuikSystem *cs);

boolean SampleCuikSystem(Tparameters *p,char *fname,Tlist *sb,
                         unsigned int nsamples,unsigned int ntries,
                         unsigned int ndof,TCuikSystem *cs);

boolean SampleCuikSystemInBox(Tparameters *p,char *fname,Tlist *sb,
                              unsigned int nsamples,unsigned int ntries,
                              unsigned int ndof,
                              Tbox *init_box,TCuikSystem *cs);

void SolveCuikSystem(Tparameters *p,unsigned int nsols,
                     boolean restart,char *fstate,
                     FILE *f_out,Tlist *sol,TCuikSystem *cs);

#if (_USE_MPI)

  void MPI_SolveCuikSystem(Tparameters *p,
                           boolean restart,char *fstate,
                           FILE *f_out,TCuikSystem *cs);
  void MPI_TreatBox(Tparameters *p,TCuikSystem *cs);
#endif


/****************************************************************************/
/*                           Auxiliary Methods.                             */
/****************************************************************************/
void GenerateInitialBox(Tbox *box,TCuikSystem *cs);

void ReGenerateOriginalBox(Tparameters *p,Tbox *boxS,Tbox *boxO,TCuikSystem *cs);

unsigned int ComputeSplitDim(Tparameters *p,Tbox *b,TCuikSystem *cs);

boolean PointInSystemBox(Tvector *v,Tbox *b,TCuikSystem *cs);

double ErrorInSolution(Tbox *b,TCuikSystem *cs);

double ErrorInInequalities(Tbox *b,TCuikSystem *cs);

boolean CoordInequalitiesHold(Tbox *b,TCuikSystem *cs);

/****************************************************************************/
/*               Output the information in the CuikSystem                   */
/****************************************************************************/

void PrintCuikSystem(Tparameters *p,FILE *f,TCuikSystem *cs);

void PrintCuikSystemWithSimplification(Tparameters *p,FILE *f,TCuikSystem *cs);


void SaveCuikSystemSimplification(Tparameters *p,FILE *f,TCuikSystem *cs);


/****************************************************************************/
/*                                 Destructor                               */
/****************************************************************************/

void DeleteCuikSystem(TCuikSystem *cs);

#endif