Added Jacobian Adjoint (conjugate as well) vector product.

Author: AndreyAPopov

src/+otp/+kuramotosivashinsky/KuramotoSivashinskyProblem.m

@@ -36,7 +36,9 @@ function onSettingsChanged(obj)
                 otp.Rhs.FieldNames.Jacobian, ...
                 @(t, u) otp.kuramotosivashinsky.jac(t,u, k, k2, k4), ...
                 otp.Rhs.FieldNames.JacobianVectorProduct, ...
-                @(t, u, v) otp.kuramotosivashinsky.jvp(t,u, k, k2, k4, v));
+                @(t, u, v) otp.kuramotosivashinsky.jvp(t,u, k, k2, k4, v), ...
+                otp.Rhs.FieldNames.JacobianAdjointVectorProduct, ...
+                @(t, u, v) otp.kuramotosivashinsky.javp(t,u, k, k2, k4, v));
 
         end
 

src/+otp/+kuramotosivashinsky/javp.m

@@ -0,0 +1,5 @@
+function jv = javp(~, u, k, k2, k4, v)
+
+jv = -k2.*v - k4.*v - fft(real(ifft(u)).*ifft(conj(k).*v));
+
+end

src/+otp/+kuramotosivashinsky/jvp.m

@@ -1,5 +1,5 @@
 function jv = jvp(~, u, k, k2, k4, v)
 
-jv = -k2.*v - k4.*v - k.*fft(real(ifft(u)).*real(ifft(v)));
+jv = -k2.*v - k4.*v - k.*fft(real(ifft(u)).*ifft(v));
 
 end