new option

NonLinearSolver/modelReduction/Tree.cpp

@@ -2138,7 +2138,7 @@ void Tree::clear()
   _root = NULL;
 };
 
-void Tree::initialize(bool rand, bool normalizedWeight)
+void Tree::initialize(bool rand, bool normalizedWeight, bool type)
 {
   // random weight for leaves
   unsigned seed = std::chrono::system_clock::now().time_since_epoch().count();
@@ -2174,8 +2174,15 @@ void Tree::initialize(bool rand, bool normalizedWeight)
     if (normalizedWeight)
     {
       double vv = (n->af->getVal(n->weight))/toltalWeightLeaves;
+      if (type)
+      {
         n->weight = n->af->getReciprocalVal(vv);      
       }
+      else
+      {
+        n->weight = vv;
+      }
+    }
     // propagate data
     double w = n->weight;
     TreeNode* parent = n->parent;

NonLinearSolver/modelReduction/Tree.h

@@ -59,7 +59,7 @@ class Tree
     void printTree() const;
     void printTreeInteraction(const std::string fname="treeInteraction.txt", bool colorMat = true, int dir=1) const;
     void clear();
-    void initialize(bool rand=true, bool normalizedWeight=true);
+    void initialize(bool rand=true, bool normalizedWeight=true, bool type = true);
     const TreeNode* getRootNode() const {return _root;};
     
     void assignMaterialLaws(int numPhases);