Mercurial > repos > public > sbplib
comparison +parametrization/old/min_h.m @ 886:8894e9c49e40 feature/timesteppers
Merge with default for latest changes
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Thu, 15 Nov 2018 16:36:21 -0800 |
| parents | 3a3cf386bb7e |
| children |
comparison
equal
deleted
inserted
replaced
| 816:b5e5b195da1e | 886:8894e9c49e40 |
|---|---|
| 1 function [d_min, i1_min, j1_min, i2_min, j2_min] = min_h(X,Y) | |
| 2 ni = size(X,1); | |
| 3 nj = size(X,2); | |
| 4 d_min = norm([X(1,1);Y(1,1)] - [X(ni,nj);Y(ni,nj)]); | |
| 5 | |
| 6 i1_min = 0; | |
| 7 j1_min = 0; | |
| 8 i2_min = 0; | |
| 9 j2_min = 0; | |
| 10 | |
| 11 D = {[-1,-1],[0,-1],[1,-1],[1,0],[1,1],[0,1],[-1,1],[-1,0]}; | |
| 12 % D = {[0,-1],[1,0],[0,1],[-1,0]}; | |
| 13 | |
| 14 for i = 1:ni | |
| 15 for j = 1:nj | |
| 16 p1 = [X(i,j); Y(i,j)]; | |
| 17 for k = 1:length(D) | |
| 18 i2 = i+D{k}(1); | |
| 19 j2 = j+D{k}(2); | |
| 20 if i2 >= 1 && i2 <= ni && j2 >= 1 && j2 <= nj | |
| 21 p2 = [X(i2,j2); Y(i2,j2)]; | |
| 22 d = norm(p2-p1); | |
| 23 if d < d_min; | |
| 24 d_min = d; | |
| 25 i1_min = i; | |
| 26 j1_min = j; | |
| 27 i2_min = i2; | |
| 28 j2_min = j2; | |
| 29 end | |
| 30 end | |
| 31 end | |
| 32 end | |
| 33 end | |
| 34 end |