int main(int argc, char *argv[])
{
int i,j,p=0;
char *filename = default_filename;
SMat sm = (SMat)malloc(sizeof(struct smat));
if (argc == 2) {
filename = argv[1];
}
sm = svdLoadSparseMatrix(filename, SVD_F_ST);
printf("col: %ld\n", sm->cols);
printf("row: %ld\n", sm->rows);
printf("val: %ld\n", sm->vals);
for (i=0;i<sm->cols;++i) {
for (j=0;j<sm->rows;++j) {
if (sm->rowind[p] == j) {
printf("%.4lf ", sm->value[p++]);
}
else {
printf(" 0 ");
}
}
printf("\n");
}
free(sm);
return 0;
}
int main(int argc, char *argv[]) {
extern char *optarg;
extern int optind;
int opt;
SVDRec R = NULL;
SMat A = NULL;
char transpose = FALSE;
int readFormat = SVD_F_ST;
int writeFormat = SVD_F_DT;
int algorithm = LAS2;
int iterations = 0;
int dimensions = 0;
char *vectorFile = NULL;
double las2end[2] = {-1.0e-30, 1.0e-30};
double kappa = 1e-6;
double exetime;
while ((opt = getopt(argc, argv, "a:c:d:e:hk:i:o:r:tv:w:")) != -1) {
switch (opt) {
case 'a':
if (!strcasecmp(optarg, "las2"))
algorithm = LAS2;
else fatalError("unknown algorithm: %s", optarg);
break;
case 'c':
if (optind != argc - 1) printUsage(argv[0]);
if (SVDVerbosity > 0) printf("Converting %s to %s\n", optarg, argv[optind]);
if (SVD_IS_SPARSE(readFormat) && SVD_IS_SPARSE(writeFormat)) {
SMat S = svdLoadSparseMatrix(optarg, readFormat);
if (!S) fatalError("failed to read sparse matrix");
if (transpose) {
if (SVDVerbosity > 0) printf(" Transposing the matrix...\n");
S = svdTransposeS(S); /* This leaks memory. */
}
svdWriteSparseMatrix(S, argv[optind], writeFormat);
} else {
DMat D = svdLoadDenseMatrix(optarg, readFormat);
if (!D) fatalError("failed to read dense matrix");
if (transpose) {
if (SVDVerbosity > 0) printf(" Transposing the matrix...\n");
D = svdTransposeD(D); /* This leaks memory. */
}
svdWriteDenseMatrix(D, argv[optind], writeFormat);
}
exit(0);
break;
case 'd':
dimensions = atoi(optarg);
if (dimensions < 0) fatalError("dimensions must be non-negative");
break;
case 'e':
las2end[1] = atof(optarg);
las2end[0] = -las2end[1];
break;
case 'h':
printUsage(argv[0]);
break;
case 'k':
kappa = atof(optarg);
break;
case 'i':
iterations = atoi(optarg);
break;
case 'o':
vectorFile = optarg;
break;
case 'r':
if (!strcasecmp(optarg, "sth")) {
readFormat = SVD_F_STH;
} else if (!strcasecmp(optarg, "st")) {
readFormat = SVD_F_ST;
} else if (!strcasecmp(optarg, "dt")) {
readFormat = SVD_F_DT;
} else if (!strcasecmp(optarg, "sb")) {
readFormat = SVD_F_SB;
} else if (!strcasecmp(optarg, "db")) {
readFormat = SVD_F_DB;
} else fatalError("bad file format: %s", optarg);
break;
case 't':
transpose = TRUE;
break;
case 'v':
SVDVerbosity = atoi(optarg);
/*if (SVDVerbosity) printf("Verbosity = %ld\n", SVDVerbosity);*/
break;
case 'w':
if (!strcasecmp(optarg, "sth")) {
writeFormat = SVD_F_STH;
} else if (!strcasecmp(optarg, "st")) {
writeFormat = SVD_F_ST;
} else if (!strcasecmp(optarg, "dt")) {
writeFormat = SVD_F_DT;
} else if (!strcasecmp(optarg, "sb")) {
writeFormat = SVD_F_SB;
} else if (!strcasecmp(optarg, "db")) {
writeFormat = SVD_F_DB;
} else fatalError("bad file format: %s", optarg);
break;
default: printUsage(argv[0]);
}
}
if (optind != argc - 1) printUsage(argv[0]);
if (SVDVerbosity > 0) printf("Loading the matrix...\n");
A = svdLoadSparseMatrix(argv[optind], readFormat);
if (!A) fatalError("failed to read sparse matrix. Did you specify the correct file type with the -r argument?");
if (transpose) {
if (SVDVerbosity > 0) printf(" Transposing the matrix...\n");
SMat T = A;
A = svdTransposeS(A);
svdFreeSMat(T);
}
if (dimensions <= 0) dimensions = imin(A->rows, A->cols);
exetime = timer();
if (SVDVerbosity > 0) printf("Computing the SVD...\n");
if (algorithm == LAS2) {
if (!(R = svdLAS2(A, dimensions, iterations, las2end, kappa)))
fatalError("error in svdLAS2");
} else {
fatalError("unknown algorithm");
}
exetime = timer() - exetime;
if (SVDVerbosity > 0) {
printf("\nELAPSED CPU TIME = %6g sec.\n", exetime);
printf("MULTIPLICATIONS BY A = %6ld\n",
(SVDCount[SVD_MXV] - R->d) / 2 + R->d);
printf("MULTIPLICATIONS BY A^T = %6ld\n",
(SVDCount[SVD_MXV] - R->d) / 2);
}
if (vectorFile) {
char filename[128];
sprintf(filename, "%s-Ut", vectorFile);
svdWriteDenseMatrix(R->Ut, filename, writeFormat);
sprintf(filename, "%s-S", vectorFile);
svdWriteDenseArray(R->S, R->d, filename, FALSE);
sprintf(filename, "%s-Vt", vectorFile);
svdWriteDenseMatrix(R->Vt, filename, writeFormat);
}
return 0;
}
