Mercurial > repos > public > sbplib

comparison +multiblock/nonConformingInterfaceOptions.m @ 913:95cd70f4b07d feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix in nonConformingInterfaceOptions
author Martin Almquist <malmquist@stanford.edu>
date Sun, 25 Nov 2018 21:08:55 -0800
parents ecfccf236af0
children
comparison
equal deleted inserted replaced
912:1cdf5ead2a16 913:95cd70f4b07d
35 elseif mi/mj == 2 35 elseif mi/mj == 2
36 % Block i is finer 36 % Block i is finer
37 37
38 switch interpOpSet 38 switch interpOpSet
39 case 'MC' 39 case 'MC'
40 interpOpSet = sbp.InterpMC(mi, mj, orders{i}, orders{j}); 40 interpOpSet = sbp.InterpMC(mj+1, mi+1, orders{j}, orders{i});
41 I_i2j_good = interpOpSet.IF2C; 41 I_i2j_good = interpOpSet.IF2C;
42 I_i2j_bad = interpOpSet.IF2C; 42 I_i2j_bad = interpOpSet.IF2C;
43 I_j2i_good = interpOpSet.IC2F; 43 I_j2i_good = interpOpSet.IC2F;
44 I_j2i_bad = interpOpSet.IC2F; 44 I_j2i_bad = interpOpSet.IC2F;
45 45
46 case 'AWW' 46 case 'AWW'
47 interpOpSetF2C = sbp.InterpAWW(mi+1, mj+1, orders{i}, orders{j}, 'F2C'); 47 interpOpSetF2C = sbp.InterpAWW(mj+1, mi+1, orders{j}, orders{i}, 'F2C');
48 interpOpSetC2F = sbp.InterpAWW(mi+1, mj+1, orders{i}, orders{j}, 'C2F'); 48 interpOpSetC2F = sbp.InterpAWW(mj+1, mi+1, orders{j}, orders{i}, 'C2F');
49 I_i2j_good = interpOpSetF2C.IF2C; 49 I_i2j_good = interpOpSetF2C.IF2C;
50 I_i2j_bad = interpOpSetC2F.IF2C; 50 I_i2j_bad = interpOpSetC2F.IF2C;
51 I_j2i_good = interpOpSetC2F.IC2F; 51 I_j2i_good = interpOpSetC2F.IC2F;
52 I_j2i_bad = interpOpSetF2C.IC2F; 52 I_j2i_bad = interpOpSetF2C.IC2F;
53 end 53 end