Jacobian -> PartialDerivative for time and parameters

Steven-Roberts · Steven-Roberts · commit d6b4e901ecf7 · 2021-05-12T21:45:13.000-04:00
diff --git a/src/+otp/+lorenz63/Lorenz63Problem.m b/src/+otp/+lorenz63/Lorenz63Problem.m
@@ -26,7 +26,7 @@ function onSettingsChanged(obj)
                 otp.Rhs.FieldNames.Jacobian, @(t, y) otp.lorenz63.jac(t, y, sigma, rho, beta), ...
                 otp.Rhs.FieldNames.JacobianVectorProduct, @(t, y, v) otp.lorenz63.jvp(t, y, v, sigma, rho, beta), ...
                 otp.Rhs.FieldNames.JacobianAdjointVectorProduct, @(t, y, v) otp.lorenz63.javp(t, y, v, sigma, rho, beta), ...
-                otp.Rhs.FieldNames.JacobianParameters, @(t, y) otp.lorenz63.jacparams(t, y, sigma, rho, beta), ...
+                otp.Rhs.FieldNames.PartialDerivativeParameters, @(t, y) otp.lorenz63.pdparams(t, y, sigma, rho, beta), ...
                 otp.Rhs.FieldNames.HessianVectorProduct, @(t, y, u, v) otp.lorenz63.hvp(t, y, u, v, sigma, rho, beta), ...
                 otp.Rhs.FieldNames.HessianAdjointVectorProduct, @(t, y, u, v) otp.lorenz63.havp(t, y, u, v, sigma, rho, beta));
         end
diff --git a/src/+otp/+lorenz63/jacparams.m b/src/+otp/+lorenz63/jacparams.m
diff --git a/src/+otp/+lorenz63/pdparams.m b/src/+otp/+lorenz63/pdparams.m
@@ -0,0 +1,7 @@
+function Jp = pdparams(~, y, ~, ~, ~)
+
+% This will be the jacobian with respect to sigma, rho, and beta
+
+Jp = diag([y(2) - y(1), y(1), -y(3)]);
+
+end
diff --git a/src/+otp/+qg/QuasiGeostrophicProblem.m b/src/+otp/+qg/QuasiGeostrophicProblem.m
@@ -150,9 +150,6 @@ function onSettingsChanged(obj)
                 ...
                 otp.Rhs.FieldNames.JacobianAdjointVectorProduct, @(t, psi, u) ...
                 otp.qg.javp(psi, u, L, RdnL, PdnL, Ddx, Ddy, Re, Ro), ...
-                ...
-                otp.Rhs.FieldNames.JacobianTime, @(t, psi) ...
-                otp.qg.jact(psi, L, RdnL, PdnL, Ddx, Ddy, ymat, Re, Ro) ...
                 );
             
 
diff --git a/src/+otp/+qg/jact.m b/src/+otp/+qg/jact.m
diff --git a/src/+otp/+quadratic/QuadraticProblem.m b/src/+otp/+quadratic/QuadraticProblem.m
@@ -18,9 +18,9 @@ function onSettingsChanged(obj)
             
             obj.Rhs = otp.Rhs(@(t, x) otp.quadratic.f(t, x, a, B, C), ...
                 otp.Rhs.FieldNames.Jacobian, @(t, x) otp.quadratic.jac(t, x, a, B, C), ...
-                'JacobianParametera',  @(t, x) otp.quadratic.jacpara(t, x, a, B, C), ...
-                'JacobianParameterB',  @(t, x) otp.quadratic.jacparaB(t, x, a, B, C), ...
-                'JacobianParameterC',  @(t, x) otp.quadratic.jacparC(t, x, a, B, C));
+                'PartialDerivativea',  @(t, x) otp.quadratic.pda(t, x, a, B, C), ...
+                'PartialDerivativeB',  @(t, x) otp.quadratic.pdb(t, x, a, B, C), ...
+                'PartialDerivativeC',  @(t, x) otp.quadratic.pdc(t, x, a, B, C));
             
             obj.RhsClosure = otp.Rhs(@(t, x, Btilde) otp.quadratic.fclosureB(t, x, a, B, C, Btilde));
         end
diff --git a/src/+otp/+quadratic/jacpara.m b/src/+otp/+quadratic/jacpara.m
diff --git a/src/+otp/+quadratic/pda.m b/src/+otp/+quadratic/pda.m
@@ -0,0 +1,5 @@
+function ja = pda(~, ~, a, ~, ~)
+
+ja = speye(length(a));
+
+end
diff --git a/src/+otp/+quadratic/pdb.m b/src/+otp/+quadratic/pdb.m
@@ -1,4 +1,4 @@
-function jB = jacparB(~, y, ~, B, ~)
+function jB = pdb(~, y, ~, B, ~)
 
 r = size(B, 1);
 
diff --git a/src/+otp/+quadratic/pdc.m b/src/+otp/+quadratic/pdc.m
@@ -1,4 +1,4 @@
-function jC = jacparC(~, y, ~, ~, C)
+function jC = pdc(~, y, ~, ~, C)
 
 r = size(C, 1);
 
diff --git a/src/+otp/Rhs.m b/src/+otp/Rhs.m
@@ -18,8 +18,8 @@
             ... % other fields
             'JacobianVectorProduct',        'JacobianVectorProduct', ...
             'JacobianAdjointVectorProduct', 'JacobianAdjointVectorProduct', ...
-            'JacobianParameters',           'JacobianParameters', ...
-            'JacobianTime',                 'JacobianTime', ...
+            'PartialDerivativeParameters',  'PartialDerivativeParameters', ...
+            'PartialDerivativeTime',        'PartialDerivativeTime', ...
             'HessianVectorProduct',         'HessianVectorProduct', ...
             'HessianAdjointVectorProduct',  'HessianAdjointVectorProduct', ...
             'OnEvent',                      'OnEvent');

-Original file line number
+Diff line change
@@ @@ -0,0 +1,5 @@ @@
 +function ja = pda(~, ~, a, ~, ~)
++
 +ja = speye(length(a));
++
 +end
Original file line number	Diff line number	Diff line change
`@@ -1,4 +1,4 @@`
`1`		`-function jB = jacparB(~, y, ~, B, ~)`
	`1`	`+function jB = pdb(~, y, ~, B, ~)`
`2`	`2`
`3`	`3`	`r = size(B, 1);`
`4`	`4`
Original file line number	Diff line number	Diff line change
`@@ -1,4 +1,4 @@`
`1`		`-function jC = jacparC(~, y, ~, ~, C)`
	`1`	`+function jC = pdc(~, y, ~, ~, C)`
`2`	`2`
`3`	`3`	`r = size(C, 1);`
`4`	`4`