You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
15 lines
380 B
15 lines
380 B
3 years ago
|
%--------------------------
|
||
|
|
% @Author: Jingqiao Hu
|
||
|
|
% @Date: 2021-05-18 22:01:15
|
||
|
|
% @LastEditTime: 2021-05-18 22:04:01
|
||
|
|
%--------------------------
|
||
|
|
function P = PK1_fast_nonlinear(u, l, F)
|
||
|
|
|
||
|
|
Fdet = det_3rank(F); % 1,m
|
||
|
|
JJ = log(Fdet);
|
||
|
|
|
||
|
|
Finv = inv_3rank(F); % 9,m
|
||
|
|
FinvT = Finv([1:3:9, 2:3:9, 3:3:9], :);
|
||
|
|
|
||
|
|
P = u(:)' .* (F - FinvT) + l(:)' .* JJ .* FinvT;
|
||
|
|
end
|