forked from ThorstenHellert/SC
-
Notifications
You must be signed in to change notification settings - Fork 0
/
SCscaleCircumference.m
66 lines (62 loc) · 1.27 KB
/
SCscaleCircumference.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
54
55
56
57
58
59
60
61
62
63
64
65
66
function RING = SCscaleCircumference(RING,circ,varargin)
% SCscaleCircumference
% ====================
%
% Name
% ----
% SCscaleCircumference - scales the circumference of a ring
%
% Synopsis
% --------
% `RING = SCscaleCircumference(RING, circ [, mode])`
%
%
% Description
% -----------
% The circumference of `RING` is changed by scaling the lengths of the drift
% spaces in the lattice.
%
%
% MODE
% ----
% `'abs'` (default)::
% The current circumference is scaled to the value `circ`.
% `'rel'`::
% The current circumference is scaled by the value `circ`.
%
%
% Return Value
% ------------
% `RING`::
% Lattice with new circumference.
% Get mode
if nargin==3
mode=varargin{1};
else
mode='abs';
end
% Length of the ring
C=findspos(RING,length(RING)+1);
% Sum of drift spaces in the ring
D = 0;
for ord=1:length(RING)
if strcmp(RING{ord}.PassMethod,'DriftPass')
D = D + RING{ord}.Length;
end
end
switch mode
case 'rel'
% Scaling factor for the drifts
Dscale = 1 - (1-circ) * C/D;
case 'abs'
Dscale = 1 - (C-circ)/D;
otherwise
error('Unsupported circumference scaling mode: ''%s''',mode)
end
% Apply scaling factor
for ord=1:length(RING)
if strcmp(RING{ord}.PassMethod,'DriftPass')
RING{ord}.Length = RING{ord}.Length * Dscale;
end
end
end