% seti contrast Fresnel Synth
% is used in tests/fresnelSynth.m

% -- setGrid
seti.dim = 2;
seti.rCD = 0.2;
seti.nCD = 256; % original

% -- fresnel data (loadData.m)
seti.expData = 'fresnel';
seti.fresnelFreq = 5*1E9; % standard 3 GHz (4 GHz)
seti.fresnelFile = 'inexpdata/fresnel_opus_1/dielTM_dec8f.exp'; % 1--8 GHz, test: single cylinder, 4 GHz
%seti.fresnelFile = 'inexpdata/fresnel_opus_1/twodielTM_8f.exp'; % 1--8 GHz
filename = seti.fresnelFile;

seti.nuMax = 7; seti.ampCalc = 1; % default (A1)
%seti.nuMax = 7; seti.ampCalc = 2; % test (A2)
%seti.nuMax = 7; seti.ampCalc = 3; % test (A3)

seti.nuMax = 10; seti.ampCalc = 1; % default (B1)
%seti.nuMax = 10; seti.ampCalc = 2; % test (B2)
%seti.nuMax = 10; seti.ampCalc = 3; % test (B3)

% 3 GHz: single dielectricum
% B1: 1.3 & 14.9 & 15.3
% B2:  1.8 & 15.0 & 15.7
% B3 (with 2E3): 1.7 & 15.0 & 15.7
% B3 (with 1E5): 1.6 & 14.9 & 15.7
% B3 (with 1E8): 1.5 & 14.9 & 15.5

% 5 GHz: single dielectricum
% B1:  3.4 & 20.1 & 22.4
% B2: 21.2 & 20.5 & 26.9
% B3 (with 2E3): 18.9 & 19.6 & 24.7
% B3 (with 1E5): 15.2 & 19.1 & 23.1

% single dielectricum with method (B1)
%1\,GHz & 0.9 & 29.8 & 39.7\\
%2\,GHz & 0.6 & 18.0 & 14.9\\
%3\,GHz & 1.3 & 14.9 & 15.3\\
%4\,GHz & 0.8 & 25.2 & 19.7\\
%5\,GHz & 3.4 & 20.1 & 22.4\\
%6\,GHz & 4.2 & 39.7 & 32.7\\
%7\,GHz & 16.4 & 22.2 & 28.1\\
%8\,GHz & 18.1 & 44.3 & 39.1

% two dielectrics with method (B1)
%1\,GHz & 1.0 & 11.5 & 64.5\\
%2\,GHz & 0.6 & 46.7 &  8.1\\
%3\,GHz & 1.3 & 13.4 & 21.6\\
%4\,GHz & 0.8 & 47.9 & 15.6\\
%5\,GHz & 3.3 & 20.3 & 27.4\\
%6\,GHz & 4.0 & 59.5 & 36.8\\
%7\,GHz & 16.1 & 25.2 & 42.9\\
%8\,GHz & 17.8 & 72.2 & 41.1

%-----

% -- setKernel
cLight = 2.99792458E8; % light velocity
seti.k = 2*pi*seti.fresnelFreq/cLight;
seti.model = 'helmholtz';

% -- setMeas (measurements)
% simulate fresnel_opus_1:
seti.incType = 'pointSource'; % used in op. 1 and 2
seti.measType = 'nearField'; % used in op. 1 and 2
seti.incNb = 36;
seti.measNb = 72;
seti.radSrc = 0.72;
seti.radMeas = 0.76;

% -- inversion process (regularization)
%seti.invNo = 3;
seti.invNo = 6;
seti.delta = 0; % to test

%seti = setiAuxRec(seti,usefindpdatau,usevaralpha);
seti.alpha = 1E8; % old: 2E4 corresponding 0.05 (in old code)
seti.beta = 0; % because we expect no sharp edges...
seti.physBounds = [-500,500,0,inf];
seti.nOut = 30;

% -- minimization with pda
seti.pdaN = 50;