diff +sbp/InterpOpsMC.m @ 952:27ce3f653aa7 feature/utux2D

Improve comments in InterpOpsXX
author Martin Almquist <malmquist@stanford.edu>
date Fri, 07 Dec 2018 11:10:31 -0800
parents ed8c98c4d479
children
line wrap: on
line diff
--- a/+sbp/InterpOpsMC.m	Thu Dec 06 10:32:02 2018 +0100
+++ b/+sbp/InterpOpsMC.m	Fri Dec 07 11:10:31 2018 -0800
@@ -1,3 +1,19 @@
+% Interpolation operators by Mattsson and Carpenter (MC), see
+% Mattsson and Carpenter,
+% "Stable and Accurate Interpolatino Operators for High-Order Multiblock Finite DIfference Methods",
+% https://epubs.siam.org/doi/pdf/10.1137/090750068
+%
+% Let ^* denote the adjoint. These operators satsify
+%
+% Iuv2 = Iv2u^*
+%
+% Both Iu2v and Iv2u have p:th order accuracy, if the interior stencil is
+% of order 2p.
+%
+% This approach leads to a reduction of the convergence rate by one order for
+% PDEs with 2nd derivatives in space, as compared to conforming interfaces.
+% To obtain full convergence rate, use the order-preserving (OP) operators in
+% InterpOpsOP.m
 classdef InterpOpsMC < sbp.InterpOps
     properties