islre.c

Go to the documentation of this file.

/* src/islre.c
 *
 * Copyright (C) 2011-2018 Dongdong Li
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundataion; either version 3 of the License, or (at
 * your option) any later version.
 *
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABLITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
 * General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
 */

#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <getopt.h>
#include <time.h>
#include <math.h>
#include <stdbool.h>

#include "islre.h"
#include "copyright.c"
#include "usage.c"
#include "version.c"
#include "clear_comments.c"
#include "load_db.c"
#include "load_exp.c"
#include "test_load.c"
#include "test_models.c"
#include "psc_model.c"
#include "rkg_model.c"
#include "gem_ipopt.c"
#include "gsl_reg.c"
#include "print_results.c"

char *l_opt_arg;  
char* const short_options = "vh";  
struct option long_options[] = {
  { "version", 0, NULL, 'v'},  
  { "help", 0, NULL, 'h'},
  { 0, 0, 0, 0},  
};

int main (int argc, char *argvs[])
{
  int c;

  print_copyright_islre ();

  if (argc == 1)
    {
      print_re_usage (0);  
      return 100;  
    }

  while((c = getopt_long (argc, argvs, short_options, long_options, NULL)) != -1)  
    {  
      switch (c)  
    {  
    case 'h':  
      print_re_usage (0);  
      return 100;  
    case 'v':
      print_re_version ();  
      return 201;
    }
    }

  if (argc < 4)
    {
      printf ("invalid arguments\n");
      exit(0);
    }

  // the first step must be load databank
  if (load_db (argvs[1], TEMP_FILE) == false)
    {
      return 300;
    }

  fitted_param_num = atoi (argvs[3]);

  //test_load();

  if (load_exp (argvs[2]) == false)
    {
      return 301;
    }

  gsl_reg ();

  return 0;
}