You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
80 lines
3.3 KiB
80 lines
3.3 KiB
11 months ago
|
function Bs=LocalBs(x,y,z,x0,y0,z0,h)
|
||
|
|
Bs=zeros(6,24,size(x,1),size(x,3),size(x,3));
|
||
|
|
Bs(1,1,:,:,:)=-((h - y + y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(1,4,:,:,:)=((h - y + y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(1,7,:,:,:)=-((y - y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(1,10,:,:,:)=((y - y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(1,13,:,:,:)=-((z - z0).*(h - y + y0))./h^3;
|
||
|
|
Bs(1,16,:,:,:)=((z - z0).*(h - y + y0))./h^3;
|
||
|
|
Bs(1,19,:,:,:)=-((y - y0).*(z - z0))./h^3;
|
||
|
|
Bs(1,22,:,:,:)=((y - y0).*(z - z0))./h^3;
|
||
|
|
|
||
|
|
Bs(2,2,:,:,:)=-((h - x + x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(2,5,:,:,:)=-((x - x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(2,8,:,:,:)=((h - x + x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(2,11,:,:,:)=((x - x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(2,14,:,:,:)=-((z - z0).*(h - x + x0))./h^3;
|
||
|
|
Bs(2,17,:,:,:)=-((x - x0).*(z - z0))./h^3;
|
||
|
|
Bs(2,20,:,:,:)=((z - z0).*(h - x + x0))./h^3;
|
||
|
|
Bs(2,23,:,:,:)=((x - x0).*(z - z0))./h^3;
|
||
|
|
|
||
|
|
Bs(3,3,:,:,:)=-((h - x + x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(3,6,:,:,:)=-((x - x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(3,9,:,:,:)=-((y - y0).*(h - x + x0))./h^3;
|
||
|
|
Bs(3,12,:,:,:)=-((x - x0).*(y - y0))./h^3;
|
||
|
|
Bs(3,15,:,:,:)=((h - x + x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(3,18,:,:,:)=((x - x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(3,21,:,:,:)=((y - y0).*(h - x + x0))./h^3;
|
||
|
|
Bs(3,24,:,:,:)=((x - x0).*(y - y0))./h^3;
|
||
|
|
|
||
|
|
Bs(4,1,:,:,:)=-((h - x + x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(4,2,:,:,:)=-((h - y + y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(4,4,:,:,:)=-((x - x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(4,5,:,:,:)=((h - y + y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(4,7,:,:,:)=((h - x + x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(4,8,:,:,:)=-((y - y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(4,10,:,:,:)=((x - x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(4,11,:,:,:)=((y - y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(4,13,:,:,:)=-((z - z0).*(h - x + x0))./h^3;
|
||
|
|
Bs(4,14,:,:,:)=-((z - z0).*(h - y + y0))./h^3;
|
||
|
|
Bs(4,16,:,:,:)=-((x - x0).*(z - z0))./h^3;
|
||
|
|
Bs(4,17,:,:,:)=((z - z0).*(h - y + y0))./h^3;
|
||
|
|
Bs(4,19,:,:,:)=((z - z0).*(h - x + x0))./h^3;
|
||
|
|
Bs(4,20,:,:,:)=-((y - y0).*(z - z0))./h^3;
|
||
|
|
Bs(4,22,:,:,:)=((x - x0).*(z - z0))./h^3;
|
||
|
|
Bs(4,23,:,:,:)=((y - y0).*(z - z0))./h^3;
|
||
|
|
|
||
|
|
Bs(5,2,:,:,:)=-((h - x + x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(5,3,:,:,:)=-((h - x + x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(5,5,:,:,:)=-((x - x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(5,6,:,:,:)=-((x - x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(5,8,:,:,:)=-((y - y0).*(h - x + x0))./h^3;
|
||
|
|
Bs(5,9,:,:,:)=((h - x + x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(5,11,:,:,:)=-((x - x0).*(y - y0))./h^3;
|
||
|
|
Bs(5,12,:,:,:)=((x - x0).*(h - z + z0))./h^3;
|
||
|
|
Bs(5,14,:,:,:)=((h - x + x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(5,15,:,:,:)=-((z - z0).*(h - x + x0))./h^3;
|
||
|
|
Bs(5,17,:,:,:)=((x - x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(5,18,:,:,:)=-((x - x0).*(z - z0))./h^3;
|
||
|
|
Bs(5,20,:,:,:)=((y - y0).*(h - x + x0))./h^3;
|
||
|
|
Bs(5,21,:,:,:)=((z - z0).*(h - x + x0))./h^3;
|
||
|
|
Bs(5,23,:,:,:)=((x - x0).*(y - y0))./h^3;
|
||
|
|
Bs(5,24,:,:,:)=((x - x0).*(z - z0))./h^3;
|
||
|
|
|
||
|
|
Bs(6,1,:,:,:)=-((h - x + x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(6,3,:,:,:)=-((h - y + y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(6,4,:,:,:)=-((x - x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(6,6,:,:,:)=((h - y + y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(6,7,:,:,:)=-((y - y0).*(h - x + x0))./h^3;
|
||
|
|
Bs(6,9,:,:,:)=-((y - y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(6,10,:,:,:)=-((x - x0).*(y - y0))./h^3;
|
||
|
|
Bs(6,12,:,:,:)=((y - y0).*(h - z + z0))./h^3;
|
||
|
|
Bs(6,13,:,:,:)=((h - x + x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(6,15,:,:,:)=-((z - z0).*(h - y + y0))./h^3;
|
||
|
|
Bs(6,16,:,:,:)=((x - x0).*(h - y + y0))./h^3;
|
||
|
|
Bs(6,18,:,:,:)=((z - z0).*(h - y + y0))./h^3;
|
||
|
|
Bs(6,19,:,:,:)=((y - y0).*(h - x + x0))./h^3;
|
||
|
|
Bs(6,21,:,:,:)=-((y - y0).*(z - z0))./h^3;
|
||
|
|
Bs(6,22,:,:,:)=((x - x0).*(y - y0))./h^3;
|
||
|
|
Bs(6,24,:,:,:)=((y - y0).*(z - z0))./h^3;
|
||
|
|
end
|