int
_nrrdCheck(const Nrrd *nrrd, int checkData, int useBiff) {
static const char me[]="_nrrdCheck";
int fi;
if (!nrrd) {
biffMaybeAddf(useBiff, NRRD, "%s: got NULL pointer", me);
return 1;
}
if (checkData) {
if (!(nrrd->data)) {
biffMaybeAddf(useBiff, NRRD, "%s: nrrd %p has NULL data pointer",
me, AIR_CVOIDP(nrrd));
return 1;
}
}
for (fi=nrrdField_unknown+1; fi<nrrdField_last; fi++) {
/* yes, this will call _nrrdFieldCheckSpaceInfo() many many times */
if (_nrrdFieldCheck[fi](nrrd, AIR_TRUE)) {
biffMaybeAddf(useBiff, NRRD, "%s: trouble with %s field", me,
airEnumStr(nrrdField, fi));
return 1;
}
}
return 0;
}
/*
******** nrrdDescribe
**
** writes verbose description of nrrd to given file
*/
void
nrrdDescribe(FILE *file, const Nrrd *nrrd) {
unsigned int ai;
char stmp[AIR_STRLEN_SMALL];
if (file && nrrd) {
fprintf(file, "Nrrd at 0x%p:\n", AIR_CVOIDP(nrrd));
fprintf(file, "Data at 0x%p is %s elements of type %s.\n", nrrd->data,
airSprintSize_t(stmp, nrrdElementNumber(nrrd)),
airEnumStr(nrrdType, nrrd->type));
if (nrrdTypeBlock == nrrd->type) {
fprintf(file, "The blocks have size %s\n",
airSprintSize_t(stmp, nrrd->blockSize));
}
if (airStrlen(nrrd->content)) {
fprintf(file, "Content = \"%s\"\n", nrrd->content);
}
fprintf(file, "%d-dimensional array, with axes:\n", nrrd->dim);
for (ai=0; ai<nrrd->dim; ai++) {
if (airStrlen(nrrd->axis[ai].label)) {
fprintf(file, "%d: (\"%s\") ", ai, nrrd->axis[ai].label);
} else {
fprintf(file, "%d: ", ai);
}
fprintf(file, "%s-centered, size=%s, ",
airEnumStr(nrrdCenter, nrrd->axis[ai].center),
airSprintSize_t(stmp, nrrd->axis[ai].size));
airSinglePrintf(file, NULL, "spacing=%lg, \n", nrrd->axis[ai].spacing);
airSinglePrintf(file, NULL, "thickness=%lg, \n",
nrrd->axis[ai].thickness);
airSinglePrintf(file, NULL, " axis(Min,Max) = (%lg,",
nrrd->axis[ai].min);
airSinglePrintf(file, NULL, "%lg)\n", nrrd->axis[ai].max);
if (airStrlen(nrrd->axis[ai].units)) {
fprintf(file, "units=%s, \n", nrrd->axis[ai].units);
}
}
/*
airSinglePrintf(file, NULL, "The min, max values are %lg",
nrrd->min);
airSinglePrintf(file, NULL, ", %lg\n", nrrd->max);
*/
airSinglePrintf(file, NULL, "The old min, old max values are %lg",
nrrd->oldMin);
airSinglePrintf(file, NULL, ", %lg\n", nrrd->oldMax);
/* fprintf(file, "hasNonExist = %d\n", nrrd->hasNonExist); */
if (nrrd->cmtArr->len) {
fprintf(file, "Comments:\n");
for (ai=0; ai<nrrd->cmtArr->len; ai++) {
fprintf(file, "%s\n", nrrd->cmt[ai]);
}
}
fprintf(file, "\n");
}
}
int
_nrrdFormatVTK_fitsInto(const Nrrd *nrrd, const NrrdEncoding *encoding,
int useBiff) {
static const char me[]="_nrrdFormatVTK_fitsInto";
if (!( nrrd && encoding )) {
biffMaybeAddf(useBiff, NRRD, "%s: got NULL nrrd (%p) or encoding (%p)",
me, AIR_CVOIDP(nrrd), AIR_CVOIDP(encoding));
return AIR_FALSE;
}
if (!( nrrdEncodingRaw == encoding || nrrdEncodingAscii == encoding)) {
biffMaybeAddf(useBiff, NRRD, "%s: encoding can only be %s or %s", me,
nrrdEncodingRaw->name, nrrdEncodingAscii->name);
return AIR_FALSE;
}
if (!( nrrdTypeUChar == nrrd->type
|| nrrdTypeChar == nrrd->type
|| nrrdTypeUShort == nrrd->type
|| nrrdTypeShort == nrrd->type
|| nrrdTypeUInt == nrrd->type
|| nrrdTypeInt == nrrd->type
|| nrrdTypeFloat == nrrd->type
|| nrrdTypeDouble == nrrd->type )) {
biffMaybeAddf(useBiff, NRRD,
"%s: type %s doesn't fit in VTK (as currently implemented)",
me, airEnumStr(nrrdType, nrrd->type));
return AIR_FALSE;
}
if (!( 3 == nrrd->dim
|| (4 == nrrd->dim && 3 == nrrd->axis[0].size)
|| (4 == nrrd->dim && 9 == nrrd->axis[0].size) )) {
biffMaybeAddf(useBiff, NRRD, "%s: nrrd didn't look like a volume of "
"scalars, vectors, or matrices", me);
return AIR_FALSE;
}
return AIR_TRUE;
}
/*
******** biffMsgStrlen
**
** returns length of string (not including null termination, as usual)
** of the error message that will be generated by biffMsgStrSet
*/
unsigned int
biffMsgStrlen(const biffMsg *msg) {
static const char me[]="biffMsgStrlen";
unsigned int ii, len;
if (biffMsgNoop == msg) {
return 0;
}
if (!( msg )) {
fprintf(stderr, "%s: PANIC got NULL msg %p\n", me, AIR_CVOIDP(msg));
return 0; /* exit(1); */
}
len = 0;
for (ii=0; ii<msg->errNum; ii++) {
len += AIR_UINT(strlen(msg->key)
+ strlen(msg->err[ii]) + strlen("[] \n"));
}
return len+1;
}
/*
** adds a given message to the given entry. The message is processed to
** convert all whitespace into ' ', and to eliminate whitespace at the
** end of the message.
*/
void
biffMsgAdd(biffMsg *msg, const char *err) {
static const char me[]="biffMsgAdd";
unsigned int idx;
if (biffMsgNoop == msg) {
return;
}
if (!( msg && err )) {
fprintf(stderr, "%s: PANIC got NULL msg (%p) or err (%p)\n", me,
AIR_VOIDP(msg), AIR_CVOIDP(err));
/* exit(1); */
}
idx = airArrayLenIncr(msg->errArr, 1);
if (!msg->err) {
fprintf(stderr, "%s: PANIC: couldn't add message to %s\n", me, msg->key);
/* exit(1); */
}
if (!( msg->err[idx] = airOneLinify(airStrdup(err)) )) {
fprintf(stderr, "%s: PANIC: couldn't alloc message to %s\n", me, msg->key);
/* exit(1); */
}
return;
}
