comparison util/ksvd utils/ompbox utils/ompcoreGabor.c @ 139:4bd6856a7128 ivand_dev

ompGabor mex version debuged and tested

   /*** status flags ***/
   
   DtX_specified = (DtX!=0);   /* indicates whether D'*x was provided */
   XtX_specified = (XtX!=0);   /* indicates whether sum(x.*x) was provided */
   
-  standardomp = (G==0);       /* batch-omp or standard omp are selected depending on availability of G */
+  standardomp = 1;//(G==0);       /* batch-omp or standard omp are selected depending on availability of G */
   batchomp = !standardomp;
   
   
   
   /*** allocate output matrix ***/

   n12  = (double*)mxMalloc(m/2*sizeof(double));
   
   vec_smult(1,D2, D1, n*m/2);
        
   for (i=0; i<m/2; i++) {
+	D1D2[i]=0;
+	n12[i]=0;
     for (j=0; j<n; j++) {
           D1D2[i] += D1[i*n+j];
     }
-    n12[i]=1/(1-D1D2[i]);
+    n12[i]=1/(1-D1D2[i]*D1D2[i]);
   }
   
   memcpy(D1, D      , n*m/2*sizeof(double));
   
   alpha = (double*)mxMalloc(m/2*sizeof(double));  /* contains D'*residual */

       for (k=0; k<m/2; k++){
         alpha[k] = (proj1[k] - D1D2[k]*proj2[k])*n12[k];
         beta[k]  = (proj2[k] - D1D2[k]*proj1[k])*n12[k];
       }
       for (k=0; k<m/2; k++){
+		  error[k]=0;
           for (j=0; j<n; j++){
-                  error[k]+= (abs(r[j])-D1[k*n+j]*alpha[k]-D2[(k+m/2)*n+j]*beta[k])*(abs(r[j])-D1[k*n+j]*alpha[k]-D2[(k+m/2)*n+j]*beta[k]);
+                  error[k]+= (abs(r[j])-D1[k*n+j]*alpha[k]-D2[k*n+j]*beta[k])*(abs(r[j])-D1[k*n+j]*alpha[k]-D2[k*n+j]*beta[k]);
           }
       }
       pos = maxabs(error, m/2);
       addproftime(&pd, MAXABS_TIME);