int
limnpu_sortMain(int argc, char **argv, char *me, hestParm *hparm) {
hestOpt *hopt = NULL;
char *err, *perr;
airArray *mop;
int pret;
limnPolyData *pld;
Nrrd *nin;
char *out;
hestOptAdd(&hopt, NULL, "input lpld", airTypeOther, 1, 1, &pld, NULL,
"input polydata filename",
NULL, NULL, limnHestPolyDataLMPD);
hestOptAdd(&hopt, NULL, "input nrrd", airTypeOther, 1, 1, &nin, NULL,
"input nrrd filename",
NULL, NULL, nrrdHestNrrd);
hestOptAdd(&hopt, NULL, "output", airTypeString, 1, 1, &out, NULL,
"output lpld filename");
mop = airMopNew();
airMopAdd(mop, hopt, (airMopper)hestOptFree, airMopAlways);
USAGE(myinfo);
PARSE();
airMopAdd(mop, hopt, (airMopper)hestParseFree, airMopAlways);
if (limnPolyDataPrimitiveSort(pld, nin)) {
airMopAdd(mop, err = biffGetDone(LIMN), airFree, airMopAlways);
fprintf(stderr, "%s: trouble:%s", me, err);
airMopError(mop);
return 1;
}
if (limnPolyDataSave(out, pld)) {
airMopAdd(mop, err = biffGetDone(LIMN), airFree, airMopAlways);
fprintf(stderr, "%s: trouble:%s", me, err);
airMopError(mop);
return 1;
}
airMopOkay(mop);
return 0;
}
int
limnpu_pselMain(int argc, char **argv, char *me, hestParm *hparm) {
hestOpt *hopt = NULL;
char *err, *perr;
airArray *mop;
int pret;
limnPolyData *pldIn, *pldOut;
unsigned int prange[2], pi;
Nrrd *nsel;
double *sel;
char *out;
size_t size[NRRD_DIM_MAX];
hestOptAdd(&hopt, "r", "range", airTypeUInt, 2, 2, prange, NULL,
"range of indices of primitives to select");
hestOptAdd(&hopt, NULL, "input", airTypeOther, 1, 1, &pldIn, NULL,
"input polydata filename",
NULL, NULL, limnHestPolyDataLMPD);
hestOptAdd(&hopt, NULL, "output", airTypeString, 1, 1, &out, NULL,
"output polydata filename");
mop = airMopNew();
airMopAdd(mop, hopt, (airMopper)hestOptFree, airMopAlways);
USAGE(myinfo);
PARSE();
airMopAdd(mop, hopt, (airMopper)hestParseFree, airMopAlways);
if (!( prange[0] <= pldIn->primNum-1 &&
prange[1] <= pldIn->primNum-1 )) {
fprintf(stderr, "%s: prange[0] %u or [1] %u outside range [0,%u]", me,
prange[0], prange[1], pldIn->primNum-1);
airMopError(mop);
return 1;
}
if (!( prange[0] <= prange[1] )) {
fprintf(stderr, "%s: need prange[0] %u <= [1] %u", me,
prange[0], prange[1]);
airMopError(mop);
return 1;
}
nsel = nrrdNew();
airMopAdd(mop, nsel, (airMopper)nrrdNuke, airMopAlways);
size[0] = pldIn->primNum;
if (nrrdMaybeAlloc_nva(nsel, nrrdTypeDouble, 1, size)) {
airMopAdd(mop, err = biffGetDone(NRRD), airFree, airMopAlways);
fprintf(stderr, "%s: trouble allocating buffer:%s", me, err);
airMopError(mop);
return 1;
}
sel = AIR_CAST(double *, nsel->data);
for (pi=prange[0]; pi<=prange[1]; pi++) {
sel[pi] = 1;
}
pldOut = limnPolyDataNew();
airMopAdd(mop, pldOut, (airMopper)limnPolyDataNix, airMopAlways);
if (limnPolyDataPrimitiveSelect(pldOut, pldIn, nsel)
|| limnPolyDataSave(out, pldOut)) {
airMopAdd(mop, err = biffGetDone(LIMN), airFree, airMopAlways);
fprintf(stderr, "%s: trouble:%s", me, err);
airMopError(mop);
return 1;
}
airMopOkay(mop);
return 0;
}
