Mercurial > repos > public > sbplib
view +sbp/+implementations/intOpMC_orders_6to6_ratio2to1.m @ 1347:ac54767ae1fb feature/poroelastic tip
Add interface, not fully compatible.
| author | Martin Almquist <martin.almquist@it.uu.se> |
|---|---|
| date | Tue, 30 Apr 2024 14:58:35 +0200 |
| parents | c923fe6197ff |
| children |
line wrap: on
line source
% Marks New interpolation operators % 6th order accurate (diagonal norm) % M=19 is the minimum amount of points on the coarse mesh function [IC2F,IF2C] = intOpMC_orders_6to6_ratio2to1(M_C) M_F=M_C*2-1; % Coarse to fine I1=zeros(M_F,M_C); t1= [6854313/6988288 , 401925/6988288 , -401925/6988288 ... , 133975/6988288 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; ... 560547/1537664 , 1201479/1537664 , -240439/1537664 ... , 16077/1537664 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; ... 203385/5552128 , 1225647/1388032 , 364155/2776064 ... , -80385/1388032 , 39385/5552128 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; ... -145919/2743808 , 721527/1371904 , 105687/171488 , ... -25/256 , 23631/2743808 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ; ... -178863/4033024 , 1178085/8066048 , ... 1658587/2016512 , 401925/4033024 , -15/512 , ... 43801/8066048 , 0 , 0 , 0 , 0 , 0 , 0 ; ... -1668147/11213056 , 4193225/11213056 , ... 375675/2803264 , 2009625/2803264 , ... -984625/11213056 , 3/256 , 0 , 0 , 0 , 0 , 0 , 0 ; ... -561187/2949120 , 831521/1474560 , -788801/1474560 ... , 412643/368640 , 39385/589824 , -43801/1474560 ... , 5/1024 , 0 , 0 , 0 , 0 , 0 ; ... 23/1024 , 23/147456 , -43435/221184 , ... 26795/36864 , 39385/73728 , -43801/442368 , ... 3/256 , 0 , 0 , 0 , 0 , 0 ; ... 79379/368640 , -1664707/2949120 , 284431/737280 , ... 26795/294912 , 606529/737280 , 43801/589824 , ... -15/512 , 5/1024 , 0 , 0 , 0 , 0 ; ... 3589/27648 , -2225/6144 , 22939/73728 , ... -26795/221184 , 39385/73728 , 43801/73728 , ... -25/256 , 3/256 , 0 , 0 , 0 , 0 ; ... -720623/14745600 , 10637/92160 , -89513/1474560 , ... -5359/147456 , 39385/589824 , 3372677/3686400 , ... 75/1024 , -15/512 , 5/1024 , 0 , 0 , 0 ; ... -6357/81920 , 55219/276480 , -8707/61440 , ... 5359/368640 , -39385/442368 , 43801/73728 , ... 75/128 , -25/256 , 3/256 , 0 , 0 , 0 ; ... -13315/884736 , 2589/65536 , -479/16384 , ... 5359/884736 , -7877/294912 , 43801/589824 , ... 231/256 , 75/1024 , -15/512 , 5/1024 , 0 , 0 ; ... 8299/737280 , -7043/245760 , 5473/276480 , 0 , ... 7877/737280 , -43801/442368 , 75/128 , ... 75/128 , -25/256 , 3/256 , 0 , 0 ; ... 11027/2949120 , -8461/884736 , 655/98304 , 0 , ... 7877/1769472 , -43801/1474560 , 75/1024 , ... 231/256 , 75/1024 , -15/512 , 5/1024 , 0 ; ... -601/614400 , 601/245760 , -601/368640 , 0 , 0 , ... 43801/3686400 , -25/256 , 75/128 , ... 75/128 , -25/256 , 3/256 , 0 ; ... -601/1474560 , 601/589824 , -601/884736 , 0 , 0 , ... 43801/8847360 , -15/512 , 75/1024 , ... 231/256 , 75/1024 , -15/512 , 5/1024 ] ; t2= [5/1024 , -15/512 , 75/1024, 231/256 , 75/1024 , -15/512 , 5/1024]; t3= [3/256 , -25/256 , 75/128 , 75/128 , -25/256 , 3/256]; I1(1:17,1:12)=t1; I1(M_F-16:M_F,M_C-11:M_C)=fliplr(flipud(t1)); I1(18,7:12)=t3; for i=19:2:M_F-18 j=(i-3)/2; I1(i,j-1:j+5)=t2; I1(i+1,j:j+5)=t3; end % Fine to coarse I2=zeros(M_C,M_F); t1=[6854313/13976576 , 2802735/3494144 , ... 1016925/27953152 , -729595/6988288 , ... -894315/13976576 , -1668147/6988288 , ... -8417805/27953152 , 15525/436768 , ... 1190685/3494144 , 89725/436768 , ... -2161869/27953152 , -858195/6988288 , ... -332875/13976576 , 124485/6988288 , ... 165405/27953152 , -5409/3494144 , -9015/13976576; ... 80385/12301312 , 1201479/3075328 , 1225647/6150656 ... , 721527/3075328 , 1178085/24602624 , ... 838645/6150656 , 60843/300032 , 345/6150656 , ... -4994121/24602624 , -100125/768832 , ... 31911/768832 , 55219/768832 , 349515/24602624 , ... -63387/6150656 , -42305/12301312 , 5409/6150656 ... , 9015/24602624 ; ... -80385/5552128 , -240439/1388032 , 364155/5552128 ... , 105687/173504 , 1658587/2776064 , 75135/694016 ... , -2366403/5552128 , -217175/1388032 , ... 853293/2776064 , 344085/1388032 , ... -268539/5552128 , -78363/694016 , -64665/2776064 ... , 5473/347008 , 29475/5552128 , -1803/1388032 , ... -3005/5552128]; t2=[ 5/2048 , 3/512 , -15/1024 , -25/512 , ... 75/2048 , 75/256 , 231/512 , 75/256 , ... 75/2048 , -25/512 , -15/1024 , 3/512 , 5/2048]; I2(1:3,1:17)=t1; I2(M_C-2:M_C,M_F-16:M_F)=fliplr(flipud(t1)); for i=4:M_C-3 j=2*(i-4)+1; I2(i,j:j+12)=t2; end IC2F = sparse(I1); IF2C = sparse(I2);