clusttool.cpp File Reference

#include "clusttool.h"
#include "const.h"
#include "danerror.h"
#include "emalloc.h"
#include "scanutils.h"
#include <stdio.h>
#include <math.h>

Go to the source code of this file.

Macros
#define	TOKENSIZE   80
#define	QUOTED_TOKENSIZE   "79"
#define	MAXSAMPLESIZE   65535

Functions
uinT16	ReadSampleSize (TFile *fp)
PARAM_DESC *	ReadParamDesc (TFile *fp, uinT16 N)
PROTOTYPE *	ReadPrototype (TFile *fp, uinT16 N)
PROTOSTYLE	ReadProtoStyle (const char *shape)
FLOAT32 *	ReadNFloats (TFile *fp, uinT16 N, FLOAT32 Buffer[])
void	WriteParamDesc (FILE *File, uinT16 N, const PARAM_DESC ParamDesc[])
void	WritePrototype (FILE *File, uinT16 N, PROTOTYPE *Proto)
void	WriteNFloats (FILE *File, uinT16 N, FLOAT32 Array[])
void	WriteProtoStyle (FILE *File, PROTOSTYLE ProtoStyle)
void	WriteProtoList (FILE *File, uinT16 N, PARAM_DESC ParamDesc[], LIST ProtoList, BOOL8 WriteSigProtos, BOOL8 WriteInsigProtos)

Macro Definition Documentation

MAXSAMPLESIZE

#define MAXSAMPLESIZE   65535

Definition at line 33 of file clusttool.cpp.

QUOTED_TOKENSIZE

#define QUOTED_TOKENSIZE   "79"

Definition at line 32 of file clusttool.cpp.

TOKENSIZE

#define TOKENSIZE   80

Definition at line 31 of file clusttool.cpp.

Function Documentation

ReadNFloats()

FLOAT32* ReadNFloats	(	TFile *	fp,
		uinT16	N,
		FLOAT32	Buffer[]
	)

This routine reads N floats from the specified text file and places them into Buffer. If Buffer is NULL, a buffer is created and passed back to the caller. If EOF is encountered before any floats can be read, NULL is returned.

Parameters

File	open text file to read floats from
N	number of floats to read
Buffer	pointer to buffer to place floats into

Returns

Pointer to buffer holding floats or NULL if EOF

Note

Globals: None

Exceptions: ILLEGALFLOAT

History: 6/6/89, DSJ, Created.

Definition at line 220 of file clusttool.cpp.

                                                            {
  const int kMaxLineSize = 1024;
  char line[kMaxLineSize];
  if (fp->FGets(line, kMaxLineSize) == nullptr) {
    tprintf("Hit EOF in ReadNFloats!\n");
    return nullptr;
  }
  bool needs_free = false;

  if (Buffer == NULL) {
    Buffer = static_cast<FLOAT32*>(Emalloc(N * sizeof(FLOAT32)));
    needs_free = true;
  }

  char *startptr = line;
  for (int i = 0; i < N; i++) {
    char *endptr;
    Buffer[i] = strtof(startptr, &endptr);
    if (endptr == startptr) {
      tprintf("Read of %d floats failed!\n", N);
      if (needs_free) Efree(Buffer);
      return nullptr;
    }
    startptr = endptr;
  }
  return Buffer;
}

ReadParamDesc()

PARAM_DESC* ReadParamDesc	(	TFile *	fp,
		uinT16	N
	)

This routine reads textual descriptions of sets of parameters which describe the characteristics of feature dimensions.

Exceptions:

• ILLEGALCIRCULARSPEC
• ILLEGALESSENTIALSPEC
• ILLEGALMINMAXSPEC

  Parameters

  File	open text file to read N parameter descriptions from
  N	number of parameter descriptions to read

  Returns

  Pointer to an array of parameter descriptors.

  Note

  Globals: None

  History: 6/6/89, DSJ, Created.

