72 for(fi=0; fi<image->
dimt; fi++) {
74 if(image->
end[fi]>0.0) {
80 if(mode==0) lambda=-lambda;
81 if(fi==image->
dimt-1 && image->
end[fi]<=0.0)
return(3);
82 cf=hlLambda2factor_float(lambda, image->
start[fi], dur);
84 if(
IMG_TEST) printf(
"applied_dc_factor[%d] := %g\n", fi+1, cf);
92 for(pi=0; pi<image->
dimz; pi++)
93 for(i=0; i<image->
dimy; i++)
94 for(j=0; j<image->
dimx; j++)
95 image->
m[pi][i][j][fi]*=cf;
111 return(hlIsotopeCode(hlIsotopeFromHalflife(img->
isotopeHalflife/60.0)));
128 float lambda, cf, dur;
135 if(mode!=0 && image->
end[image->
dimt-1]<=0.0)
return(3);
138 for(fi=0; fi<image->
dimt; fi++) {
142 dur=image->
end[fi]-image->
start[fi];
143 if(image->
end[fi]>0.0) {
145 cf=hlLambda2factor_float(lambda, image->
start[fi], dur);
#define IMG_STATUS_OCCUPIED
char * imgIsotope(IMG *img)
int imgSetDecayCorrFactors(IMG *image, int mode)
int imgDecayCorrection(IMG *image, int mode)