cacatoo

<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"> <style> .KEYW {color: #933;} .COMM {color: #bbb; font-style: italic;} .NUMB {color: #393;} .STRN {color: #393;} .REGX {color: #339;} .line {border-right: 1px dotted #666; color: #666; font-style: normal;} </style></head><body><pre> 1 'use strict'; 2 3 /** 4 * Gridpoint is what Gridmodels are made of. Contains everything that may happen in 1 locality. 5 */ 6 7 class Gridpoint { 8 /** 9 * The constructor function for a @Gridpoint object. Takes an optional template to copy primitives from. (NOTE!! Other types of objects are NOT deep copied by default) 10 * If you need synchronous updating with complex objects (for whatever reason), replate line 18 with line 19. This will slow things down quite a bit, so ony use this 11 * if you really need it. A better option is to use asynchronous updating so you won't have to worry about this at all :) 12 * @param {Gridpoint} template Optional template to make a new @Gridpoint from 13 */ 14 constructor(template) { 15 for (var prop in template) 16 this[prop] = template[prop]; // Shallow copy. It's fast, but be careful with syncronous updating! 17 // this[prop] = copy(template[prop]) // Deep copy. Takes much more time, and you'll likely end up copying much more than necessary. Use only if you're sure you need it! 18 } 19 } 20 21 /** 22 * Graph is a wrapper-class for a Dygraph element (see https://dygraphs.com/). It is attached to the DOM-windows, and stores all values to be plotted, colours, title, axis names, etc. 23 */ 24 25 class Graph { 26 /** 27 * The constructor function for a @Canvas object. 28 * @param {Array} labels array of strings containing the labels for datapoints (e.g. for the legend) 29 * @param {Array} values Array of floats to plot (here plotted over time) 30 * @param {Array} colours Array of colours to use for plotting 31 * @param {String} title Title of the plot 32 * @param {Object} opts dictionary-style list of opts to pass onto dygraphs 33 */ 34 constructor(labels, values, colours, title, opts) { 35 36 if (typeof window == undefined) throw "Using dygraphs with cashJS only works in browser-mode" 37 this.labels = labels; 38 this.data = [values]; 39 this.title = title; 40 this.num_dps = this.labels.length; // number of data points for this graphs 41 this.elem = document.createElement("div"); 42 this.elem.className = "graph-holder"; 43 this.colours = []; 44 45 for (let v in colours) { 46 if (v == "Time") continue 47 else if (colours[v] == undefined) this.colours.push("#000000"); 48 else if (colours[v][0] + colours[v][1] + colours[v][2] == 765) this.colours.push("#dddddd"); 49 else this.colours.push(rgbToHex(colours[v][0], colours[v][1], colours[v][2])); 50 } 51 52 document.body.appendChild(this.elem); 53 document.getElementById("graph_holder").appendChild(this.elem); 54 this.g = new Dygraph(this.elem, this.data, 55 { 56 title: this.title, 57 showRoller: false, 58 ylabel: this.labels.length == 2 ? this.labels[1] : "", 59 width: 500, 60 height: 200, 61 xlabel: this.labels[0], 62 drawPoints: opts && opts.drawPoints || false, 63 pointSize: opts ? (opts.pointSize ? opts.pointSize : 0) : 0, 64 strokePattern: opts ? (opts.strokePattern ? opts.strokePattern : null) : null, 65 dateWindow: [0, 100], 66 axisLabelFontSize: 10, 67 valueRange: [0.000,], 68 strokeWidth: opts ? opts.strokeWidth : 3, 69 colors: this.colours, 70 labels: this.labels 71 }); 72 } 73 74 75 /** Push data to your graph-element 76 * @param {array} array of floats to be added to the dygraph object (stored in 'data') 77 */ 78 push_data(data_array) { 79 this.data.push(data_array); 80 } 81 82 reset_plot() { 83 this.data = [this.data[0]]; 84 this.g.updateOptions( 85 { 86 'file': [this.data] 87 }); 88 } 89 90 /** 91 * Update the graph axes 92 */ 93 update() { 94 let max_x = 0; 95 let min_x = 999999999999; 96 for (let i of this.data) { 97 if (i[0] > max_x) max_x = i[0]; 98 if (i[0] < min_x) min_x = i[0]; 99 } 100 this.g.updateOptions( 101 { 102 'file': this.data, 103 dateWindow: [min_x, max_x] 104 }); 105 106 } 107 } 108 109 /* 110 Functions below are to make sure dygraphs understands the colours used by Cacatoo (converts to hex) 111 */ 112 function componentToHex(c) { 113 var hex = c.toString(16); 114 return hex.length == 1 ? "0" + hex : hex; 115 } 116 117 function rgbToHex(r, g, b) { 118 return "#" + componentToHex(r) + componentToHex(g) + componentToHex(b); 119 } 120 121 /** 122 * The ODE class is used to call the odex.js library and numerically solve ODEs 123 */ 124 125 class ODE { 126 /** 127 * The constructor function for a @ODE object. 128 * @param {function} eq Function that describes the ODE (see examples starting with ode) 129 * @param {Array} state_vector Initial state vector 130 * @param {Array} pars Array of parameters for the ODEs 131 * @param {Array} diff_rates Array of rates at which each state diffuses to neighbouring grid point (Has to be less than 0.25!) 132 * @param {String} ode_name Name of this ODE 133 */ 134 constructor(eq, state_vector, pars, diff_rates, ode_name, acceptable_error) { 135 this.name = ode_name; 136 this.eq = eq; 137 this.state = state_vector; 138 this.diff_rates = diff_rates; 139 this.pars = pars; 140 this.solver = new Solver(state_vector.length); 141 if (acceptable_error !== undefined) this.solver.absoluteTolerance = this.solver.relativeTolerance = acceptable_error; 142 } 143 144 /** 145 * Numerically solve the ODE 146 * @param {float} delta_t Step size 147 * @param {bool} opt_pos When enabled, negative values are set to 0 automatically 148 */ 149 solveTimestep(delta_t = 0.1, pos = false) { 150 let newstate = this.solver.solve( 151 this.eq(...this.pars), // function to solve and its pars (... unlists the array as a list of args) 152 0, // Initial x value 153 this.state, // Initial y value(s) 154 delta_t // Final x value 155 ).y; 156 if (pos) for (var i = 0; i < newstate.length; i++) if (newstate[i] < 0.000001) newstate[i] = 0.0; 157 this.state = newstate; 158 } 159 /** 160 * Prints the current state to the console 161 */ 162 print_state() { 163 console.log(this.state); 164 } 165 } 166 167 /* 168 I've wrapped Makoto Matsumoto and Takuji Nishimura's code in a namespace 169 so it's better encapsulated. Now you can have multiple random number generators 170 and they won't stomp all over eachother's state. 171 172 If you want to use this as a substitute for Math.random(), use the random() 173 method like so: 174 175 var m = new MersenneTwister(); 176 var randomNumber = m.random(); 177 178 You can also call the other genrand_{foo}() methods on the instance. 179 If you want to use a specific seed in order to get a repeatable random 180 sequence, pass an integer into the constructor: 181 var m = new MersenneTwister(123); 182 and that will always produce the same random sequence. 183 Sean McCullough (banksean@gmail.com) 184 */ 185 186 /* 187 A C-program for MT19937, with initialization improved 2002/1/26. 188 Coded by Takuji Nishimura and Makoto Matsumoto. 189 190 Before using, initialize the state by using init_genrand(seed) 191 or init_by_array(init_key, key_length). 192 193 Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura, 194 All rights reserved. 195 196 Redistribution and use in source and binary forms, with or without 197 modification, are permitted provided that the following conditions 198 are met: 199 200 1. Redistributions of source code must retain the above copyright 201 notice, this list of conditions and the following disclaimer. 202 203 2. Redistributions in binary form must reproduce the above copyright 204 notice, this list of conditions and the following disclaimer in the 205 documentation and/or other materials provided with the distribution. 206 207 3. The names of its contributors may not be used to endorse or promote 208 products derived from this software without specific prior written 209 permission. 210 211 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 212 "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT 213 LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR 214 A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR 215 CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, 216 EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 217 PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 218 PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF 219 LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING 220 NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS 221 SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 222 223 224 Any feedback is very welcome. 225 http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html 226 email: m-mat @ math.sci.hiroshima-u.ac.jp (remove space) 227 */ 228 229 function MersenneTwister(seed) { 230 if (seed == undefined) { 231 seed = new Date().getTime(); 232 } 233 /* Period parameters */ 234 this.N = 624; 235 this.M = 397; 236 this.MATRIX_A = 0x9908b0df; /* constant vector a */ 237 this.UPPER_MASK = 0x80000000; /* most significant w-r bits */ 238 this.LOWER_MASK = 0x7fffffff; /* least significant r bits */ 239 240 this.mt = new Array(this.N); /* the array for the state vector */ 241 this.mti=this.N+1; /* mti==N+1 means mt[N] is not initialized */ 242 243 this.init_genrand(seed); 244 } 245 246 /* initializes mt[N] with a seed */ 247 MersenneTwister.prototype.init_genrand = function(s) { 248 this.mt[</s