Definition at line 73 of file clusttool.cpp.

                                               {
  PARAM_DESC *ParamDesc;
  char linear_token[TOKENSIZE], essential_token[TOKENSIZE];

  ParamDesc = (PARAM_DESC *) Emalloc (N * sizeof (PARAM_DESC));
  for (int i = 0; i < N; i++) {
    const int kMaxLineSize = TOKENSIZE * 4;
    char line[kMaxLineSize];
    if (fp->FGets(line, kMaxLineSize) == nullptr ||
        sscanf(line, "%" QUOTED_TOKENSIZE "s %" QUOTED_TOKENSIZE "s %f %f",
               linear_token, essential_token, &ParamDesc[i].Min,
               &ParamDesc[i].Max) != 4)
      DoError(ILLEGALCIRCULARSPEC, "Illegal Parameter specification");
    if (linear_token[0] == 'c')
      ParamDesc[i].Circular = TRUE;
    else
      ParamDesc[i].Circular = FALSE;

    if (linear_token[0] == 'e')
      ParamDesc[i].NonEssential = FALSE;
    else
      ParamDesc[i].NonEssential = TRUE;
    ParamDesc[i].Range = ParamDesc[i].Max - ParamDesc[i].Min;
    ParamDesc[i].HalfRange = ParamDesc[i].Range / 2;
    ParamDesc[i].MidRange = (ParamDesc[i].Max + ParamDesc[i].Min) / 2;
  }
  return (ParamDesc);
}

ReadProtoStyle()

PROTOSTYLE ReadProtoStyle

(

const char *

shape

)

This routine reads an single token from the specified text file and interprets it as a prototype specification.

Parameters

File	open text file to read prototype style from

Returns

Prototype style read from text file

Note

Globals: None

Exceptions: ILLEGALSTYLESPEC illegal prototype style specification

History: 6/8/89, DSJ, Created.

Definition at line 191 of file clusttool.cpp.

                                             {
  switch (shape[0]) {
    case 's':
      return spherical;
    case 'e':
      return elliptical;
    case 'a':
      return automatic;
    default:
      break;
  }
  tprintf("Invalid prototype style specification:%s\n", shape);
  return elliptical;
}

ReadPrototype()

PROTOTYPE* ReadPrototype	(	TFile *	fp,
		uinT16	N
	)

This routine reads a textual description of a prototype from the specified file.

Exceptions:

• ILLEGALSIGNIFICANCESPEC
• ILLEGALSAMPLECOUNT
• ILLEGALMEANSPEC
• ILLEGALVARIANCESPEC
• ILLEGALDISTRIBUTION

  Parameters

  File	open text file to read prototype from
  N	number of dimensions used in prototype

  Returns

  List of prototypes

  Note

  Globals: None

  History: 6/6/89, DSJ, Created.

Definition at line 118 of file clusttool.cpp.

                                              {
  char sig_token[TOKENSIZE], shape_token[TOKENSIZE];
  PROTOTYPE *Proto;
  int SampleCount;
  int i;

  const int kMaxLineSize = TOKENSIZE * 4;
  char line[kMaxLineSize];
  if (fp->FGets(line, kMaxLineSize) == nullptr ||
      sscanf(line, "%" QUOTED_TOKENSIZE "s %" QUOTED_TOKENSIZE "s %d",
             sig_token, shape_token, &SampleCount) != 3) {
    tprintf("Invalid prototype: %s\n", line);
    return nullptr;
  }
  Proto = (PROTOTYPE *)Emalloc(sizeof(PROTOTYPE));
  Proto->Cluster = NULL;
  if (sig_token[0] == 's')
    Proto->Significant = TRUE;
  else
    Proto->Significant = FALSE;

  Proto->Style = ReadProtoStyle(shape_token);

  if (SampleCount < 0) DoError(ILLEGALSAMPLECOUNT, "Illegal sample count");
  Proto->NumSamples = SampleCount;

  Proto->Mean = ReadNFloats(fp, N, NULL);
  if (Proto->Mean == NULL) DoError(ILLEGALMEANSPEC, "Illegal prototype mean");

  switch (Proto->Style) {
    case spherical:
      if (ReadNFloats(fp, 1, &(Proto->Variance.Spherical)) == NULL)
        DoError(ILLEGALVARIANCESPEC, "Illegal prototype variance");
      Proto->Magnitude.Spherical =
          1.0 / sqrt((double)(2.0 * PI * Proto->Variance.Spherical));
      Proto->TotalMagnitude = pow(Proto->Magnitude.Spherical, (float)N);
      Proto->LogMagnitude = log((double)Proto->TotalMagnitude);
      Proto->Weight.Spherical = 1.0 / Proto->Variance.Spherical;
      Proto->Distrib = NULL;
      break;
    case elliptical:
      Proto->Variance.Elliptical = ReadNFloats(fp, N, NULL);
      if (Proto->Variance.Elliptical == NULL)
        DoError(ILLEGALVARIANCESPEC, "Illegal prototype variance");
      Proto->Magnitude.Elliptical = (FLOAT32 *)Emalloc(N * sizeof(FLOAT32));
      Proto->Weight.Elliptical = (FLOAT32 *)Emalloc(N * sizeof(FLOAT32));
      Proto->TotalMagnitude = 1.0;
      for (i = 0; i < N; i++) {
        Proto->Magnitude.Elliptical[i] =
            1.0 / sqrt((double)(2.0 * PI * Proto->Variance.Elliptical[i]));
        Proto->Weight.Elliptical[i] = 1.0 / Proto->Variance.Elliptical[i];
        Proto->TotalMagnitude *= Proto->Magnitude.Elliptical[i];
      }
      Proto->LogMagnitude = log((double)Proto->TotalMagnitude);
      Proto->Distrib = NULL;
      break;
    default:
      Efree(Proto);
      tprintf("Invalid prototype style\n");
      return nullptr;
  }
  return Proto;
}

ReadSampleSize()

uinT16 ReadSampleSize

(

TFile *

fp

)

This routine reads a single integer from the specified file and checks to ensure that it is between 0 and MAXSAMPLESIZE.

Parameters

File	open text file to read sample size from

Returns

Sample size

Note

Globals: None

Exceptions: ILLEGALSAMPLESIZE illegal format or range

History: 6/6/89, DSJ, Created.

Definition at line 47 of file clusttool.cpp.

                                 {
  int SampleSize = 0;

  const int kMaxLineSize = 100;
  char line[kMaxLineSize];
  if (fp->FGets(line, kMaxLineSize) == nullptr ||
      sscanf(line, "%d", &SampleSize) != 1 || (SampleSize < 0) ||
      (SampleSize > MAXSAMPLESIZE))
    DoError (ILLEGALSAMPLESIZE, "Illegal sample size");
  return (SampleSize);
}

WriteNFloats()

void WriteNFloats	(	FILE *	File,
		uinT16	N,
		FLOAT32	Array[]
	)

This routine writes a text representation of N floats from an array to a file. All of the floats are placed on one line.

Parameters

File	open text file to write N floats to
N	number of floats to write
Array	array of floats to write

Returns

None

Note

Globals: None

Exceptions: None

History: 6/6/89, DSJ, Created.

Definition at line 338 of file clusttool.cpp.

                                                          {
  for (int i = 0; i < N; i++)
    fprintf(File, " %9.6f", Array[i]);
  fprintf(File, "\n");
}

WriteParamDesc()

void WriteParamDesc	(	FILE *	File,
		uinT16	N,
		const PARAM_DESC	ParamDesc[]
	)

This routine writes an array of dimension descriptors to the specified text file.

Parameters

File	open text file to write param descriptors to
N	number of param descriptors to write
ParamDesc	array of param descriptors to write

Returns

None

Note

Globals: None

Exceptions: None

History: 6/6/89, DSJ, Created.

