linny-r 2.1.4 → 2.1.6

package/static/scripts/linny-r-model.js

@@ -542,7 +542,7 @@ class LinnyRModel {
         e = this.objectByName(en);
     if(!e) return `Unknown model entity "${en}"`;
     const
-        ao = ea[1].split('@'),
+        ao = (ea.length > 1 ? ea[1].split('@') : ['']),
         a = ao[0].trim();
     // Valid if no attribute, as all entity types have a default attribute. 
     if(!a) return true;
@@ -1037,11 +1037,11 @@ class LinnyRModel {
     // NOTE: A dimension is a list of one or more relevant selectors.
     this.dimensions.length = 0;
     // NOTE: Ignore the equations dataset.
-    for(let d in this.datasets) if(this.datasets.hasOwnProperty(d) &&
-        this.datasets[d] !== this.equations_dataset) {
+    for(let k in this.datasets) if(this.datasets.hasOwnProperty(k) &&
+        this.datasets[k] !== this.equations_dataset) {
       // Get the selector list for this dataset.
       // NOTE: Ignore wildcard selectors!
-      this.processSelectorList(this.datasets[d].plainSelectors);
+      this.processSelectorList(this.datasets[k].plainSelectors);
     }
     // Analyze constraint bound lines in the same way.
     for(let k in this.constraints) if(this.constraints.hasOwnProperty(k)) {
@@ -1636,8 +1636,7 @@ class LinnyRModel {
     const ci = this.indexOfChart(title);
     if(ci >= 0) return this.charts[ci];
     // Otherwise, add it. NOTE: Unlike datasets, charts are not "entities".
-    let c = new Chart();
-    c.title = title;
+    let c = new Chart(title);
     if(node) c.initFromXML(node);
     this.charts.push(c);
     // Sort the chart titles alphabetically...
@@ -3075,10 +3074,12 @@ class LinnyRModel {
       if(n) {
         // NOTE: Use a "dummy experiment object" as parent for SA runs.
         const dummy = {title: SENSITIVITY_ANALYSIS.experiment_title};
+        let r = 0;
         for(const c of n.childNodes) if(c.nodeName === 'experiment-run') {
-          const xr = new ExperimentRun(dummy, i);
+          const xr = new ExperimentRun(dummy, r);
           xr.initFromXML(c);
           this.sensitivity_runs.push(xr);
+          r++;
         }
       }
       n = childNodeByTag(node, 'experiments');
@@ -3321,7 +3322,7 @@ class LinnyRModel {
     vbls.sort((a, b) => UI.compareFullNames(a.displayName, b.displayName));
     // Create a new chart as dummy, so without adding it to this model.
     const
-        c = new Chart(),
+        c = new Chart('__d_u_m_m_y__c_h_a_r_t__'),
         wcdm = [];
     for(const v of vbls) {
       // NOTE: Prevent adding wildcard dataset modifiers more than once.
@@ -3337,9 +3338,9 @@ class LinnyRModel {
     return [c.dataAsString, c.statisticsAsString];
   }
     
-  get listOfAllSelectors() {
-    // Returns list of all dataset modifier selectors as a "dictionary"
-    // like so: {selector_1: [list of datasets], ...}
+  get dictOfAllSelectors() {
+    // Returns "dictionary" of all dataset modifier selectors like so:
+    // {selector_1: [list of datasets], ...}
     const ds_dict = {};
     for(let k in this.datasets) if(this.datasets.hasOwnProperty(k)) {
       const ds = this.datasets[k];
@@ -6240,19 +6241,19 @@ class Cluster extends NodeBox {
   }
   
   get nestingLevel() {
-    // Return the "depth" of this cluster in the cluster hierarchy
+    // Return the "depth" of this cluster in the cluster hierarchy.
     if(this.cluster) return this.cluster.nestingLevel + 1; // recursion!
     return 0;
   }
   
   get toBeIgnored() {
-    // Return TRUE if this cluster or some parent cluster is set to be ignored
+    // Return TRUE if this cluster or some parent cluster is set to be ignored.
     return this.ignore || MODEL.ignoreClusterInThisRun(this) ||
         (this.cluster && this.cluster.toBeIgnored); // recursion!
   }
   
   get blackBoxed() {
-    // Return TRUE if this cluster or some parent cluster is marked as black box
+    // Return TRUE if this cluster or some parent cluster is marked as black box.
     return this.black_box ||
         (this.cluster && this.cluster.blackBoxed); // recursion!
   }
@@ -7928,6 +7929,7 @@ class Process extends Node {
     const a = {name: this.displayName};
     a.LB = this.lower_bound.asAttribute;
     a.UB = (this.equal_bounds ? a.LB : this.upper_bound.asAttribute);
+    if(this.grid) a.LB = -a.UB;
     a.IL = this.initial_level.asAttribute;
     a.LCF = this.pace_expression.asAttribute;
     if(MODEL.solved) {
@@ -8861,13 +8863,13 @@ class Link {
   }
 
   get identifier() {
-    // NOTE: link IDs are based on the node codes rather than IDs, as this
-    // prevents problems when nodes are renamed
+    // NOTE: Link IDs are based on the node codes rather than IDs, as this
+    // prevents problems when nodes are renamed.
     return this.from_node.code + '___' + this.to_node.code;
   }
 
   get attributes() {
-    // NOTE: link is named by its tab-separated node names
+    // NOTE: Link is named by its tab-separated node names.
     const a = {name: this.from_node.displayName + '\t' + this.to_node.displayName};
     a.R = this.relative_rate.asAttribute;
     if(MODEL.infer_cost_prices) a.SOC = this.share_of_cost;
@@ -9220,7 +9222,7 @@ class Dataset {
   }
   
   get attributes() {
-    // NOTE: modifiers are appended as additional lines of text
+    // NOTE: Modifiers are appended as additional lines of text.
     const a = {name: this.displayName};
     a.D = '\t' + (this.vector ? this.vector[MODEL.t] : this.default_value);
     for(let k in this.modifiers) if(this.modifiers.hasOwnProperty(k)) {
@@ -9755,13 +9757,14 @@ class ChartVariable {
     this.wildcard_index = false;
   }
   
-  setProperties(obj, attr, stck, clr, sf=1, lw=1, vis=true, sort='not') {
+  setProperties(obj, attr, stck, clr, sf=1, abs=false, lw=1, vis=true, sort='not') {
     // Sets the defining properties for this chart variable.
     this.object = obj;
     this.attribute = attr;
     this.stacked = stck;
     this.color = clr;
     this.scale_factor = sf;
+    this.absolute = abs;
     this.line_width = lw;
     this.visible = vis;
     this.sorted = sort;
@@ -9777,9 +9780,11 @@ class ChartVariable {
     // Returns the display name for this variable. This is the name of
     // the Linny-R entity and its attribute, followed by its scale factor
     // unless it equals 1 (no scaling).
-    const sf = (this.scale_factor === 1 ? '' :
-        // NOTE: Pass tiny = TRUE to permit very small scaling factors.
-        ` (x${VM.sig4Dig(this.scale_factor, true)})`);
+    const
+        bar = (this.absolute ? '\u2503' : ''),
+        sf = (this.scale_factor === 1 ? '' :
+            // NOTE: Pass tiny = TRUE to permit very small scaling factors.
+            ` (x${VM.sig4Dig(this.scale_factor, true)})`);
     // Display name of equation is just the equations dataset selector. 
     if(this.object instanceof DatasetModifier) {
       let eqn = this.object.selector;
@@ -9795,7 +9800,7 @@ class ChartVariable {
         // method name (leading colon replaced by the prefixer ": ").
         eqn = this.chart.prefix + UI.PREFIXER + eqn.substring(1);
       }
-      return eqn + sf;
+      return bar + eqn + bar + sf;
     }
     // NOTE: Same holds for "dummy variables" added for wildcard
     // dataset selectors.
@@ -9804,11 +9809,13 @@ class ChartVariable {
       if(this.wildcard_index !== false) {
         eqn = eqn.replace('??', this.wildcard_index);
       }
-      return eqn + sf;
+      return bar + eqn + bar + sf;
     }
     // NOTE: Do not display the vertical bar if no attribute is specified.
-    if(!this.attribute) return this.object.displayName + sf;
-    return this.object.displayName + UI.OA_SEPARATOR + this.attribute + sf;
+    if(!this.attribute) {
+      return bar + this.object.displayName + bar + sf;
+    }
+    return bar + this.object.displayName + '|' + this.attribute + bar + sf;
   }
   
   get asXML() {
@@ -9818,7 +9825,9 @@ class ChartVariable {
     if(MODEL.black_box_entities.hasOwnProperty(id)) {
       id = UI.nameToID(MODEL.black_box_entities[id]);
     }
-    const xml = ['<chart-variable', (this.stacked ? ' stacked="1"' : ''),
+    const xml = ['<chart-variable',
+        (this.stacked ? ' stacked="1"' : ''),
+        (this.absolute ? ' absolute="1"' : ''),
         (this.visible ? ' visible="1"' : ''),
         (this.wildcard_index !== false ?
             ` wildcard-index="${this.wildcard_index}"` : ''),
@@ -9876,6 +9885,7 @@ class ChartVariable {
         nodeParameterValue(node, 'stacked') === '1',
         nodeContentByTag(node, 'color'),
         safeStrToFloat(nodeContentByTag(node, 'scale-factor')),
+        nodeParameterValue(node, 'absolute') === '1',
         safeStrToFloat(nodeContentByTag(node, 'line-width')),
         nodeParameterValue(node, 'visible') === '1',
         nodeParameterValue(node, 'sorted') || 'not');
@@ -9970,7 +9980,10 @@ class ChartVariable {
         v = 0;
       }
       // Scale the value unless run result (these are already scaled!).
-      if(!rr) v *= this.scale_factor;
+      if(!rr) {
+        if(this.absolute) v = Math.abs(v);
+        v *= this.scale_factor;
+      }
       this.vector.push(v);
       // Do not include values for t = 0 in statistics.
       if(t > 0) {
@@ -10184,7 +10197,7 @@ class Chart {
     return '#c00000';
   }
 
-  addVariable(n, a) {
+  addVariable(n, a='') {
     // Add variable [entity name `n`|attribute `a`] to the chart unless
     // it is already in the variable list.
     let dn = n + UI.OA_SEPARATOR + a;
@@ -10220,7 +10233,7 @@ class Chart {
       }
     } else {
       const v = new ChartVariable(this);
-      v.setProperties(obj, a, false, this.nextAvailableDefaultColor, 1, 1);
+      v.setProperties(obj, a, false, this.nextAvailableDefaultColor);
       this.variables.push(v);
     }
     return this.variables.length - 1;
@@ -11197,9 +11210,9 @@ class ExperimentRunResult {
       this.x_variable = true;
       this.object_id = v.object.identifier;
       this.attribute = v.attribute;
-      this.was_ignored = MODEL.ignored_entities[this.object_id];
+      this.was_ignored = MODEL.ignored_entities[this.object_id] || false;
       if(this.was_ignored) {
-        // Chart variable entity was ignored => all results are undefined
+        // Chart variable entity was ignored => all results are undefined.
         this.vector = [];
         this.N = VM.UNDEFINED;
         this.sum = VM.UNDEFINED;
@@ -11236,6 +11249,11 @@ class ExperimentRunResult {
         // statistic.
         this.last = (this.vector.length > 0 ?
             this.vector[this.vector.length - 1] : VM.UNDEFINED);
+        // NOTE: For sensitivity analyses, the vector is NOT stored because
+        // the SA reports only the descriptive statistics.
+        if(this.run.experiment === SENSITIVITY_ANALYSIS.experiment) {
+          this.vector.length = 0;
+        }
       }
     } else if(v instanceof Dataset) {
       // This dataset will be an "outcome" dataset => store statistics only
@@ -11329,12 +11347,12 @@ class ExperimentRunResult {
     }
     // The vector MAY need to be scaled to model time by different methods,
     // but since this is likely to be rare, such scaling is performed
-    // "lazily", so the method-specific vectors are initially set to NULL.
+    // "lazily", so the method-specific vectors are initially empty.
     this.resetScaledVectors();
   }
 
   resetScaledVectors() {
-    // Set the special vectors to null, so they will be recalculated.
+    // Clear the special vectors, so they will be recalculated.
     this.scaled_vectors = {'NEAREST': [], 'MEAN': [], 'SUM': [], 'MAX': []};
   }
   
@@ -11612,8 +11630,12 @@ class ExperimentRun {
     // NOTE: All equations are also considered to be outcomes EXCEPT
     // methods (selectors starting with a colon).
     this.eq_list = [];
-    const eml = Object.keys(MODEL.equations_dataset.modifiers);
-    for(const em of eml) if(!em.startsWith(':')) this.eq_list.push(em);
+    // NOTE: For sensitivity analyses, equations are NOT outcomes, as all
+    // SA outcomes must be specified explicitly.
+    if(this.experiment !== SENSITIVITY_ANALYSIS.experiment) {
+      const eml = Object.keys(MODEL.equations_dataset.modifiers);
+      for(const em of eml) if(!em.startsWith(':')) this.eq_list.push(em);
+    }
     const
         cv = this.experiment.variables.length,
         oc = this.oc_list.length,
@@ -11918,7 +11940,7 @@ class Experiment {
     }
     return index;
   }
-
+  
   isDimensionSelector(s) {
     // Return TRUE if `s` is a dimension selector in this experiment.
     for(const dim of this.dimensions) if(dim.indexOf(s) >= 0) return true;
@@ -11928,6 +11950,15 @@ class Experiment {
     return false;
   }
   
+  get allDimensionSelectors() {
+    // Return list of all dimension selectors in this experiment.
+    const
+        dict = MODEL.dictOfAllSelectors,
+        dims = [];
+    for(let s in dict) if(this.isDimensionSelector(s)) dims.push(s);
+    return dims;
+  }
+
   get asXML() {
     let d = '';
     for(const dim of this.dimensions) {
@@ -12152,12 +12183,6 @@ class Experiment {
     }
   }
   
-  get allDimensionSelectors() {
-    const sl = Object.keys(MODEL.listOfAllSelectors);
-    // Add selectors of actor, iterator and settings dimensions.
-    return sl;
-  }
-
   orthogonalSelectors(c) {
     // Return TRUE iff the selectors in set `c` all are elements of
     // different experiment dimensions.
@@ -13026,7 +13051,7 @@ class Constraint {
   }
 
   get typeLetter() {
-    return 'C';
+    return 'B';
   }
 
   get identifier() {

package/static/scripts/linny-r-utils.js

@@ -189,11 +189,16 @@ function uniformDecimals(data) {
     }
     maxi = Math.max(maxi, ss[0].length);
   }
-  // STEP 2: Convert the data to a uniform format
+  // STEP 2: Convert the data to a uniform format.
+  const special = ['\u221E', '-\u221E', '\u2047', '\u00A2'];
   for(let i = 0; i < data.length; i++) {
-    const f = parseFloat(data[i]);
+    const
+        v = data[i],
+        f = parseFloat(v);
     if(isNaN(f)) {
-      data[i] = '\u26A0'; // Unicode warning sign
+      // Keep special values such as infinity, and replace error values
+      // by Unicode warning sign.
+      if(special.indexOf(v) < 0) data[i] = '\u26A0'; 
     } else if(maxe > 0) {
     // Convert ALL numbers to exponential notation with two decimals (1.23e+7)
       const v = f.toExponential(2);
@@ -405,51 +410,59 @@ function patternList(str) {
 
 function patternMatch(str, patterns) {
   // Returns TRUE when `str` matches the &|^-pattern.
-  // NOTE: If a pattern starts with equals sign = then `str` must
-  // equal the rest of the pattern to match; if it starts with a tilde
-  // ~ then `str` must start with the rest of the pattern to match.
+  // NOTE: If a pattern starts with an opening bracket [ then `str` must
+  // start with the rest of the pattern to match. If it ends with a closing
+  // bracket ] then `str` must end with the first part of the pattern.
+  // In this way, [pattern] denotes that `str` should exactly match
   for(let i = 0; i < patterns.length; i++) {
     const p = patterns[i];
     // NOTE: `p` is an OR sub-pattern that tests for a set of "plus"
     // sub-sub-patterns (all of which should match) and a set of "min"
     // sub-sub-patters (all should NOT match)
     let pm,
+        swob,
+        ewcb,
         re,
         match = true;
     for(let j = 0; match && j < p.plus.length; j++) {
       pm = p.plus[j];
-      if(pm.startsWith('=')) {
-        match = (str === pm.substring(1));
-      } else if(pm.startsWith('~')) {
+      swob = pm.startsWith('[');
+      ewcb = pm.endsWith(']');
+      if(swob && ewcb) {
+        match = (str === pm.slice(1, -1));
+      } else if(swob) {
         match = str.startsWith(pm.substring(1));
+      } else if(ewcb) {
+        match = str.endsWith(pm.slice(0, -1));
       } else {
         match = (str.indexOf(pm) >= 0);
       }
-      // If no match, check whether pattern contains wildcards 
+      // If no match, check whether pattern contains wildcards.
       if(!match && pm.indexOf('#') >= 0) {
         // If so, rematch using regular expression that tests for a
-        // number or a ?? wildcard
+        // number or a ?? wildcard.
         let res = pm.split('#');
         for(let i = 0; i < res.length; i++) {
           res[i] = escapeRegex(res[i]);
         }
         res = res.join('(\\d+|\\?\\?)');
-        if(pm.startsWith('=')) {
-          res = '^' + res + '$';
-        } else if(pm.startsWith('~')) {
-          res = '^' + res;
-        }
+        if(swob) res = '^' + res;
+        if(ewcb) res += '$';
         re = new RegExp(res, 'g');
         match = re.test(str);
       }
     }
-    // Any "min" match indicates NO match for this sub-pattern,
+    // Any "min" match indicates NO match for this sub-pattern.
     for(let j = 0; match && j < p.min.length; j++) {
       pm = p.min[j];
-      if(pm.startsWith('=')) {
-        match = (str !== pm.substring(1));
-      } else if(pm.startsWith('~')) {
+      swob = pm.startsWith('[');
+      ewcb = pm.endsWith(']');
+      if(swob && ewcb) {
+        match = (str !== pm.slice(1, -1));
+      } else if(swob) {
         match = !str.startsWith(pm.substring(1));
+      } else if(ewcb) {
+        match = !str.endsWith(pm.slice(0, -1));
       } else {
         match = (str.indexOf(pm) < 0);
       }
@@ -461,11 +474,8 @@ function patternMatch(str, patterns) {
           res[i] = escapeRegex(res[i]);
         }
         res = res.join('(\\d+|\\?\\?)');
-        if(pm.startsWith('=')) {
-          res = '^' + res + '$';
-        } else if(pm.startsWith('~')) {
-          res = '^' + res;
-        }
+        if(swob) res = '^' + res;
+        if(ewcb) res += '$';
         re = new RegExp(res, 'g');
         match = !re.test(str);
       }
@@ -973,18 +983,20 @@ function nameToLines(name, actor_name = '') {
   // the node box.
   let m = actor_name.length;
   const
-      d = Math.floor(Math.sqrt(0.3 * name.length)),
+      d = Math.floor(Math.sqrt(0.25 * name.length)),
       // Do not wrap strings shorter than 13 characters (about 50 pixels).
       limit = Math.max(Math.ceil(name.length / d), m, 13),
-      a = name.split(' ');
-  // Split words at '-' when wider than limit
+      // NOTE: Do not split on spaces followed by a number or a single
+      // capital letter.
+      a = name.split(/\s(?!\d+:|\d+$|[A-Z]\W)/);
+  // Split words at '-' when wider than limit.
   for(let j = 0; j < a.length; j++) {
     if(a[j].length > limit) {
       const sw = a[j].split('-');
       if(sw.length > 1) {
-        // Replace j-th word by last fragment of split string
+        // Replace j-th word by last fragment of split string.
         a[j] = sw.pop();
-        // Insert remaining fragments before
+        // Insert remaining fragments before.
         while(sw.length > 0) a.splice(j, 0, sw.pop() + '-');
       }
     }

package/static/scripts/linny-r-vm.js

@@ -885,7 +885,7 @@ class ExpressionParser {
             // of commas, semicolons and spaces.
             x.r = run_spec.split(/[\,\;\/\s]+/g);
             // NOTE: The VMI instruction accepts `x.r` to be a list of selectors
-            // or an integer number.
+            // or an integer number. 
           } else {
             // If the specifier does start with a #, trim it...
             run_spec = run_spec.substring(1);
@@ -933,6 +933,20 @@ class ExpressionParser {
             x.x = MODEL.experiments[n];
           }
         }
+        // If run specifier `x.r` is a list, check whether all elements in the
+        // list are selectors in a dimension of experiment `x.x` (if specified).
+        // If experiment is unknown, check against the list of all selectors
+        // defined in the model.
+        if(Array.isArray(x.r)) {
+          const
+              sl = (x.x instanceof Experiment ? x.x.allDimensionSelectors :
+                  Object.keys(MODEL.dictOfAllSelectors)),
+              unknown = complement(x.r, sl);
+          if(unknown.length) {
+            msg = pluralS(unknown.length, 'unknown selector') + ': <tt>' +
+                unknown.join(' ') + '</tt>';
+          }
+        }
         // END of code for parsing an experiment result specifier.
         // Now proceed with parsing the variable name.
         
@@ -1350,7 +1364,10 @@ class ExpressionParser {
       const
           parts = name.split(UI.PREFIXER),
           tail = parts.pop();
-      if(parts.length > 0) {
+      if(!tail && parts.length) {
+        // Prefix without its trailing colon+space could identify an entity.
+        obj = MODEL.objectByID(UI.nameToID(parts.join(UI.PREFIXER)));
+      } else if(parts.length > 0) {
         // Name contains at least one prefix => last part *could* be a
         // method name, so look it up after adding a leading colon.
         const method = MODEL.equationByID(UI.nameToID(':' + tail));
@@ -3377,13 +3394,19 @@ class VirtualMachine {
     // Infer cycle basis for combined power grids for which Kirchhoff's
     // voltage law must be enforced.
     if(MODEL.with_power_flow) {
-      MONITOR.logMessage(1, 'POWER FLOW: ' +
+      this.logMessage(1, 'POWER FLOW: ' +
           pluralS(Object.keys(MODEL.power_grids).length, 'grid'));
+      if(MODEL.ignore_grid_capacity) this.logMessage(1,
+          'NOTE: Assuming infinite grid line cacity');
+      if(MODEL.ignore_KVL) this.logMessage(1,
+          'NOTE: Disregarding Kirchhoff\'s Voltage Law');
+      if(MODEL.ignore_power_losses) this.logMessage(1,
+          'NOTE: Disregarding transmission losses');
       POWER_GRID_MANAGER.inferCycleBasis();
       if(POWER_GRID_MANAGER.messages.length > 1) {
         UI.warn('Check monitor for power grid warnings');
       }
-      MONITOR.logMessage(1, POWER_GRID_MANAGER.messages.join('\n'));
+      this.logMessage(1, POWER_GRID_MANAGER.messages.join('\n'));
       if(POWER_GRID_MANAGER.cycle_basis.length) this.logMessage(1,
           'Enforcing Kirchhoff\'s voltage law for ' +
           POWER_GRID_MANAGER.cycleBasisAsString);
@@ -6190,7 +6213,7 @@ Solver status = ${json.status}`);
     } catch(err) {
       const msg = `ERROR while processing solver data for block ${bnr}: ${err}`;
       console.log(msg);      
-      MONITOR.logMessage(bnr, msg);
+      this.logMessage(bnr, msg);
       UI.alert(msg);
       this.stopSolving();
       this.halted = true;
@@ -7269,12 +7292,12 @@ function VMI_push_run_result(x, args) {
               rr = r.results[rri],
               tsd = r.time_step_duration,
               // Get the delta-t multiplier: divide model time step duration
-              // by time step duration of the experiment run if they differ 
+              // by time step duration of the experiment run if they differ.
               dtm = (Math.abs(tsd - model_dt) < VM.NEAR_ZERO ? 1 : model_dt / tsd);
           let stat = rrspec.s;
-          // For outcome datasets without specific statistic, default to LAST
+          // For outcome datasets without specific statistic, default to LAST.
           if(!(stat || rr.x_variable)) stat = 'LAST';
-          // For a valid experiment variable, the default value is 0
+          // For a valid experiment variable, the default value is 0.
           v = 0;
           if(stat) {
             if(stat === 'LAST') {
@@ -7302,14 +7325,14 @@ function VMI_push_run_result(x, args) {
               console.log(trc.join(''));
             }
           } else {
-            // No statistic => return the vector for local time step
+            // No statistic => return the vector for local time step,
             // using here, too, the delta-time-modifier to adjust the offsets
             // for different time steps per experiment.
             const tot = twoOffsetTimeStep(x.step[x.step.length - 1],
                 args[1], args[2], args[3], args[4], dtm, x);
             // Scale the (midpoint) time step (at current model run time scale)
             // to the experiment run time scale and get the run result value.
-            // NOTE: the .m property specifies the time scaling method, and
+            // NOTE: The .m property specifies the time scaling method, and
             // the .p property whether the run result vector should be used as
             // a periodic time series.
             v = rr.valueAtModelTime(tot[0], model_dt, rrspec.m, rrspec.p);
@@ -7325,6 +7348,11 @@ function VMI_push_run_result(x, args) {
       }
     }
   }
+  // Truncate near-zero values.
+  if(v && Math.abs(v) < VM.SIG_DIF_FROM_ZERO) {
+    console.log('NOTE: Truncated experiment run result', v, 'to zero');
+    v = 0;
+  }
   x.push(v);
 }
 
@@ -8413,11 +8441,18 @@ function VMI_set_bounds(args) {
       VM.variables[vi - 1][0],'] t = ', VM.t, ' LB = ', VM.sig4Dig(l),
       ', UB = ', VM.sig4Dig(u), fixed].join(''), l, u, inf_val);
   } else if(u < l) {
-    // Warn that "impossible" bounds would have been set...
-    const vk = vbl.displayName;
-    if(!VM.bound_issues[vk]) VM.bound_issues[vk] = [];
-    VM.bound_issues[vk].push(VM.t);
-    // ... and set LB to UB, so that lowest value is bounding.
+    // Check the difference, as this may be negligible.
+    if(u - l < VM.SIG_DIF_FROM_ZERO) {
+      u = Math.round(u * 1e5) / 1e5;
+      // NOTE: This may result in -0 (minus zero) => then set to 0.
+      if(u < 0 && u > -VM.NEAR_ZERO) u = 0;
+    } else {
+      // If substantial, warn that "impossible" bounds would have been set.
+      const vk = vbl.displayName;
+      if(!VM.bound_issues[vk]) VM.bound_issues[vk] = [];
+      VM.bound_issues[vk].push(VM.t);
+    }
+    // Set LB to UB, so that lowest value is bounding.
     l = u;
   }
   // NOTE: Since the VM vectors for lower bounds and upper bounds are