tblrdc.c File Reference

Reduce look-up table size using recursive, error-controlled interpolation. More...

#include "jurassic.h"

Go to the source code of this file.

Functions
int	main (int argc, char *argv[])

Detailed Description

Reduce look-up table size using recursive, error-controlled interpolation.

Definition in file tblrdc.c.

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 62 of file tblrdc.c.

                {
  ctl_t ctl;
 
  /* Print usage information... */
  USAGE;
 
  /* Check arguments... */
  if (argc < 6)
    ERRMSG("Missing or invalid command-line arguments.\n\n"
           "Usage: tblrdc <ctl> <tblbase_in> <tblfmt_in> <tblbase_out> <tblfmt_out> [KEY VALUE ...]\n\n"
           "Use -h for full help.");
 
  /* Read control parameters... */
  read_ctl(argc, argv, &ctl);
  const double atol = scan_ctl(argc, argv, "ATOL", -1, "1e-9", NULL);
  const double rtol = scan_ctl(argc, argv, "RTOL", -1, "0.01", NULL);
 
  /* Read tables... */
  sprintf(ctl.tblbase, "%s", argv[2]);
  ctl.tblfmt = atoi(argv[3]);
  tbl_t *tbl = read_tbl(&ctl);
 
  /* Loop over trace gases and channels... */
  for (int id = 0; id < ctl.nd; id++)
    for (int ig = 0; ig < ctl.ng; ig++) {
 
      /* Init counters... */
      int total_input = 0;
      int total_output = 0;
 
      /* Loop over (p, T) combinations... */
      for (int ip = 0; ip < tbl->np[id][ig]; ip++) {
        for (int it = 0; it < tbl->nt[id][ig][ip]; it++) {
 
          const int num_points = tbl->nu[id][ig][ip][it];
 
          /* First entry is implicitly kept (in-place reduction)... */
          int write_idx = 1;
 
          /* Reduce Middle-Entries if more than 2 entries... */
          if (num_points > 2)
            reduction(tbl, id, ig, ip, it, 0, num_points - 1, &write_idx,
                      atol, rtol);
 
          /* Keep last entry if more than 1 entry... */
          if (num_points > 1) {
            copy_density(tbl, id, ig, ip, it, write_idx, num_points - 1);
            write_idx++;
          }
 
          /* Update number of kept points... */
          tbl->nu[id][ig][ip][it] = write_idx;
          total_input += num_points;
          total_output += write_idx;
        }
      }
 
      /* Write info... */
      LOG(1, "Reduction: %s | %.4f cm^-1 | n= %d | CR= %g", ctl.emitter[ig],
          ctl.nu[id], total_input,
          total_output > 0 ? (double) total_input / total_output : NAN);
 
    }
 
  /* Write tables... */
  sprintf(ctl.tblbase, "%s", argv[4]);
  ctl.tblfmt = atoi(argv[5]);
  write_tbl(&ctl, tbl);
 
  /* Free... */
  tbl_free(&ctl, tbl);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: