void test1()
{
matrix *mx, *u, *v;
dvector *w;
NewMatrix(&mx, 4, 5);
MatrixSet(mx, 0.f);
mx->data[0][0] = 1; mx->data[0][4] = 2;
mx->data[1][2] = 3;
mx->data[3][1] = 4;
PrintMatrix(mx);
initMatrix(&u);
initMatrix(&v);
initDVector(&w);
svd(mx, &u, &w, &v);
puts("U");
PrintMatrix(u);
puts("W");
PrintDVector(w);
puts("V");
PrintMatrix(v);
DelMatrix(&u);
DelMatrix(&v);
DelDVector(&w);
DelMatrix(&mx);
}
void LogDVector(double *inVector,
int inRL,
int inRH,
char *inMessage)
{
int i;
FILE *theLogFilePtr;
if ((theLogFilePtr = fopen(gLogFileName, "at")) == NULL)
{
fprintf(stderr, "LogDVector -- Cannot open output file.\n");
return;
}
fprintf(theLogFilePtr,"%s\n",inMessage);
fprintf(theLogFilePtr,"\n");
for(i=inRL;i<=inRH;i++)
{
fprintf(theLogFilePtr,"%1.9f\n",inVector[i]);
}
fprintf(theLogFilePtr,"\n");
fclose(theLogFilePtr);
printf("%s\n",inMessage);
PrintDVector(inVector,inRL,inRH);
printf("LogDVector: -- wrote to file %s\n",gLogFileName);
}
void PrintLDAModel(LDAMODEL* m)
{
puts("Eigenvalues"); PrintDVector(m->eval);
puts("Eigenvectors"); PrintMatrix(m->evect);
puts("Features");
PrintArray(m->features);
puts("Multivariate Normal Distribution");
PrintArray(m->mnpdf);
puts("Class Average. Each row represent a class");
PrintMatrix(m->mu);
puts("Validation...");
puts("Senstivity"); PrintDVector(m->sens);
puts("Specificity"); PrintDVector(m->spec);
puts("Positive Predicted Value"); PrintDVector(m->ppv);
puts("Negative Predicted Value"); PrintDVector(m->npv);
puts("Accuracy"); PrintDVector(m->acc);
}