Definition at line 259 of file clusttool.cpp.

                                                                        {
  int i;

  for (i = 0; i < N; i++) {
    if (ParamDesc[i].Circular)
      fprintf (File, "circular ");
    else
      fprintf (File, "linear   ");

    if (ParamDesc[i].NonEssential)
      fprintf (File, "non-essential ");
    else
      fprintf (File, "essential     ");

    fprintf (File, "%10.6f %10.6f\n", ParamDesc[i].Min, ParamDesc[i].Max);
  }
}

WriteProtoList()

void WriteProtoList	(	FILE *	File,
		uinT16	N,
		PARAM_DESC	ParamDesc[],
		LIST	ProtoList,
		BOOL8	WriteSigProtos,
		BOOL8	WriteInsigProtos
	)

This routine writes a textual description of each prototype in the prototype list to the specified file. It also writes a file header which includes the number of dimensions in feature space and the descriptions for each dimension.

Parameters

File	open text file to write prototypes to
N	number of dimensions in feature space
ParamDesc	descriptions for each dimension
ProtoList	list of prototypes to be written
WriteSigProtos	TRUE to write out significant prototypes
WriteInsigProtos	TRUE to write out insignificants

Note

Globals: None

Returns

None

Note

Exceptions: None

History: 6/12/89, DSJ, Created.

Definition at line 389 of file clusttool.cpp.

                                            {
  PROTOTYPE *Proto;

  /* write file header */
  fprintf(File,"%0d\n",N);
  WriteParamDesc(File,N,ParamDesc);

  /* write prototypes */
  iterate(ProtoList)
    {
      Proto = (PROTOTYPE *) first_node ( ProtoList );
      if ((Proto->Significant && WriteSigProtos) ||
          (!Proto->Significant && WriteInsigProtos))
        WritePrototype(File, N, Proto);
    }
}

WriteProtoStyle()

void WriteProtoStyle	(	FILE *	File,
		PROTOSTYLE	ProtoStyle
	)

This routine writes to the specified text file a word which represents the ProtoStyle. It does not append a carriage return to the end.

Parameters

File	open text file to write prototype style to
ProtoStyle	prototype style to write

Returns

None

Note

Globals: None

Exceptions: None

History: 6/8/89, DSJ, Created.

Definition at line 355 of file clusttool.cpp.

                                                        {
  switch (ProtoStyle) {
    case spherical:
      fprintf (File, "spherical");
      break;
    case elliptical:
      fprintf (File, "elliptical");
      break;
    case mixed:
      fprintf (File, "mixed");
      break;
    case automatic:
      fprintf (File, "automatic");
      break;
  }
}

WritePrototype()

void WritePrototype	(	FILE *	File,
		uinT16	N,
		PROTOTYPE *	Proto
	)

This routine writes a textual description of a prototype to the specified text file.

Parameters

File	open text file to write prototype to
N	number of dimensions in feature space
Proto	prototype to write out

Returns

None

Note

Globals: None

Exceptions: None

History: 6/12/89, DSJ, Created.

Definition at line 288 of file clusttool.cpp.

                                                            {
  int i;

  if (Proto->Significant)
    fprintf (File, "significant   ");
  else
    fprintf (File, "insignificant ");
  WriteProtoStyle (File, (PROTOSTYLE) Proto->Style);
  fprintf (File, "%6d\n\t", Proto->NumSamples);
  WriteNFloats (File, N, Proto->Mean);
  fprintf (File, "\t");

  switch (Proto->Style) {
    case spherical:
      WriteNFloats (File, 1, &(Proto->Variance.Spherical));
      break;
    case elliptical:
      WriteNFloats (File, N, Proto->Variance.Elliptical);
      break;
    case mixed:
      for (i = 0; i < N; i++)
      switch (Proto->Distrib[i]) {
        case normal:
          fprintf (File, " %9s", "normal");
          break;
        case uniform:
          fprintf (File, " %9s", "uniform");
          break;
        case D_random:
          fprintf (File, " %9s", "random");
          break;
        case DISTRIBUTION_COUNT:
          ASSERT_HOST(!"Distribution count not allowed!");
      }
      fprintf (File, "\n\t");
      WriteNFloats (File, N, Proto->Variance.Elliptical);
  }
}