source: PSPA/madxPSPA/tools/numdiff/src/utils.c @ 430

/*
 o---------------------------------------------------------------------o
 |
 | Numdiff
 |
 | Copyright (c) 2012+ laurent.deniau@cern.ch
 | Gnu General Public License
 |
 o---------------------------------------------------------------------o
  
   Purpose:
     provides utilities
 
 o---------------------------------------------------------------------o
*/

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <assert.h>
#include <time.h>

#include "error.h"
#include "utils.h"
#include "args.h"

static const double pow10_tbl[2*99+1] = { 
  1e-99, 1e-98, 1e-97, 1e-96, 1e-95, 1e-94, 1e-93, 1e-92, 1e-91, 1e-90,
  1e-89, 1e-88, 1e-87, 1e-86, 1e-85, 1e-84, 1e-83, 1e-82, 1e-81, 1e-80,
  1e-79, 1e-78, 1e-77, 1e-76, 1e-75, 1e-74, 1e-73, 1e-72, 1e-71, 1e-70,
  1e-69, 1e-68, 1e-67, 1e-66, 1e-65, 1e-64, 1e-63, 1e-62, 1e-61, 1e-60,
  1e-59, 1e-58, 1e-57, 1e-56, 1e-55, 1e-54, 1e-53, 1e-52, 1e-51, 1e-50,
  1e-49, 1e-48, 1e-47, 1e-46, 1e-45, 1e-44, 1e-43, 1e-42, 1e-41, 1e-40,
  1e-39, 1e-38, 1e-37, 1e-36, 1e-35, 1e-34, 1e-33, 1e-32, 1e-31, 1e-30,
  1e-29, 1e-28, 1e-27, 1e-26, 1e-25, 1e-24, 1e-23, 1e-22, 1e-21, 1e-20,
  1e-19, 1e-18, 1e-17, 1e-16, 1e-15, 1e-14, 1e-13, 1e-12, 1e-11, 1e-10,
  1e-09, 1e-08, 1e-07, 1e-06, 1e-05, 1e-04, 1e-03, 1e-02, 1e-01,

  1e+00, 1e+01, 1e+02, 1e+03, 1e+04, 1e+05, 1e+06, 1e+07, 1e+08, 1e+09,
  1e+10, 1e+11, 1e+12, 1e+13, 1e+14, 1e+15, 1e+16, 1e+17, 1e+18, 1e+19,
  1e+20, 1e+21, 1e+22, 1e+23, 1e+24, 1e+25, 1e+26, 1e+27, 1e+28, 1e+29,
  1e+30, 1e+31, 1e+32, 1e+33, 1e+34, 1e+35, 1e+36, 1e+37, 1e+38, 1e+39,
  1e+40, 1e+41, 1e+42, 1e+43, 1e+44, 1e+45, 1e+46, 1e+47, 1e+48, 1e+49,
  1e+50, 1e+51, 1e+52, 1e+53, 1e+54, 1e+55, 1e+56, 1e+57, 1e+58, 1e+59,
  1e+60, 1e+61, 1e+62, 1e+63, 1e+64, 1e+65, 1e+66, 1e+67, 1e+68, 1e+69,
  1e+70, 1e+71, 1e+72, 1e+73, 1e+74, 1e+75, 1e+76, 1e+77, 1e+78, 1e+79,
  1e+80, 1e+81, 1e+82, 1e+83, 1e+84, 1e+85, 1e+86, 1e+87, 1e+88, 1e+89,
  1e+90, 1e+91, 1e+92, 1e+93, 1e+94, 1e+95, 1e+96, 1e+97, 1e+98, 1e+99,
};

const double *const pow10_table99 = &pow10_tbl[99];

