THMMY.gr

Error using parallel.gpu.CUDAKernel
An unexpected error occurred during CUDA execution. The CUDA error was:
CUDA_ERROR_UNKNOWN

Error in sample_kernel (line 24)
  k = parallel.gpu.CUDAKernel( '../cuda/sampleAddKernel.ptx', ...
 

patchCube = @(X,w) ...
      permute( ...
          reshape( ...
              im2col( ...
                  padarray( ...
                      X, ...
                      (w-1)./2, 'symmetric'), ...
                  w, 'sliding' ), ...
              [prod(w) size(X)] ), ...
          [2 3 1] );

  % gaussian patch
  H = fspecial('gaussian',patchSize, patchSigma);
  H = H(:) ./ max(H(:));
  

1)
..
  % create 3D cube
  B = patchCube(I, patchSize); %ανοίγει τον πίνακα ανάλογα με το patchSize πρακτίκα εάν έχεις μια I 3x3 και %patch size 3x3 κάνεις ένα padding γύρω από κάθε pixel άρα κάθε z της εικόνας 3x3x9 είναι το i στοιχείο με το %padding του γύρω γύρω. 
 [m, n, d] = size( B );
  B = reshape(B, [ m*n d ] );
% κάνει τον m x n x z πίνακα σε  m*n x z δηλαδή μια γραμμή είναι το i pixel μαζί με το padd του 


2)
% φτίαχνει την gaussian σε 3x3 (patchSize)
  % gaussian patch
  H = fspecial('gaussian',patchSize, patchSigma);
  H = H(:) ./ max(H(:));
% την ανοίγει σε 9x1  και την πολλαπλασιάζει με κάθε i pixel μαζί με το pad του (εδώ εάν κάνεις αθροισμά έχεις δυσδιάστατη συνελιξή )
  % apply gaussian patch on 3D cube
  B = bsxfun( @times, B, H' );

μέσος όρος διαφοράς
....
B[icenter*rows+jcenter]+=abs(J[(icenter+i)*rows+(jcenter+j)]-J[(xcenter+i)*rows+(ycenter+j)]);
...
B[icenter*rows+jcenter]=exp(-((B[icenter*rows+jcenter])/(z))/(filtSigma));



αποσταση πινακα γειτονιας απο μια αλλη
...
B[icenter*rows+jcenter]+=(J[(icenter+i)*rows+(jcenter+j)]-J[(xcenter+i)*rows+(ycenter+j)])*(J[(icenter+i)*rows+(jcenter+j)]-J[(xcenter+i)*rows+(ycenter+j)])
...
B[icenter*rows+jcenter]=exp(-sqrt(B[icenter*rows+jcenter])/(filtSigma));

μεσος ορος καθε περιοχης μειων μεσο ορο της αλλης περιοχης
...
mi+=J[(icenter+i)*rows+(jcenter+j)];
mx+=J[(xcenter+i)*rows+(ycenter+j)];
...
mi=mi/z;
mx=mx/z;
...
B[icenter*rows+jcenter]=exp(-(abs(mi-mx)*abs(mi-mx))/(filtSigma));

% gaussian patch
  H = fspecial('gaussian',patchSize, patchSigma);
  H = H(:) ./ max(H(:));

((J[(icenter+ki)*rows+(jcenter+lj)]-J[(xcenter+kx)*rows+(ycenter+ly)])*(J[(icenter+ki)*rows+(jcenter+lj)]-J[(xcenter+kx)*rows+(ycenter+ly)]))

chameleon= imread('chameleon.jpg');
chameleon = rgb2gray(chameleon);
chameleon=im2double(chameleon);
save('chameleon.mat', 'chameleon');

fprintf('PSNR : %0.3f dB    filtSigma : %0.5f \n',PSNR,filtSigma);
  
  if q==0
      maxpeak=0;
  else
      if PSNR>maxpeak
          maxpeak=PSNR;
          bestfiltsigma=filtSigma;
      end
  end
  if filtSigma<0.11
      filtSigma=filtSigma+0.001;
      q=q+1;
      pipeline_non_local_means
  end

Error using gpuArray/gather
An unexpected error occurred during CUDA execution. The CUDA error was:
CUDA_ERROR_ILLEGAL_ADDRESS

Error in dispInternal>iTransferPortionDense (line 36)
data = gather( subsref( x, s ) );

Error in parallel.internal.shared.buildDisplayHelper>iFirstNNumericDisplayHelper (line 72)
        maybeTruncatedValue = transferDenseFcn( x, rangeStruct );

Error in parallel.internal.shared.buildDisplayHelper>iBuildDisplayHelper (line 33)
    dh = iFirstNNumericDisplayHelper( ...

Error in parallel.internal.shared.buildDisplayHelper (line 24)
    dh = iBuildDisplayHelper( x, transferDenseFcn, transferSparseFcn, xClassName, xName, N );

Error in dispInternal (line 13)
    dh = parallel.internal.shared.buildDisplayHelper( ...

Error in gpuArray/display (line 21)
dh = dispInternal( obj, thisClassName, objName );