test_c.c

Go to the documentation of this file.

#include "amns_interface.h"

int main(int argc, char *argv[])
{
  void* amns_handle = NULL; 
  amns_c_error_type error_stat = DEFAULT_AMNS_C_ERROR_TYPE;
  void* reactants_handle = NULL;
  amns_c_reactant_type species1  = {.ZN=6, .ZA=1, .MI=12, .LR=0}; 
  amns_c_reactant_type species2  = {.ZN=1, .ZA=0, .MI=2 , .LR=0}; 
  amns_c_reactant_type species3  = {.ZN=6, .ZA=0, .MI=12, .LR=1}; 
  amns_c_reactant_type species4  = {.ZN=1, .ZA=1, .MI=2 , .LR=1}; 
  amns_c_reaction_type xx_rx  = {.string = "CX"};
  void* amns_cx_handle;
  double rate;

  ITM_AMNS_CC_SETUP(&amns_handle, &error_stat);
  printf("error = %s: %s\n", error_stat.flag ? "true" : "false", error_stat.string);
  ITM_AMNS_CC_SETUP_REACTANTS(&reactants_handle, "", 0, 4);
  ITM_AMNS_CC_SET_REACTANT(reactants_handle, 1, &species1);
  ITM_AMNS_CC_SET_REACTANT(reactants_handle, 2, &species2);
  ITM_AMNS_CC_SET_REACTANT(reactants_handle, 3, &species3);
  ITM_AMNS_CC_SET_REACTANT(reactants_handle, 4, &species4);
  ITM_AMNS_CC_SETUP_TABLE(amns_handle, &xx_rx, reactants_handle, &amns_cx_handle, &error_stat);
  printf("error = %s: %s\n", error_stat.flag ? "true" : "false", error_stat.string);
  ITM_AMNS_CC_RX_0_B(amns_cx_handle, &rate, 100.0, 1e20, &error_stat);
  printf("error = %s: %s\n", error_stat.flag ? "true" : "false", error_stat.string);
  printf("Rate = %25.15e\n", rate);
  ITM_AMNS_CC_FINISH_TABLE(&amns_cx_handle, &error_stat);
  printf("error = %s: %s\n", error_stat.flag ? "true" : "false", error_stat.string);
  ITM_AMNS_CC_FINISH_REACTANTS(&reactants_handle);
  ITM_AMNS_CC_FINISH(&amns_handle, &error_stat);
  printf("error = %s: %s\n", error_stat.flag ? "true" : "false", error_stat.string);
  return 0;
}