Mercurial > repos > public > sbplib

view +sbp/+implementations/intOpMC_orders_2to2_ratio2to1.m @ 1118:07d0caf915e4 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Introduce optFlag so that one can choose not to build the optimization-related cell of matrices called B. It is too computationally costly and should probably be done in a different way.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 05 May 2019 19:05:31 -0700
parents c923fe6197ff
children
line wrap: on
line source

function [IC2F,IF2C] = intOpMC_orders_2to2_ratio2to1(mc)

mf = 2*(mc-1) + 1;

stencil_F2C = [1.0./4.0,1.0./2.0,1.0./4.0];
stencil_width = length(stencil_F2C);
stencil_halfwidth = (stencil_width-1)/2;

BC_F2C = [1.0./2.0,1.0./2.0];

Hc = speye(mc,mc);
Hc(1,1) = 1/2;
Hc(end,end) = 1/2;

Hf = 1/2*speye(mf,mf);
Hf(1,1) = 1/4;
Hf(end,end) = 1/4; 

IF2C = sparse(mc,mf);
[BCrows, BCcols] = size(BC_F2C);
IF2C(1:BCrows, 1:BCcols) = BC_F2C;
IF2C(mc-BCrows+1:mc, mf-BCcols+1:mf) = rot90(BC_F2C,2);

for i = BCrows+1 : mc-BCrows
	IF2C(i,(2*i-stencil_halfwidth-1) :(2*i+stencil_halfwidth-1))...
		 = stencil_F2C;
end

IC2F = Hf\(IF2C'*Hc);