Add simple test of mex_ifsh.

This commit is contained in:
Jostein R. Natvig 2010-09-22 11:48:08 +00:00
parent c81336db05
commit 297b79d4f0

9
test_mex_ifsh.m Normal file
View File

@ -0,0 +1,9 @@
G = computeGeometry(cartGrid([30,30,1]));
src = [];
src = addSource(src, 1, 1);
src = addSource(src, G.cells.num, -1);
bc = [];
rock.perm = ones(G.cells.num, 1);
rock.poro = ones(G.cells.num, 1);
x = initResSol(G, 0, 0);
x = mex_ifsh(x, G, rock, bc, src)