Changeset 3263cf5 in tol

Timestamp:

Oct 21, 2008, 11:04:48 AM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

e3f5af6

Parents:

3c7e287

Message:

New method ModSes.Mix.Constraint.Create

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@43 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/dynhlm/_db_api.tools.modses.mix.tol (modified) (1 diff)

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/dynhlm/_db_api.tools.modses.mix.tol

@@ -158 +158 @@
   Eval("ModSes.Mix.OrderRelation.Create.ok=ok")
 };
+
+//////////////////////////////////////////////////////////////////////////////
+Text _.autodoc.member.ModSes.Mix.Constraint.Create =
+"Inserts a new arbitrary constraint
+
+  leftBound <= linComb[1] + ··· + linComb[k] <= rightBound
+
+where 
+  leftBound could be -1/0 or ? if there are no lower limit
+  rightBound could be +1/0 or ? if there are no upper limit
+  Each element of linComb is a NameBlock like this
+  [[
+    Text _.id_node; 
+    Text _.id_parameter; 
+    Real _.coef
+  ]]
+
+Example: In order to insert this inequation
+
+  3.5*nodo1::beta2 -2.1*nodo2::beta5 <= 7.1*nodo1::beta3
+
+Syntax could be 
+
+  Real { ModSes.Mix.Constraint.Create(
+    my_dbName, my_model,my_session,
+    Real leftBound = -1/0,
+    Set  linComb = SetOfNameBlock
+    ([[
+      Text _.id_node = \"nodo1\"; 
+      Text _.id_parameter = \"beta2\"; 
+      Real _.coef = 3.5
+    ]],[[
+      Text _.id_node = \"nodo2\"; 
+      Text _.id_parameter = \"beta5\"; 
+      Real _.coef = -2.1
+    ]],[[
+      Text _.id_node = \"nodo1\"; 
+      Text _.id_parameter = \"beta3\"; 
+      Real _.coef = -7.1
+    ]]),
+    Real rightBound = 0)};
+";
+Real ModSes.Mix.Constraint.Create(Text dbName,
+                                  Text model,
+                                  Text session,
+                                  Real leftBound,
+                                  Set  linComb,
+                                  Real rightBound)
+//////////////////////////////////////////////////////////////////////////////
+{
+  NameBlock _table = BysMcmc::Bsr::DynHlm::DBApi::_.table;
+  Real numInequation = 
+  {
+    Real aux = BysMcmc::Bsr::DynHlm::DBApi::SqlTableTrace(
+    "SELECT max(nu_inequation) FROM "+dbName+"."+_table::ConstrBorder+"\n"
+    "WHERE id_model='"+model+"' AND\n"
+    "      id_session='"+session+"';","")[1][1];
+    If(IsUnknown(aux),1,aux+1)
+  };
+  Set nodes = Classify(linComb,Real(NameBlock a, NameBlock b)
+  {
+    Compare(a::_.id_node, b::_.id_node)
+  });
+  Text id_node = If(Card(nodes)==1, (linComb[1])::_.id_node,
+  {
+    Real If(Not(BysMcmc::Bsr::DynHlm::DBApi::SqlTableTrace(
+    "SELECT count(*) FROM "+dbName+"."+_table::Node+"\n"
+    "WHERE id_model='"+model+"' AND\n"
+    "      id_session='"+session+"' AND\n"
+    "      id_node='MIXTURE';","")[1][1]),
+    {
+      Real BysMcmc::Bsr::DynHlm::DBApi::SqlExecTrace(
+        "INSERT INTO "+_table::Node+"\n"
+        "VALUES('"+model+"','"+session+"','MIXTURE','MIX',999999999,-1,'')\n")
+    });
+    "MIXTURE"
+  });
+  Real ok.1 = BysMcmc::Bsr::DynHlm::DBApi::SqlExecTrace(
+  "INSERT INTO "+dbName+"."+_table::ConstrBorder+"\n"+
+  "VALUES('"+model+"','"+session+"',"+
+  IntText(numInequation)+","+
+  "'"+id_node+"',"+
+  BysMcmc::Bsr::DynHlm::DBApi::SqlFormatReal(leftBound)+","+
+  BysMcmc::Bsr::DynHlm::DBApi::SqlFormatReal(rightBound)+");");
+  Real ok.2 = SetMin(EvalSet(linComb, Real(NameBlock a)
+  {
+    BysMcmc::Bsr::DynHlm::DBApi::SqlExecTrace(
+    "INSERT INTO "+dbName+"."+_table::ConstrLinComb+"\n"+
+    "VALUES('"+model+"','"+session+"',"+
+    IntText(numInequation)+","+
+    "'"+a::_.id_node+"',"+
+    "'"+a::_.id_parameter+"',"+
+    BysMcmc::Bsr::DynHlm::DBApi::SqlFormatReal(a::_.coef)+");")
+  }));
+  Eval("ModSes.Mix.Constraint.Create.ok=ok.1 & ok.2")
+};
+PutDescription(_.autodoc.member.ModSes.Mix.Constraint.Create,
+               ModSes.Mix.Constraint.Create);