Adaptive Gaussian Filtering SVN

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

#include "agf_defs.h"

using namespace std;
using namespace libagf;

int main(int argc, char **argv) {
  FILE *fs;  

  real_a *con;
  cls_ta *cls;
  nel_ta n1, n2;

  char *confile;
  char *clsfile;

  nel_ta n;
  dim_ta nvar;

  real_a *r;

  if (argc < 2) {
    printf("Converts classes and confidence ratings into a \n");
    printf("a difference in conditional probabilities:\n");
    printf("     R = 2*C*(c - 1);   (C = |R|)\n\n");
    printf("usage: C2R basename outfile\n\n");
    printf("    where:\n");
    printf("basename = base name of binary files:\n");
    printf("           .cls for classes, .con for confidence ratings\n");
    printf("outfile  = name of output file\n");
    exit(-1);
  }

  //generate the input file names:
  clsfile=new char[strlen(argv[1])+5];
  strcpy(clsfile, argv[1]);
  strcat(clsfile, ".cls");

  confile=new char[strlen(argv[1])+5];
  strcpy(confile, argv[1]);
  strcat(confile, ".con");

  //printf("%s; %s\n", confile, clsfile);

  //keep it as a stand-alone utility (no libraries to include):
  fs=fopen(clsfile, "r");
  fseek(fs, 0, SEEK_END);
  n1 = ftell(fs)/sizeof(cls_ta);
  rewind(fs);
  cls=new cls_ta[n1];
  fread(cls, sizeof(cls_ta), n1, fs);
  fclose(fs);

  fs=fopen(confile, "r");
  fseek(fs, 0, SEEK_END);
  n2 = ftell(fs)/sizeof(real_a);
  rewind(fs);
  con=new real_a[n2];
  fread(con, sizeof(real_a), n2, fs);
  fclose(fs);

  assert(n1 == n2);

  r=new real_a[n1];
  for (nel_ta i=0; i<n1; i++) {
    if (cls[i] == 0) {
      r[i]=con[i];
    } else {
      r[i]=-con[i];
    }
    //printf("%f\n", r[i]);
  }

  //printf("ii\n");
  fs=fopen(argv[2], "w");
  fwrite(r, sizeof(real_a), n1, fs);
  fclose(fs);

  //printf("iii\n");
  delete [] con;
  delete [] cls;
  delete [] r;

}