3 % product of two sinus waves in both coordinate directions
5 % computes a product of axis-dependent waves of values between
'c_init_min'
8 % required fields of params:
9 % c_init_max : maximum value in field
10 % c_init_min : minimum value in field
11 % c_init_freq_x : frequency in x-direction
12 % c_init_freq_y : frequency in y-direction
13 % c_init_phase_x : phase shift
14 % c_init_phase_y : phase shift
18 if ~isempty(glob) && size(glob,1) < size(glob,2)
19 warning('coordinates in variable glob are given row-wise, but expected them to be column-wise');
26 error('wrong number of parameters!');
29 decomp_mode = params.decomp_mode;
31 % if ismember('c_init_freq_x',params.mu_names) | ...
32 % ismember('c_init_freq_y',params.mu_names) | ...
33 % ismember('c_init_min',params.mu_names) | ...
34 % ismember('c_init_max',params.mu_names) | ...
35 % ismember('c_init_phase_x',params.mu_names) | ...
36 % ismember('c_init_phase_y',params.mu_names)
37 % error('affine decomp with respect to mu_names not possible!');
46 Uinit = params.c_init_min + ...
47 (params.c_init_max-params.c_init_min) * ...
48 0.5 * (sin_sym(params.c_init_freq_x * X(:) + params.c_init_phase_x) ...
49 .* sin_sym(params.c_init_freq_y * Y(:) + params.c_init_phase_y) ...
53 elseif decomp_mode == 1
function U0 = init_values_waveproduct(glob, params)
product of two sinus waves in both coordinate directions