optcrit.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
/*****************************************************************************/
#include "tpcextensions.h"
/*****************************************************************************/
#include "tpcmodels.h"
/*****************************************************************************/
 
/*****************************************************************************/
typedef struct TPC_OPTCRIT {
  char *code;
  optimality_criterion id;
  char *desc;
} TPC_OPTCRIT;
/*****************************************************************************/
static TPC_OPTCRIT tpc_optcrit[]={
  {"unknown", OPTCRIT_UNKNOWN,  "Unknown optimality criterion"},
  {"OLS",     OPTCRIT_OLS,      "Ordinary Least Squares"},
  {"LMS",     OPTCRIT_LMS,      "Least Median of Squares"},
  {"LAD",     OPTCRIT_LAD,      "Least Absolute Deviations"},
  {"MAD",     OPTCRIT_MAD,      "Median Absolute Deviation"},
  {"ODR",     OPTCRIT_ODR,      "Orthogonal deviations"},
  {"",        0,                ""}
};
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int optcritNr()
{
  unsigned int i=0;
  while(tpc_optcrit[i].code[0]!='\0') i++;
  return(i);
}
/*****************************************************************************/
 
/*****************************************************************************/
char *optcritCode(
  optimality_criterion id
) {
  for(unsigned int i=0; i<optcritNr(); i++) if(id==tpc_optcrit[i].id) return(tpc_optcrit[i].code);
  return((char*)NULL);
}
/*****************************************************************************/
 
/*****************************************************************************/
char *optcritDesc(
  optimality_criterion id
) {
  for(unsigned int i=0; i<optcritNr(); i++) if(id==tpc_optcrit[i].id) return(tpc_optcrit[i].desc);
  return((char*)NULL);
}
/*****************************************************************************/
 
/*****************************************************************************/
optimality_criterion optcritId(
  const char *s
) {
  if(s==NULL || *s=='\0') return(OPTCRIT_UNKNOWN);
  for(unsigned int i=0; i<optcritNr(); i++) 
    if(!strcasecmp(tpc_optcrit[i].code, s)) return(tpc_optcrit[i].id);
  /* Try also some alternative names */
  if(!strcasecmp("SS", s)) return(OPTCRIT_OLS);
  if(!strcasecmp("WSS", s)) return(OPTCRIT_OLS);
  if(!strcasecmp("SUMSQR", s)) return(OPTCRIT_OLS);
  if(!strcasecmp("LAE", s)) return(OPTCRIT_LAD);
  if(!strcasecmp("LAV", s)) return(OPTCRIT_LAD);
  if(!strcasecmp("LAR", s)) return(OPTCRIT_LAD);
  if(!strcasecmp("Deming", s)) return(OPTCRIT_ODR);
  return(OPTCRIT_UNKNOWN);
}
/*****************************************************************************/
 
/*****************************************************************************/