FILE*
open_indexedFile(const char* str, int idx, const char *ext, int optext, int required)
{
  char buf[FILENAME_MAX+100];

  assert(str && ext);

  // copy filename
  strncpy(buf, str, sizeof buf);

  // find extension, if any
  const char *dot = strrchr(buf, '.');
  int pos = (int)strlen(buf);

  // remove extension if it matches ext
  if (dot && !strcmp(dot, ext)) {
    pos = dot-buf; 
    buf[pos] = 0;
  }

  // add formatted index
  if (idx > 0) pos += sprintf(buf+pos, option.fmt, idx);

  // copy filename into option for further reporting
  strncpy(option.indexed_filename, buf, sizeof option.indexed_filename);

  // add extension
  strncat(buf+pos, ext, sizeof buf - pos);

  // try to open
  FILE *fp = fopen(buf, "r");

  // try again upon failure if extension is optional
  if (!fp && optext) {
    buf[pos] = 0;
    fp = fopen(buf, "r");
  }

  if (!fp && required) ensure(fp, "failed to open %s", buf);

  // resize buffer for faster read
  if (fp && BUFSIZ < 65536 && setvbuf(fp, 0, _IOFBF, 65536)) {
    fclose(fp);
    error("unable to resize the stream buffer size");
  }

  // debug information
  if (fp) {
    if (optext) inform("processing %s", option.indexed_filename);
    debug("file %s open for reading", buf);
  } else
    trace("<-open_indexedFile: unable to open file %s for reading", buf);

  return fp;
}

void
accum_summary(int total, int failed, long lines, long numbers)
{
  if (!option.acc) return;

  int n;
//  double tz;
  struct tm tm;
  time_t now = time(0);
  int total_tests=0, total_passed=0, total_failed=0;
  long total_lines=0, total_numbers=0;
  double total_ndtime=0;

  FILE *fp = fopen(option.acc, "r+");
  if (fp) {
    // read time stamps
    n = fscanf(fp, " = tests summary (started at %d.%d.%d %d:%d:%d)\n",
                   &tm.tm_year, &tm.tm_mon, &tm.tm_mday,
                   &tm.tm_hour, &tm.tm_min, &tm.tm_sec);
    ensure(n == 6, "invalid summary file format %s", option.acc);
    tm.tm_year -= 1900;
    tm.tm_mon  -= 1;

    // correct for TZ shift (i.e. emulate non-standard timegm)
    struct tm tm2 = tm;
    option.dat_t0 = mktime(&tm2);
//    tm2 = *gmtime(&now);
//    tz = difftime(now, mktime(&tm2));

    // read diff time, line and number counts
    n = fscanf(fp, "   total diff time %lf s  -  total lines %ld  -  total numbers %ld\n",
                   &total_ndtime, &total_lines, &total_numbers);
    ensure(n == 3, "invalid summary file format %s (lines)", option.acc);
    // read tests counts
    n = fscanf(fp, "   total run  time %*f s  -  total files %d - PASSED %d - FAILED %d\n",
                   &total_tests, &total_passed, &total_failed);
    ensure(n == 3, "invalid summary file format %s (files)", option.acc);
    ensure(total_tests == total_passed+total_failed, "invalid summary count in file %s", option.acc);

    // reset file
    rewind(fp);
  }
  else {
    // create the file
    fp = fopen(option.acc, "w+");
    ensure(fp, "failed to create or read summary file %s", option.acc);

    // init the time stamp
    tm = *localtime(&option.dat_t0);
//    tz = 0;
  }

  double total_time = difftime(now, option.dat_t0); // + tz;
  total_ndtime  += (option.clk_t1 - option.clk_t0) / CLOCKS_PER_SEC;
  total_lines   += lines;
  total_numbers += numbers;

  fprintf(fp, " = tests summary (started at %04d.%02d.%02d %02d:%02d:%02d)\n",
          tm.tm_year+1900, tm.tm_mon+1, tm.tm_mday, tm.tm_hour, tm.tm_min, tm.tm_sec);

  fprintf(fp, "   total diff time %6.2lf s  -  total lines %6ld  -  total numbers %8ld\n",
              total_ndtime, total_lines, total_numbers);

  fprintf(fp, "   total run  time %6.0f s  -  total files %6d  -  PASSED %4d  -  FAILED %4d\n",
              total_time, total_tests+total, total_passed+(total-failed), total_failed+failed);

  fclose(fp);
}