Changeset 32725 for trunk/psModules/src/objects/pmPCM_MinimizeChisq.c

Timestamp:

Nov 20, 2011, 3:55:38 PM (14 years ago)

Author:

eugene

Message:

create a pre-alloc version of image smooth to avoid excessive allocs

File:

• trunk/psModules/src/objects/pmPCM_MinimizeChisq.c (modified) (3 diffs)

trunk/psModules/src/objects/pmPCM_MinimizeChisq.c

@@ -324 +324 @@
         // * threading takes place above this level
         pcm->modelConvFlux = psImageCopy (pcm->modelConvFlux, pcm->modelFlux, pcm->modelFlux->type.type);
-        psImageSmooth (pcm->modelConvFlux, pcm->sigma, pcm->nsigma);
+        psImageSmooth_PreAlloc_F32 (pcm->modelConvFlux, pcm->smdata);
+        // psImageSmooth (pcm->modelConvFlux, pcm->sigma, pcm->nsigma);
     } else {
         psImageConvolveKernel (pcm->modelConvFlux, pcm->modelFlux, NULL, 0, pcm->psfFFT);
@@ -343 +344 @@
             // * threading takes place above this level
             dmodelConv = psImageCopy (dmodelConv, dmodel, dmodel->type.type);
-            psImageSmooth (dmodelConv, pcm->sigma, pcm->nsigma);
+            psImageSmooth_PreAlloc_F32 (dmodelConv, pcm->smdata);
+            // psImageSmooth (dmodelConv, pcm->sigma, pcm->nsigma);
         } else {
             psImageConvolveKernel (dmodelConv, dmodel, NULL, 0, pcm->psfFFT);
@@ -363 +365 @@
             // * threading takes place above this level
             dmodelConv = psImageCopy (dmodelConv, dmodel, dmodel->type.type);
-            psImageSmooth (dmodelConv, pcm->sigma, pcm->nsigma);
+            psImageSmooth_PreAlloc_F32 (dmodelConv, pcm->smdata);
+            // psImageSmooth (dmodelConv, pcm->sigma, pcm->nsigma);
         } else {
             psImageConvolveFFT (dmodelConv, dmodel, NULL, 0, pcm->psf);