-
Notifications
You must be signed in to change notification settings - Fork 0
/
tapas_ehgf_jget_transp.m
37 lines (32 loc) · 1.68 KB
/
tapas_ehgf_jget_transp.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
function [pvec, pstruct] = tapas_ehgf_jget_transp(r, ptrans)
% --------------------------------------------------------------------------------------------------
% Copyright (C) 2013-2020 Christoph Mathys, TNU, UZH & ETHZ
%
% This file is part of the HGF toolbox, which is released under the terms of the GNU General Public
% Licence (GPL), version 3. You can redistribute it and/or modify it under the terms of the GPL
% (either version 3 or, at your option, any later version). For further details, see the file
% COPYING or <http://www.gnu.org/licenses/>.
pvec = NaN(1,length(ptrans));
pstruct = struct;
l = r.c_prc.n_levels;
pvec(1:l) = ptrans(1:l); % mux_0
pstruct.mux_0 = pvec(1:l);
pvec(l+1:2*l) = exp(ptrans(l+1:2*l)); % sax_0
pstruct.sax_0 = pvec(l+1:2*l);
pvec(2*l+1:3*l) = ptrans(2*l+1:3*l); % mua_0
pstruct.mua_0 = pvec(2*l+1:3*l);
pvec(3*l+1:4*l) = exp(ptrans(3*l+1:4*l)); % saa_0
pstruct.saa_0 = pvec(3*l+1:4*l);
pvec(4*l+1) = exp(ptrans(4*l+1)); % kau
pstruct.kau = pvec(4*l+1);
pvec(4*l+2:5*l) = exp(ptrans(4*l+2:5*l)); % kax
pstruct.kax = pvec(4*l+2:5*l);
pvec(5*l+1:6*l-1) = exp(ptrans(5*l+1:6*l-1)); % kaa
pstruct.kaa = pvec(5*l+1:6*l-1);
pvec(6*l) = ptrans(6*l); % omu
pstruct.omu = pvec(6*l);
pvec(6*l+1:7*l) = ptrans(6*l+1:7*l); % omx
pstruct.omx = pvec(6*l+1:7*l);
pvec(7*l+1:8*l) = ptrans(7*l+1:8*l); % oma
pstruct.oma = pvec(7*l+1:8*l);
end