Fusion is_polynomial in check_model

models.pl

@@ -309,6 +309,7 @@ check_model :-
   biocham_command,
   doc('Check if the current model is well formed and strict, (see TCS15_FGS).'),
   doc('If it is not the case, proposes a witness.'),
+  doc('Check if the current ODE system is polynomial.'),
   single_model(ModelId),
   findall(
     Molecule,
@@ -338,8 +339,21 @@ check_model :-
     ;
       format("The model is NOT strict.~n", [])
     )
+  ),
+  (
+    with_current_ode_system(
+      forall(
+        ode(_Var, Expr),
+        is_polynomial(Expr)
+      )
+    )
+  ->
+    format('The current ode system is polynomial.~n', [])
+  ;
+    format('The current ode system is NOT polynomial.~n', [])
   ).
 
+
 :- devdoc('\\section{Public API}').
 
 

ode.pl

@@ -23,7 +23,6 @@
     add_influences_from_ode_system/0,
     (init)/1,
     remove_fraction/0,
-    is_polynomial/0,
 %    load_ode_system/1,
     load_reactions_from_ode/1,
 %    add_ode_system/1,
@@ -454,21 +453,6 @@ remove_fraction :-
   select_ode_system(without_fraction_ode).
 
 
-is_polynomial :-
-  biocham_command,
-  doc('Check if the current ODE system is polynomial and print YES or NO accordingly.'),
-  (
-    with_current_ode_system(
-    forall(
-      ode(_Var, Expr),
-      is_polynomial(Expr)
-    )
-    )
-  ->
-    format('Yes, the current ode system is polynomial.~n', [])
-  ;
-    format('No, the current ode system is not polynomial.~n', [])
-  ).
 
 
 :- devdoc('\\section{Public API}').