Steven's changes all implemented.

Author: AndreyAPopov

src/+otp/+kuramotosivashinsky/KuramotoSivashinskyProblem.m

@@ -7,7 +7,7 @@
     end
 
     methods
-        function soly = solution2real(soly)
+        function soly = convert2grid(soly)
 
             if isstruct(soly)
                 soly.y = real(ifft(soly.y));
@@ -31,14 +31,15 @@ function onSettingsChanged(obj)
             k = (1i*[0:(N/2 - 1), 0, (-N/2 + 1):-1].'/div);
             k2 = k.^2;
             k4 = k.^4;
+            k24 = k2 + k4;
 
-            obj.Rhs = otp.Rhs(@(t, u) otp.kuramotosivashinsky.f(t, u, k, k2, k4), ...
+            obj.Rhs = otp.Rhs(@(t, u) otp.kuramotosivashinsky.f(t, u, k, k24), ...
                 otp.Rhs.FieldNames.Jacobian, ...
-                @(t, u) otp.kuramotosivashinsky.jac(t,u, k, k2, k4), ...
+                @(t, u) otp.kuramotosivashinsky.jac(t,u, k, k24), ...
                 otp.Rhs.FieldNames.JacobianVectorProduct, ...
-                @(t, u, v) otp.kuramotosivashinsky.jvp(t, u, v, k, k2, k4), ...
+                @(t, u, v) otp.kuramotosivashinsky.jvp(t, u, v, k, k24), ...
                 otp.Rhs.FieldNames.JacobianAdjointVectorProduct, ...
-                @(t, u, v) otp.kuramotosivashinsky.javp(t, u, v, k, k2, k4));
+                @(t, u, v) otp.kuramotosivashinsky.javp(t, u, v, k, k24));
 
         end
 

src/+otp/+kuramotosivashinsky/f.m

@@ -1,7 +1,7 @@
-function ut = f(~, u, k, k2, k4)
+function ut = f(~, u, k, k24)
 
 u2 = fft(real(ifft(u)).^2);
 
-ut = - k2.*u - k4.*u -(k/2).*u2;
+ut = -k24.*u - (k/2).*u2;
 
 end

src/+otp/+kuramotosivashinsky/jac.m

@@ -1,5 +1,5 @@
-function j = jac(~, u, k, k2, k4)
+function j = jac(~, u, k, k24)
 
-j = -diag(k2) - diag(k4) - diag(k)*ifft(fft(diag(ifft(u))).').';
+j = -diag(k24) - k.*ifft(fft(diag(ifft(u))).').';
 
 end

src/+otp/+kuramotosivashinsky/javp.m

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

src/+otp/+kuramotosivashinsky/jvp.m

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