void
_tenDwiGageFilter(gageContext *ctx, gagePerVolume *pvl) {
char me[]="_tenDwiGageFilter";
gage_t *fw00, *fw11, *fw22, *dwi;
int fd;
tenDwiGageKindData *kindData;
unsigned int J, dwiNum;
fd = 2*ctx->radius;
dwi = pvl->directAnswer[tenDwiGageAll];
kindData = AIR_CAST(tenDwiGageKindData *, pvl->kind->data);
dwiNum = kindData->num;
if (!ctx->parm.k3pack) {
fprintf(stderr, "!%s: sorry, 6pack filtering not implemented\n", me);
return;
}
fw00 = ctx->fw + fd*3*gageKernel00;
fw11 = ctx->fw + fd*3*gageKernel11;
fw22 = ctx->fw + fd*3*gageKernel22;
/* HEY: these will have to be updated if there is ever any use for
gradients in DWIs: the second-to-last argument would change */
switch (fd) {
case 2:
for (J=0; J<dwiNum; J++) {
gageScl3PFilter2(pvl->iv3 + J*8, pvl->iv2 + J*4, pvl->iv1 + J*2,
fw00, fw11, fw22,
dwi + J, NULL, NULL,
pvl->needD[0], AIR_FALSE, AIR_FALSE);
}
break;
case 4:
for (J=0; J<dwiNum; J++) {
gageScl3PFilter4(pvl->iv3 + J*64, pvl->iv2 + J*16, pvl->iv1 + J*4,
fw00, fw11, fw22,
dwi + J, NULL, NULL,
pvl->needD[0], AIR_FALSE, AIR_FALSE);
}
break;
default:
for (J=0; J<dwiNum; J++) {
gageScl3PFilterN(fd, pvl->iv3 + J*fd*fd*fd,
pvl->iv2 + J*fd*fd, pvl->iv1 + J*fd,
fw00, fw11, fw22,
dwi + J, NULL, NULL,
pvl->needD[0], AIR_FALSE, AIR_FALSE);
}
break;
}
return;
}
void
_gageVecFilter(gageContext *ctx, gagePerVolume *pvl) {
char me[]="_gageVecFilter";
double *fw00, *fw11, *fw22, *vec, *jac, *hes;
int fd;
gageScl3PFilter_t *filter[5] = {NULL, gageScl3PFilter2, gageScl3PFilter4,
gageScl3PFilter6, gageScl3PFilter8};
unsigned int valIdx;
fd = 2*ctx->radius;
vec = pvl->directAnswer[gageVecVector];
jac = pvl->directAnswer[gageVecJacobian];
hes = pvl->directAnswer[gageVecHessian];
if (!ctx->parm.k3pack) {
fprintf(stderr, "!%s: sorry, 6pack filtering not implemented\n", me);
return;
}
fw00 = ctx->fw + fd*3*gageKernel00;
fw11 = ctx->fw + fd*3*gageKernel11;
fw22 = ctx->fw + fd*3*gageKernel22;
/* perform the filtering */
if (fd <= 8) {
for (valIdx=0; valIdx<3; valIdx++) {
filter[ctx->radius](ctx->shape,
pvl->iv3 + valIdx*fd*fd*fd,
pvl->iv2 + valIdx*fd*fd,
pvl->iv1 + valIdx*fd,
fw00, fw11, fw22,
vec + valIdx, jac + valIdx*3, hes + valIdx*9,
pvl->needD[0], pvl->needD[1], pvl->needD[2]);
}
} else {
for (valIdx=0; valIdx<3; valIdx++) {
gageScl3PFilterN(ctx->shape, fd,
pvl->iv3 + valIdx*fd*fd*fd,
pvl->iv2 + valIdx*fd*fd,
pvl->iv1 + valIdx*fd,
fw00, fw11, fw22,
vec + valIdx, jac + valIdx*3, hes + valIdx*9,
pvl->needD[0], pvl->needD[1], pvl->needD[2]);
}
}
return;
}
void
_tenDwiGageFilter(gageContext *ctx, gagePerVolume *pvl) {
char me[]="_tenDwiGageFilter";
double *fw00, *fw11, *fw22, *dwi;
int fd;
tenDwiGageKindData *kindData;
gageScl3PFilter_t *filter[5] = {NULL, gageScl3PFilter2, gageScl3PFilter4,
gageScl3PFilter6, gageScl3PFilter8};
unsigned int J, dwiNum;
fd = 2*ctx->radius;
dwi = pvl->directAnswer[tenDwiGageAll];
kindData = AIR_CAST(tenDwiGageKindData *, pvl->kind->data);
dwiNum = pvl->kind->valLen;
if (!ctx->parm.k3pack) {
fprintf(stderr, "!%s: sorry, 6pack filtering not implemented\n", me);
return;
}
fw00 = ctx->fw + fd*3*gageKernel00;
fw11 = ctx->fw + fd*3*gageKernel11;
fw22 = ctx->fw + fd*3*gageKernel22;
/* HEY: these will have to be updated if there is ever any use for
derivatives in DWIs: can't pass NULL pointers for gradient info,
and the needD[1] has to be passed */
if (fd <= 8) {
for (J=0; J<dwiNum; J++) {
filter[ctx->radius](ctx->shape,
pvl->iv3 + J*8, pvl->iv2 + J*4, pvl->iv1 + J*2,
fw00, fw11, fw22,
dwi + J, NULL, NULL,
pvl->needD[0], AIR_FALSE, AIR_FALSE);
}
} else {
for (J=0; J<dwiNum; J++) {
gageScl3PFilterN(ctx->shape, fd, pvl->iv3 + J*fd*fd*fd,
pvl->iv2 + J*fd*fd, pvl->iv1 + J*fd,
fw00, fw11, fw22,
dwi + J, NULL, NULL,
pvl->needD[0], AIR_FALSE, AIR_FALSE);
}
}
return;
}
void
_gageSclFilter(gageContext *ctx, gagePerVolume *pvl) {
char me[]="_gageSclFilter";
int fd;
double *fw00, *fw11, *fw22;
gageScl3PFilter_t *filter[5] = {NULL, gageScl3PFilter2, gageScl3PFilter4,
gageScl3PFilter6, gageScl3PFilter8};
fd = 2*ctx->radius;
if (!ctx->parm.k3pack) {
fprintf(stderr, "!%s: sorry, 6-pack filtering not implemented\n", me);
return;
}
fw00 = ctx->fw + fd*3*gageKernel00;
fw11 = ctx->fw + fd*3*gageKernel11;
fw22 = ctx->fw + fd*3*gageKernel22;
/* perform the filtering */
if (fd <= 8) {
filter[ctx->radius](ctx->shape, pvl->iv3, pvl->iv2, pvl->iv1,
fw00, fw11, fw22,
pvl->directAnswer[gageSclValue],
pvl->directAnswer[gageSclGradVec],
pvl->directAnswer[gageSclHessian],
pvl->needD);
} else {
gageScl3PFilterN(ctx->shape, fd,
pvl->iv3, pvl->iv2, pvl->iv1,
fw00, fw11, fw22,
pvl->directAnswer[gageSclValue],
pvl->directAnswer[gageSclGradVec],
pvl->directAnswer[gageSclHessian],
pvl->needD);
}
return;
}
