-
Notifications
You must be signed in to change notification settings - Fork 1
/
generate_weightFun.m
53 lines (48 loc) · 1.4 KB
/
generate_weightFun.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
%% generate_weights
% Generates the weight function .
%
% INPUT:
% weightFun : weight function
% dim : weight function
%
% OUTPUT:
% omega : weight function
function omega = generate_weightFun( weightFun, dim)
%% 1
if strcmp( weightFun, '1')
if dim == 1
omega = @(x) 1 + 0*x;
elseif dim == 2
omega = @(x,y) 1 + 0*x.*y;
elseif dim == 3
omega = @(x,y,z) 1 + 0*x.*y.*z;
else
error('Desired dimension not yet implemented!')
end
%% C2k
elseif strcmp( weightFun, 'C2k')
if dim == 1
omega = @(x) sqrt(1 - x.^2);
elseif dim == 2
omega = @(x,y) sqrt(1 - x.^2).*sqrt(1 - y.^2);
elseif dim == 3
omega = @(x,y,z) sqrt(1 - x.^2).*sqrt(1 - y.^2).*sqrt(1 - z.^2);
else
error('Desired dimension not yet implemented!')
end
%% sqrt(r)
elseif strcmp( weightFun, 'sqrt(r)')
if dim == 1
omega = @(x) sqrt( sqrt( x.^2 ) );
elseif dim == 2
omega = @(x,y) sqrt( sqrt( x.^2 + y.^2 ) );
elseif dim == 3
omega = @(x,y,z) sqrt( sqrt( x.^2 + y.^2 + z.^2 ) );
else
error('Desired dimension not yet implemented!')
end
%% else
else
error('Desired weight function not yet implemented!')
end
end