/** @file reaction.c

 *  @brief Header file to define reactions.

 */

#ifndef _REACTION_H_

#include "react_table.h"

typedef struct reaction_t reaction_t;

/* Maximum number of reactants that can enter into a reaction. */

#define REACTION_MAX_IN 4

#define REACTION_MAX_OUT 6

struct reaction_t {

  /* Photoionization is an special reaction where all the other parameters

   * are ignored. Be careful to put photoionization in the correct position

   * of the reaction list.

   */

  int is_photo;

  /* Number of species in and out. */

  int nin, nout;

  int input[REACTION_MAX_IN];

  int output[REACTION_MAX_OUT];

  const char *inname[REACTION_MAX_IN];

  const char *outname[REACTION_MAX_OUT];

  void (*f) (double *in, int nin, double *out, int nout, double k, double dens, react_table *rt);

  /* For reactions where k(E) is given by a reaction table. Such reactions

   * require 'f' to be f_react_table. 'k' is ignored for these reactions.

   * 'rt' is initialized with 'NULL', while 'tablefile' contains the filename

   * with the reaction table.

   *

   * This file is read in react_add().

   */

  react_table *rt;

  char *tablefile;

  double k;

  /* This allows us to define many simultaneous reactions. */

  reaction_t *next;

};

extern reaction_t *reactions_list;

extern reaction_t *reaction_index[];

#define _REACTION_H_

#endif