casa6
Public
Commits
| 2430 2430 | end[2]=0; |
| 2431 2431 | Matrix<Complex> convPlane(wtconv(begin,end).reform(IPosition(2, convSize, convSize))); |
| 2432 2432 | Float maxAbsConvFunc, minAbsConvFunc; |
| 2433 2433 | IPosition minpos, maxpos; |
| 2434 2434 | minMax(minAbsConvFunc, maxAbsConvFunc, minpos, maxpos, amplitude(convPlane)); |
| 2435 2435 | Bool found=false; |
| 2436 2436 | Int trial=0; |
| 2437 2437 | Float suplimit = 5e-3; |
| 2438 2438 | if (isLarge) |
| 2439 2439 | suplimit = 1e-2; |
| 2440 - | |
| 2440 + | //cerr << "###SUPLIMIT " << suplimit << endl; |
| 2441 2441 | for (trial = convSize / 2 - 2; trial > 0; trial--) { |
| 2442 2442 | // Searching down a diagonal |
| 2443 2443 | if (abs(convPlane(convSize / 2 - trial, convSize / 2 - trial)) > |
| 2444 2444 | (suplimit * maxAbsConvFunc)) { |
| 2445 2445 | found = true; |
| 2446 2446 | trial = Int(sqrt(2.0 * Float(trial * trial))); |
| 2447 2447 | |
| 2448 2448 | break; |
| 2449 2449 | } |
| 2450 2450 | } |