Commit 3a89e6d9 authored by guillaume le guludec's avatar guillaume le guludec

Add four species cycle by means of sigmoidal absence indicator

parent 82d2337d
......@@ -13,7 +13,8 @@
test_auto/0,
test_hill2/0,
test_sigma/0
test_sigma/0,
test_cycle/1
]).
......@@ -295,6 +296,28 @@ sigma_absence(X, G, Y, Z, P) :-
% DEBUG - TEST ------------------------------------------------------------
test_cycle(C) :-
biocham_command,
sigma_absence(x1, g1, y1, z1, p),
sigma_absence(x2, g2, y2, z2, p),
sigma_absence(x3, g3, y3, z3, p),
sigma_absence(x4, g4, y4, z4, p),
present([g1], 0.5), present([y1], 2), present([z1], 0.5),
present([g2], 0.5), present([y2], 2), present([z2], 0.5),
present([g3], 0.5), present([y3], 2), present([z3], 0.5),
present([g4], 0.5), present([y4], 2), present([z4], 0.5),
present([p], 5),
present([x1], C),
S = [x1=[z4]=>x2,
x2=[z1]=>x3,
x3=[z2]=>x4,
x4=[z3]=>x1],
compile(S).
test_sigma :-
biocham_command,
X = 'x', G = 'g', Y = 'y', Z = 'z', P = 'p',
......@@ -302,8 +325,7 @@ test_sigma :-
present([G], 0.5),
present([Y], 2),
present([Z], 0.5),
present([P], 5),
add_reaction(0.1*X for X=>_).
present([P], 5).
test_compile(String) :-
biocham_command,
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment