int
main(int argc, char *argv[]) {
int A[SIZE];
int n;
/* Read in the array. */
n = read_int_array(A, SIZE);
/* Print the most frequent value. */
printf("The most frequent value is %d.\n", most_frequent(A, n));
return 0;
}
HBM *hbm_read_mat(const char *fname)
{ struct dsa _dsa, *dsa = &_dsa;
HBM *hbm = NULL;
dsa->fname = fname;
xprintf("hbm_read_mat: reading matrix from '%s'...\n",
dsa->fname);
dsa->fp = fopen(dsa->fname, "r");
if (dsa->fp == NULL)
{ xprintf("hbm_read_mat: unable to open '%s' - %s\n",
dsa->fname, strerror(errno));
goto fail;
}
dsa->seqn = 0;
hbm = xmalloc(sizeof(HBM));
memset(hbm, 0, sizeof(HBM));
/* read the first heading card */
if (read_card(dsa)) goto fail;
memcpy(hbm->title, dsa->card, 72), hbm->title[72] = '\0';
strtrim(hbm->title);
xprintf("%s\n", hbm->title);
memcpy(hbm->key, dsa->card+72, 8), hbm->key[8] = '\0';
strspx(hbm->key);
xprintf("key = %s\n", hbm->key);
/* read the second heading card */
if (read_card(dsa)) goto fail;
if (scan_int(dsa, "totcrd", 0, 14, &hbm->totcrd)) goto fail;
if (scan_int(dsa, "ptrcrd", 14, 14, &hbm->ptrcrd)) goto fail;
if (scan_int(dsa, "indcrd", 28, 14, &hbm->indcrd)) goto fail;
if (scan_int(dsa, "valcrd", 42, 14, &hbm->valcrd)) goto fail;
if (scan_int(dsa, "rhscrd", 56, 14, &hbm->rhscrd)) goto fail;
xprintf("totcrd = %d; ptrcrd = %d; indcrd = %d; valcrd = %d; rhsc"
"rd = %d\n", hbm->totcrd, hbm->ptrcrd, hbm->indcrd,
hbm->valcrd, hbm->rhscrd);
/* read the third heading card */
if (read_card(dsa)) goto fail;
memcpy(hbm->mxtype, dsa->card, 3), hbm->mxtype[3] = '\0';
if (strchr("RCP", hbm->mxtype[0]) == NULL ||
strchr("SUHZR", hbm->mxtype[1]) == NULL ||
strchr("AE", hbm->mxtype[2]) == NULL)
{ xprintf("%s:%d: matrix type '%s' not recognised\n",
dsa->fname, dsa->seqn, hbm->mxtype);
goto fail;
}
if (scan_int(dsa, "nrow", 14, 14, &hbm->nrow)) goto fail;
if (scan_int(dsa, "ncol", 28, 14, &hbm->ncol)) goto fail;
if (scan_int(dsa, "nnzero", 42, 14, &hbm->nnzero)) goto fail;
if (scan_int(dsa, "neltvl", 56, 14, &hbm->neltvl)) goto fail;
xprintf("mxtype = %s; nrow = %d; ncol = %d; nnzero = %d; neltvl ="
" %d\n", hbm->mxtype, hbm->nrow, hbm->ncol, hbm->nnzero,
hbm->neltvl);
/* read the fourth heading card */
if (read_card(dsa)) goto fail;
memcpy(hbm->ptrfmt, dsa->card, 16), hbm->ptrfmt[16] = '\0';
strspx(hbm->ptrfmt);
memcpy(hbm->indfmt, dsa->card+16, 16), hbm->indfmt[16] = '\0';
strspx(hbm->indfmt);
memcpy(hbm->valfmt, dsa->card+32, 20), hbm->valfmt[20] = '\0';
strspx(hbm->valfmt);
memcpy(hbm->rhsfmt, dsa->card+52, 20), hbm->rhsfmt[20] = '\0';
strspx(hbm->rhsfmt);
xprintf("ptrfmt = %s; indfmt = %s; valfmt = %s; rhsfmt = %s\n",
hbm->ptrfmt, hbm->indfmt, hbm->valfmt, hbm->rhsfmt);
/* read the fifth heading card (optional) */
if (hbm->rhscrd <= 0)
{ strcpy(hbm->rhstyp, "???");
hbm->nrhs = 0;
hbm->nrhsix = 0;
}
else
{ if (read_card(dsa)) goto fail;
memcpy(hbm->rhstyp, dsa->card, 3), hbm->rhstyp[3] = '\0';
if (scan_int(dsa, "nrhs", 14, 14, &hbm->nrhs)) goto fail;
if (scan_int(dsa, "nrhsix", 28, 14, &hbm->nrhsix)) goto fail;
xprintf("rhstyp = '%s'; nrhs = %d; nrhsix = %d\n",
hbm->rhstyp, hbm->nrhs, hbm->nrhsix);
}
/* read matrix structure */
hbm->colptr = xcalloc(1+hbm->ncol+1, sizeof(int));
if (read_int_array(dsa, "colptr", hbm->ptrfmt, hbm->ncol+1,
hbm->colptr)) goto fail;
hbm->rowind = xcalloc(1+hbm->nnzero, sizeof(int));
if (read_int_array(dsa, "rowind", hbm->indfmt, hbm->nnzero,
hbm->rowind)) goto fail;
/* read matrix values */
if (hbm->valcrd <= 0) goto done;
if (hbm->mxtype[2] == 'A')
{ /* assembled matrix */
hbm->values = xcalloc(1+hbm->nnzero, sizeof(double));
if (read_real_array(dsa, "values", hbm->valfmt, hbm->nnzero,
hbm->values)) goto fail;
}
else
{ /* elemental (unassembled) matrix */
hbm->values = xcalloc(1+hbm->neltvl, sizeof(double));
if (read_real_array(dsa, "values", hbm->valfmt, hbm->neltvl,
hbm->values)) goto fail;
}
/* read right-hand sides */
if (hbm->nrhs <= 0) goto done;
if (hbm->rhstyp[0] == 'F')
{ /* dense format */
hbm->nrhsvl = hbm->nrow * hbm->nrhs;
hbm->rhsval = xcalloc(1+hbm->nrhsvl, sizeof(double));
if (read_real_array(dsa, "rhsval", hbm->rhsfmt, hbm->nrhsvl,
hbm->rhsval)) goto fail;
}
else if (hbm->rhstyp[0] == 'M' && hbm->mxtype[2] == 'A')
{ /* sparse format */
/* read pointers */
hbm->rhsptr = xcalloc(1+hbm->nrhs+1, sizeof(int));
if (read_int_array(dsa, "rhsptr", hbm->ptrfmt, hbm->nrhs+1,
hbm->rhsptr)) goto fail;
/* read sparsity pattern */
hbm->rhsind = xcalloc(1+hbm->nrhsix, sizeof(int));
if (read_int_array(dsa, "rhsind", hbm->indfmt, hbm->nrhsix,
hbm->rhsind)) goto fail;
/* read values */
hbm->rhsval = xcalloc(1+hbm->nrhsix, sizeof(double));
if (read_real_array(dsa, "rhsval", hbm->rhsfmt, hbm->nrhsix,
hbm->rhsval)) goto fail;
}
else if (hbm->rhstyp[0] == 'M' && hbm->mxtype[2] == 'E')
{ /* elemental format */
hbm->rhsval = xcalloc(1+hbm->nrhsvl, sizeof(double));
if (read_real_array(dsa, "rhsval", hbm->rhsfmt, hbm->nrhsvl,
hbm->rhsval)) goto fail;
}
else
{ xprintf("%s:%d: right-hand side type '%c' not recognised\n",
dsa->fname, dsa->seqn, hbm->rhstyp[0]);
goto fail;
}
/* read starting guesses */
if (hbm->rhstyp[1] == 'G')
{ hbm->nguess = hbm->nrow * hbm->nrhs;
hbm->sguess = xcalloc(1+hbm->nguess, sizeof(double));
if (read_real_array(dsa, "sguess", hbm->rhsfmt, hbm->nguess,
hbm->sguess)) goto fail;
}
/* read solution vectors */
if (hbm->rhstyp[2] == 'X')
{ hbm->nexact = hbm->nrow * hbm->nrhs;
hbm->xexact = xcalloc(1+hbm->nexact, sizeof(double));
if (read_real_array(dsa, "xexact", hbm->rhsfmt, hbm->nexact,
hbm->xexact)) goto fail;
}
done: /* reading has been completed */
xprintf("hbm_read_mat: %d cards were read\n", dsa->seqn);
fclose(dsa->fp);
return hbm;
fail: /* something wrong in Danish kingdom */
if (hbm != NULL)
{ if (hbm->colptr != NULL) xfree(hbm->colptr);
if (hbm->rowind != NULL) xfree(hbm->rowind);
if (hbm->rhsptr != NULL) xfree(hbm->rhsptr);
if (hbm->rhsind != NULL) xfree(hbm->rhsind);
if (hbm->values != NULL) xfree(hbm->values);
if (hbm->rhsval != NULL) xfree(hbm->rhsval);
if (hbm->sguess != NULL) xfree(hbm->sguess);
if (hbm->xexact != NULL) xfree(hbm->xexact);
xfree(hbm);
}
if (dsa->fp != NULL) fclose(dsa->fp);
return NULL;
}
