evalrational.js (4281B)
1 /** 2 * @license Apache-2.0 3 * 4 * Copyright (c) 2018 The Stdlib Authors. 5 * 6 * Licensed under the Apache License, Version 2.0 (the "License"); 7 * you may not use this file except in compliance with the License. 8 * You may obtain a copy of the License at 9 * 10 * http://www.apache.org/licenses/LICENSE-2.0 11 * 12 * Unless required by applicable law or agreed to in writing, software 13 * distributed under the License is distributed on an "AS IS" BASIS, 14 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 15 * See the License for the specific language governing permissions and 16 * limitations under the License. 17 */ 18 19 /* 20 * This script compiles modules for evaluating polynomial functions. If any polynomial coefficients change, this script should be rerun to update the compiled files. 21 */ 22 'use strict'; 23 24 // MODULES // 25 26 var resolve = require( 'path' ).resolve; 27 var writeFileSync = require( '@stdlib/fs/write-file' ).sync; 28 var licenseHeader = require( '@stdlib/_tools/licenses/header' ); 29 var compile = require( './../../../../base/tools/evalrational-compile' ); 30 31 32 // VARIABLES // 33 34 // Polynomial coefficients ordered in ascending degree... 35 36 // S(x) for small x: 37 var PS = [ 38 3.18016297876567817986E11, 39 -4.42979518059697779103E10, 40 2.54890880573376359104E9, 41 -6.29741486205862506537E7, 42 7.08840045257738576863E5, 43 -2.99181919401019853726E3, 44 0.0 45 ]; 46 47 var QS = [ 48 6.07366389490084639049E11, 49 2.24411795645340920940E10, 50 4.19320245898111231129E8, 51 5.17343888770096400730E6, 52 4.55847810806532581675E4, 53 2.81376268889994315696E2, 54 1.00000000000000000000E0 55 ]; 56 57 // C(x) for small x: 58 var PC = [ 59 9.99999999999999998822E-1, 60 -2.05525900955013891793E-1, 61 1.88843319396703850064E-2, 62 -6.45191435683965050962E-4, 63 9.50428062829859605134E-6, 64 -4.98843114573573548651E-8, 65 0.0 66 ]; 67 68 var QC = [ 69 1.00000000000000000118E0, 70 4.12142090722199792936E-2, 71 8.68029542941784300606E-4, 72 1.22262789024179030997E-5, 73 1.25001862479598821474E-7, 74 9.15439215774657478799E-10, 75 3.99982968972495980367E-12 76 ]; 77 78 // Auxiliary function f(x): 79 var PF = [ 80 3.76329711269987889006E-20, 81 1.34283276233062758925E-16, 82 1.72010743268161828879E-13, 83 1.02304514164907233465E-10, 84 3.05568983790257605827E-8, 85 4.63613749287867322088E-6, 86 3.45017939782574027900E-4, 87 1.15220955073585758835E-2, 88 1.43407919780758885261E-1, 89 4.21543555043677546506E-1, 90 0.0 91 ]; 92 93 var QF = [ 94 1.25443237090011264384E-20, 95 4.52001434074129701496E-17, 96 5.88754533621578410010E-14, 97 3.60140029589371370404E-11, 98 1.12699224763999035261E-8, 99 1.84627567348930545870E-6, 100 1.55934409164153020873E-4, 101 6.44051526508858611005E-3, 102 1.16888925859191382142E-1, 103 7.51586398353378947175E-1, 104 1.00000000000000000000E0 105 ]; 106 107 // Auxiliary function g(x): 108 var PG = [ 109 1.86958710162783235106E-22, 110 8.36354435630677421531E-19, 111 1.37555460633261799868E-15, 112 1.08268041139020870318E-12, 113 4.45344415861750144738E-10, 114 9.82852443688422223854E-8, 115 1.15138826111884280931E-5, 116 6.84079380915393090172E-4, 117 1.87648584092575249293E-2, 118 1.97102833525523411709E-1, 119 5.04442073643383265887E-1, 120 0.0 121 ]; 122 123 var QG = [ 124 1.86958710162783236342E-22, 125 8.39158816283118707363E-19, 126 1.38796531259578871258E-15, 127 1.10273215066240270757E-12, 128 4.60680728146520428211E-10, 129 1.04314589657571990585E-7, 130 1.27545075667729118702E-5, 131 8.14679107184306179049E-4, 132 2.53603741420338795122E-2, 133 3.37748989120019970451E-1, 134 1.47495759925128324529E0, 135 1.00000000000000000000E0 136 ]; 137 138 // Header to add to output files: 139 var header = licenseHeader( 'Apache-2.0', 'js', { 140 'year': ( new Date() ).getFullYear(), 141 'copyright': 'The Stdlib Authors' 142 }); 143 header += '\n/* This is a generated file. Do not edit directly. */\n'; 144 145 146 // MAIN // 147 148 /** 149 * Main execution sequence. 150 * 151 * @private 152 */ 153 function main() { 154 var fpath; 155 var opts; 156 var str; 157 158 opts = { 159 'encoding': 'utf8' 160 }; 161 162 fpath = resolve( __dirname, '..', 'lib', 'rational_psqs.js' ); 163 str = header + compile( PS, QS ); 164 writeFileSync( fpath, str, opts ); 165 166 fpath = resolve( __dirname, '..', 'lib', 'rational_pcqc.js' ); 167 str = header + compile( PC, QC ); 168 writeFileSync( fpath, str, opts ); 169 170 fpath = resolve( __dirname, '..', 'lib', 'rational_pfqf.js' ); 171 str = header + compile( PF, QF ); 172 writeFileSync( fpath, str, opts ); 173 174 fpath = resolve( __dirname, '..', 'lib', 'rational_pgqg.js' ); 175 str = header + compile( PG, QG ); 176 writeFileSync( fpath, str, opts ); 177 } 178 179 main();