star-cmap

Mapping values onto color ramps
git clone git://git.meso-star.fr/star-cmap.git
Log | Files | Refs | README | LICENSE

scmap_palettes.c (61505B)


      1 /* Copyright (C) 2020, 2021, 2023 |Méso|Star> (contact@meso-star.com)
      2  *
      3  * This program is free software: you can redistribute it and/or modify
      4  * it under the terms of the GNU General Public License as published by
      5  * the Free Software Foundation, either version 3 of the License, or
      6  * (at your option) any later version.
      7  *
      8  * This program is distributed in the hope that it will be useful,
      9  * but WITHOUT ANY WARRANTY; without even the implied warranty of
     10  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
     11  * GNU General Public License for more details.
     12  *
     13  * You should have received a copy of the GNU General Public License
     14  * along with this program. If not, see <http://www.gnu.org/licenses/>. */
     15 
     16 #include "scmap.h"
     17 #include <string.h>
     18 
     19 static const double accent[] = {
     20   0.498039, 0.788235, 0.498039,
     21   0.745098, 0.682353, 0.831373,
     22   0.992157, 0.752941, 0.525490,
     23   1.000000, 1.000000, 0.600000,
     24   0.219608, 0.423529, 0.690196,
     25   0.941176, 0.007843, 0.498039,
     26   0.749020, 0.356863, 0.090196,
     27   0.400000, 0.400000, 0.400000
     28 };
     29 
     30 static const double blues[] = {
     31   0.968627, 0.984314, 1.000000,
     32   0.870588, 0.921569, 0.968627,
     33   0.776471, 0.858824, 0.937255,
     34   0.619608, 0.792157, 0.882353,
     35   0.419608, 0.682353, 0.839216,
     36   0.258824, 0.572549, 0.776471,
     37   0.129412, 0.443137, 0.709804,
     38   0.031373, 0.270588, 0.580392
     39 };
     40 
     41 static const double brbg[] = {
     42   0.549020, 0.317647, 0.039216,
     43   0.749020, 0.505882, 0.176471,
     44   0.874510, 0.760784, 0.490196,
     45   0.964706, 0.909804, 0.764706,
     46   0.780392, 0.917647, 0.898039,
     47   0.501961, 0.803922, 0.756863,
     48   0.207843, 0.592157, 0.560784,
     49   0.003922, 0.400000, 0.368627
     50 };
     51 
     52 static const double bugn[] = {
     53   0.968627, 0.988235, 0.992157,
     54   0.898039, 0.960784, 0.976471,
     55   0.800000, 0.925490, 0.901961,
     56   0.600000, 0.847059, 0.788235,
     57   0.400000, 0.760784, 0.643137,
     58   0.254902, 0.682353, 0.462745,
     59   0.137255, 0.545098, 0.270588,
     60   0.000000, 0.345098, 0.141176
     61 };
     62 
     63 static const double bupu[] = {
     64   0.968627, 0.988235, 0.992157,
     65   0.878431, 0.925490, 0.956863,
     66   0.749020, 0.827451, 0.901961,
     67   0.619608, 0.737255, 0.854902,
     68   0.549020, 0.588235, 0.776471,
     69   0.549020, 0.419608, 0.694118,
     70   0.533333, 0.254902, 0.615686,
     71   0.431373, 0.003922, 0.419608
     72 };
     73 
     74 static const double chromajs[] = {
     75   1.000000, 1.000000, 0.878431,
     76   1.000000, 0.874510, 0.721569,
     77   1.000000, 0.737255, 0.580392,
     78   1.000000, 0.592157, 0.466667,
     79   1.000000, 0.411765, 0.384314,
     80   0.933333, 0.258824, 0.337255,
     81   0.823529, 0.121569, 0.278431,
     82   0.690196, 0.023529, 0.172549,
     83   0.545098, 0.000000, 0.000000
     84 };
     85 
     86 static const double dark2[] = {
     87   0.105882, 0.619608, 0.466667,
     88   0.850980, 0.372549, 0.007843,
     89   0.458824, 0.439216, 0.701961,
     90   0.905882, 0.160784, 0.541176,
     91   0.400000, 0.650980, 0.117647,
     92   0.901961, 0.670588, 0.007843,
     93   0.650980, 0.462745, 0.113725,
     94   0.400000, 0.400000, 0.400000
     95 };
     96 
     97 static const double gnbu[] = {
     98   0.968627, 0.988235, 0.941176,
     99   0.878431, 0.952941, 0.858824,
    100   0.800000, 0.921569, 0.772549,
    101   0.658824, 0.866667, 0.709804,
    102   0.482353, 0.800000, 0.768627,
    103   0.305882, 0.701961, 0.827451,
    104   0.168627, 0.549020, 0.745098,
    105   0.031373, 0.345098, 0.619608
    106 };
    107 
    108 static const double gnpu[] = {
    109   0.223529, 0.388235, 0.325490,
    110   0.050980, 0.694118, 0.294118,
    111   0.427451, 0.752941, 0.403922,
    112   0.670588, 0.839216, 0.607843,
    113   0.854902, 0.917647, 0.756863,
    114   0.874510, 0.800000, 0.894118,
    115   0.780392, 0.698039, 0.839216,
    116   0.580392, 0.454902, 0.705882,
    117   0.458824, 0.250980, 0.596078,
    118   0.313725, 0.286275, 0.443137
    119 };
    120 
    121 static const double greens[] = {
    122   0.968627, 0.988235, 0.960784,
    123   0.898039, 0.960784, 0.878431,
    124   0.780392, 0.913725, 0.752941,
    125   0.631373, 0.850980, 0.607843,
    126   0.454902, 0.768627, 0.462745,
    127   0.254902, 0.670588, 0.364706,
    128   0.137255, 0.545098, 0.270588,
    129   0.000000, 0.352941, 0.196078
    130 };
    131 
    132 static const double greys[] = {
    133   1.000000, 1.000000, 1.000000,
    134   0.941176, 0.941176, 0.941176,
    135   0.850980, 0.850980, 0.850980,
    136   0.741176, 0.741176, 0.741176,
    137   0.588235, 0.588235, 0.588235,
    138   0.450980, 0.450980, 0.450980,
    139   0.321569, 0.321569, 0.321569,
    140   0.145098, 0.145098, 0.145098
    141 };
    142 
    143 static const double oranges[] = {
    144   1.000000, 0.960784, 0.921569,
    145   0.996078, 0.901961, 0.807843,
    146   0.992157, 0.815686, 0.635294,
    147   0.992157, 0.682353, 0.419608,
    148   0.992157, 0.552941, 0.235294,
    149   0.945098, 0.411765, 0.074510,
    150   0.850980, 0.282353, 0.003922,
    151   0.549020, 0.176471, 0.015686
    152 };
    153 
    154 static const double orrd[] = {
    155   1.000000, 0.968627, 0.925490,
    156   0.996078, 0.909804, 0.784314,
    157   0.992157, 0.831373, 0.619608,
    158   0.992157, 0.733333, 0.517647,
    159   0.988235, 0.552941, 0.349020,
    160   0.937255, 0.396078, 0.282353,
    161   0.843137, 0.188235, 0.121569,
    162   0.600000, 0.000000, 0.000000
    163 };
    164 
    165 static const double paired[] = {
    166   0.650980, 0.807843, 0.890196,
    167   0.121569, 0.470588, 0.705882,
    168   0.698039, 0.874510, 0.541176,
    169   0.200000, 0.627451, 0.172549,
    170   0.984314, 0.603922, 0.600000,
    171   0.890196, 0.101961, 0.109804,
    172   0.992157, 0.749020, 0.435294,
    173   1.000000, 0.498039, 0.000000
    174 };
    175 
    176 static const double parula[] = {
    177   0.207843, 0.164706, 0.529412,
    178   0.011765, 0.388235, 0.882353,
    179   0.078431, 0.521569, 0.831373,
    180   0.023529, 0.654902, 0.776471,
    181   0.219608, 0.725490, 0.619608,
    182   0.572549, 0.749020, 0.450980,
    183   0.850980, 0.729412, 0.337255,
    184   0.988235, 0.807843, 0.180392,
    185   0.976471, 0.984314, 0.054902
    186 };
    187 
    188 static const double pastel1[] = {
    189   0.984314, 0.705882, 0.682353,
    190   0.701961, 0.803922, 0.890196,
    191   0.800000, 0.921569, 0.772549,
    192   0.870588, 0.796078, 0.894118,
    193   0.996078, 0.850980, 0.650980,
    194   1.000000, 1.000000, 0.800000,
    195   0.898039, 0.847059, 0.741176,
    196   0.992157, 0.854902, 0.925490
    197 };
    198 
    199 static const double pastel2[] = {
    200   0.701961, 0.886275, 0.803922,
    201   0.992157, 0.803922, 0.674510,
    202   0.803922, 0.709804, 0.909804,
    203   0.956863, 0.792157, 0.894118,
    204   0.839216, 0.960784, 0.788235,
    205   1.000000, 0.949020, 0.682353,
    206   0.945098, 0.886275, 0.800000,
    207   0.800000, 0.800000, 0.800000
    208 };
    209 
    210 static const double piyg[] = {
    211   0.772549, 0.105882, 0.490196,
    212   0.870588, 0.466667, 0.682353,
    213   0.945098, 0.713725, 0.854902,
    214   0.992157, 0.878431, 0.937255,
    215   0.901961, 0.960784, 0.815686,
    216   0.721569, 0.882353, 0.525490,
    217   0.498039, 0.737255, 0.254902,
    218   0.301961, 0.572549, 0.129412
    219 };
    220 
    221 static const double prgn[] = {
    222   0.462745, 0.164706, 0.513725,
    223   0.600000, 0.439216, 0.670588,
    224   0.760784, 0.647059, 0.811765,
    225   0.905882, 0.831373, 0.909804,
    226   0.850980, 0.941176, 0.827451,
    227   0.650980, 0.858824, 0.627451,
    228   0.352941, 0.682353, 0.380392,
    229   0.105882, 0.470588, 0.215686
    230 };
    231 
    232 static const double pubu[] = {
    233   1.000000, 0.968627, 0.984314,
    234   0.925490, 0.905882, 0.949020,
    235   0.815686, 0.819608, 0.901961,
    236   0.650980, 0.741176, 0.858824,
    237   0.454902, 0.662745, 0.811765,
    238   0.211765, 0.564706, 0.752941,
    239   0.019608, 0.439216, 0.690196,
    240   0.011765, 0.305882, 0.482353
    241 };
    242 
    243 static const double pubugn[] = {
    244   1.000000, 0.968627, 0.984314,
    245   0.925490, 0.905882, 0.941176,
    246   0.815686, 0.819608, 0.901961,
    247   0.650980, 0.741176, 0.858824,
    248   0.403922, 0.662745, 0.811765,
    249   0.211765, 0.564706, 0.752941,
    250   0.007843, 0.505882, 0.541176,
    251   0.003922, 0.396078, 0.250980
    252 };
    253 
    254 static const double puor[] = {
    255   0.701961, 0.345098, 0.023529,
    256   0.878431, 0.509804, 0.078431,
    257   0.992157, 0.721569, 0.388235,
    258   0.996078, 0.878431, 0.713725,
    259   0.847059, 0.854902, 0.921569,
    260   0.698039, 0.670588, 0.823529,
    261   0.501961, 0.450980, 0.674510,
    262   0.329412, 0.152941, 0.533333
    263 };
    264 
    265 static const double purd[] = {
    266   0.968627, 0.956863, 0.976471,
    267   0.905882, 0.882353, 0.937255,
    268   0.831373, 0.725490, 0.854902,
    269   0.788235, 0.580392, 0.780392,
    270   0.874510, 0.396078, 0.690196,
    271   0.905882, 0.160784, 0.541176,
    272   0.807843, 0.070588, 0.337255,
    273   0.568627, 0.000000, 0.247059
    274 };
    275 
    276 static const double purples[] = {
    277   0.988235, 0.984314, 0.992157,
    278   0.937255, 0.929412, 0.960784,
    279   0.854902, 0.854902, 0.921569,
    280   0.737255, 0.741176, 0.862745,
    281   0.619608, 0.603922, 0.784314,
    282   0.501961, 0.490196, 0.729412,
    283   0.415686, 0.317647, 0.639216,
    284   0.290196, 0.078431, 0.525490
    285 };
    286 
    287 static const double rdbu[] = {
    288   0.698039, 0.094118, 0.168627,
    289   0.839216, 0.376471, 0.301961,
    290   0.956863, 0.647059, 0.509804,
    291   0.992157, 0.858824, 0.780392,
    292   0.819608, 0.898039, 0.941176,
    293   0.572549, 0.772549, 0.870588,
    294   0.262745, 0.576471, 0.764706,
    295   0.129412, 0.400000, 0.674510
    296 };
    297 
    298 static const double rdgy[] = {
    299   0.698039, 0.094118, 0.168627,
    300   0.839216, 0.376471, 0.301961,
    301   0.956863, 0.647059, 0.509804,
    302   0.992157, 0.858824, 0.780392,
    303   0.878431, 0.878431, 0.878431,
    304   0.729412, 0.729412, 0.729412,
    305   0.529412, 0.529412, 0.529412,
    306   0.301961, 0.301961, 0.301961
    307 };
    308 
    309 static const double rdpu[] = {
    310   1.000000, 0.968627, 0.952941,
    311   0.992157, 0.878431, 0.866667,
    312   0.988235, 0.772549, 0.752941,
    313   0.980392, 0.623529, 0.709804,
    314   0.968627, 0.407843, 0.631373,
    315   0.866667, 0.203922, 0.592157,
    316   0.682353, 0.003922, 0.494118,
    317   0.478431, 0.003922, 0.466667
    318 };
    319 
    320 static const double rdylbu[] = {
    321   0.843137, 0.188235, 0.152941,
    322   0.956863, 0.427451, 0.262745,
    323   0.992157, 0.682353, 0.380392,
    324   0.996078, 0.878431, 0.564706,
    325   0.878431, 0.952941, 0.972549,
    326   0.670588, 0.850980, 0.913725,
    327   0.454902, 0.678431, 0.819608,
    328   0.270588, 0.458824, 0.705882
    329 };
    330 
    331 static const double rdylgn[] = {
    332   0.843137, 0.188235, 0.152941,
    333   0.956863, 0.427451, 0.262745,
    334   0.992157, 0.682353, 0.380392,
    335   0.996078, 0.878431, 0.545098,
    336   0.850980, 0.937255, 0.545098,
    337   0.650980, 0.850980, 0.415686,
    338   0.400000, 0.741176, 0.388235,
    339   0.101961, 0.596078, 0.313725
    340 };
    341 
    342 static const double reds[] = {
    343   1.000000, 0.960784, 0.941176,
    344   0.996078, 0.878431, 0.823529,
    345   0.988235, 0.733333, 0.631373,
    346   0.988235, 0.572549, 0.447059,
    347   0.984314, 0.415686, 0.290196,
    348   0.937255, 0.231373, 0.172549,
    349   0.796078, 0.094118, 0.113725,
    350   0.600000, 0.000000, 0.050980
    351 };
    352 
    353 static const double sand[] = {
    354   0.376471, 0.282353, 0.376471,
    355   0.470588, 0.282353, 0.376471,
    356   0.658824, 0.376471, 0.376471,
    357   0.752941, 0.470588, 0.376471,
    358   0.941176, 0.658824, 0.282353,
    359   0.972549, 0.792157, 0.549020,
    360   0.996078, 0.925490, 0.682353,
    361   1.000000, 0.956863, 0.760784,
    362   1.000000, 0.968627, 0.858824,
    363   1.000000, 0.988235, 0.964706
    364 };
    365 
    366 static const double set1[] = {
    367   0.894118, 0.101961, 0.109804,
    368   0.215686, 0.494118, 0.721569,
    369   0.301961, 0.686275, 0.290196,
    370   0.596078, 0.305882, 0.639216,
    371   1.000000, 0.498039, 0.000000,
    372   1.000000, 1.000000, 0.200000,
    373   0.650980, 0.337255, 0.156863,
    374   0.968627, 0.505882, 0.749020
    375 };
    376 
    377 static const double set2[] = {
    378   0.400000, 0.760784, 0.647059,
    379   0.988235, 0.552941, 0.384314,
    380   0.552941, 0.627451, 0.796078,
    381   0.905882, 0.541176, 0.764706,
    382   0.650980, 0.847059, 0.329412,
    383   1.000000, 0.850980, 0.184314,
    384   0.898039, 0.768627, 0.580392,
    385   0.701961, 0.701961, 0.701961
    386 };
    387 
    388 static const double set3[] = {
    389   0.552941, 0.827451, 0.780392,
    390   1.000000, 1.000000, 0.701961,
    391   0.745098, 0.729412, 0.854902,
    392   0.984314, 0.501961, 0.447059,
    393   0.501961, 0.694118, 0.827451,
    394   0.992157, 0.705882, 0.384314,
    395   0.701961, 0.870588, 0.411765,
    396   0.988235, 0.803922, 0.898039
    397 };
    398 
    399 static const double spectral[] = {
    400   0.835294, 0.243137, 0.309804,
    401   0.956863, 0.427451, 0.262745,
    402   0.992157, 0.682353, 0.380392,
    403   0.996078, 0.878431, 0.545098,
    404   0.901961, 0.960784, 0.596078,
    405   0.670588, 0.866667, 0.643137,
    406   0.400000, 0.760784, 0.647059,
    407   0.196078, 0.533333, 0.741176
    408 };
    409 
    410 static const double whgnbu[] = {
    411   1.000000, 1.000000, 1.000000,
    412   0.878431, 0.952941, 0.858824,
    413   0.800000, 0.921569, 0.772549,
    414   0.658824, 0.866667, 0.709804,
    415   0.482353, 0.800000, 0.768627,
    416   0.305882, 0.701961, 0.827451,
    417   0.168627, 0.549020, 0.745098,
    418   0.031373, 0.345098, 0.619608
    419 };
    420 
    421 static const double whylrd[] = {
    422   1.000000, 1.000000, 1.000000,
    423   1.000000, 0.933333, 0.000000,
    424   1.000000, 0.439216, 0.000000,
    425   0.933333, 0.000000, 0.000000,
    426   0.498039, 0.000000, 0.000000
    427 };
    428 
    429 static const double ylgn[] = {
    430   1.000000, 1.000000, 0.898039,
    431   0.968627, 0.988235, 0.725490,
    432   0.850980, 0.941176, 0.639216,
    433   0.678431, 0.866667, 0.556863,
    434   0.470588, 0.776471, 0.474510,
    435   0.254902, 0.670588, 0.364706,
    436   0.137255, 0.517647, 0.262745,
    437   0.000000, 0.352941, 0.196078
    438 };
    439 
    440 static const double ylgnbu[] = {
    441   1.000000, 1.000000, 0.850980,
    442   0.929412, 0.972549, 0.694118,
    443   0.780392, 0.913725, 0.705882,
    444   0.498039, 0.803922, 0.733333,
    445   0.254902, 0.713725, 0.768627,
    446   0.113725, 0.568627, 0.752941,
    447   0.133333, 0.368627, 0.658824,
    448   0.047059, 0.172549, 0.517647
    449 };
    450 
    451 static const double ylorbr[] = {
    452   1.000000, 1.000000, 0.898039,
    453   1.000000, 0.968627, 0.737255,
    454   0.996078, 0.890196, 0.568627,
    455   0.996078, 0.768627, 0.309804,
    456   0.996078, 0.600000, 0.160784,
    457   0.925490, 0.439216, 0.078431,
    458   0.800000, 0.298039, 0.007843,
    459   0.549020, 0.176471, 0.015686
    460 };
    461 
    462 static const double ylorrd[] = {
    463   1.000000, 1.000000, 0.800000,
    464   1.000000, 0.929412, 0.627451,
    465   0.996078, 0.850980, 0.462745,
    466   0.996078, 0.698039, 0.298039,
    467   0.992157, 0.552941, 0.235294,
    468   0.988235, 0.305882, 0.164706,
    469   0.890196, 0.101961, 0.109804,
    470   0.694118, 0.000000, 0.149020
    471 };
    472 
    473 static const double ylrd[] = {
    474   1.000000, 0.933333, 0.000000,
    475   1.000000, 0.439216, 0.000000,
    476   0.933333, 0.000000, 0.000000,
    477   0.498039, 0.000000, 0.000000
    478 };
    479 
    480 static const double inferno[] = {
    481   0.001462, 0.000466, 0.013866,
    482   0.002267, 0.001270, 0.018570,
    483   0.003299, 0.002249, 0.024239,
    484   0.004547, 0.003392, 0.030909,
    485   0.006006, 0.004692, 0.038558,
    486   0.007676, 0.006136, 0.046836,
    487   0.009561, 0.007713, 0.055143,
    488   0.011663, 0.009417, 0.063460,
    489   0.013995, 0.011225, 0.071862,
    490   0.016561, 0.013136, 0.080282,
    491   0.019373, 0.015133, 0.088767,
    492   0.022447, 0.017199, 0.097327,
    493   0.025793, 0.019331, 0.105930,
    494   0.029432, 0.021503, 0.114621,
    495   0.033385, 0.023702, 0.123397,
    496   0.037668, 0.025921, 0.132232,
    497   0.042253, 0.028139, 0.141141,
    498   0.046915, 0.030324, 0.150164,
    499   0.051644, 0.032474, 0.159254,
    500   0.056449, 0.034569, 0.168414,
    501   0.061340, 0.036590, 0.177642,
    502   0.066331, 0.038504, 0.186962,
    503   0.071429, 0.040294, 0.196354,
    504   0.076637, 0.041905, 0.205799,
    505   0.081962, 0.043328, 0.215289,
    506   0.087411, 0.044556, 0.224813,
    507   0.092990, 0.045583, 0.234358,
    508   0.098702, 0.046402, 0.243904,
    509   0.104551, 0.047008, 0.253430,
    510   0.110536, 0.047399, 0.262912,
    511   0.116656, 0.047574, 0.272321,
    512   0.122908, 0.047536, 0.281624,
    513   0.129285, 0.047293, 0.290788,
    514   0.135778, 0.046856, 0.299776,
    515   0.142378, 0.046242, 0.308553,
    516   0.149073, 0.045468, 0.317085,
    517   0.155850, 0.044559, 0.325338,
    518   0.162689, 0.043554, 0.333277,
    519   0.169575, 0.042489, 0.340874,
    520   0.176493, 0.041402, 0.348111,
    521   0.183429, 0.040329, 0.354971,
    522   0.190367, 0.039309, 0.361447,
    523   0.197297, 0.038400, 0.367535,
    524   0.204209, 0.037632, 0.373238,
    525   0.211095, 0.037030, 0.378563,
    526   0.217949, 0.036615, 0.383522,
    527   0.224763, 0.036405, 0.388129,
    528   0.231538, 0.036405, 0.392400,
    529   0.238273, 0.036621, 0.396353,
    530   0.244967, 0.037055, 0.400007,
    531   0.251620, 0.037705, 0.403378,
    532   0.258234, 0.038571, 0.406485,
    533   0.264810, 0.039647, 0.409345,
    534   0.271347, 0.040922, 0.411976,
    535   0.277850, 0.042353, 0.414392,
    536   0.284321, 0.043933, 0.416608,
    537   0.290763, 0.045644, 0.418637,
    538   0.297178, 0.047470, 0.420491,
    539   0.303568, 0.049396, 0.422182,
    540   0.309935, 0.051407, 0.423721,
    541   0.316282, 0.053490, 0.425116,
    542   0.322610, 0.055634, 0.426377,
    543   0.328921, 0.057827, 0.427511,
    544   0.335217, 0.060060, 0.428524,
    545   0.341500, 0.062325, 0.429425,
    546   0.347771, 0.064616, 0.430217,
    547   0.354032, 0.066925, 0.430906,
    548   0.360284, 0.069247, 0.431497,
    549   0.366529, 0.071579, 0.431994,
    550   0.372768, 0.073915, 0.432400,
    551   0.379001, 0.076253, 0.432719,
    552   0.385228, 0.078591, 0.432955,
    553   0.391453, 0.080927, 0.433109,
    554   0.397674, 0.083257, 0.433183,
    555   0.403894, 0.085580, 0.433179,
    556   0.410113, 0.087896, 0.433098,
    557   0.416331, 0.090203, 0.432943,
    558   0.422549, 0.092501, 0.432714,
    559   0.428768, 0.094790, 0.432412,
    560   0.434987, 0.097069, 0.432039,
    561   0.441207, 0.099338, 0.431594,
    562   0.447428, 0.101597, 0.431080,
    563   0.453651, 0.103848, 0.430498,
    564   0.459875, 0.106089, 0.429846,
    565   0.466100, 0.108322, 0.429125,
    566   0.472328, 0.110547, 0.428334,
    567   0.478558, 0.112764, 0.427475,
    568   0.484789, 0.114974, 0.426548,
    569   0.491022, 0.117179, 0.425552,
    570   0.497257, 0.119379, 0.424488,
    571   0.503493, 0.121575, 0.423356,
    572   0.509730, 0.123769, 0.422156,
    573   0.515967, 0.125960, 0.420887,
    574   0.522206, 0.128150, 0.419549,
    575   0.528444, 0.130341, 0.418142,
    576   0.534683, 0.132534, 0.416667,
    577   0.540920, 0.134729, 0.415123,
    578   0.547157, 0.136929, 0.413511,
    579   0.553392, 0.139134, 0.411829,
    580   0.559624, 0.141346, 0.410078,
    581   0.565854, 0.143567, 0.408258,
    582   0.572081, 0.145797, 0.406369,
    583   0.578304, 0.148039, 0.404411,
    584   0.584521, 0.150294, 0.402385,
    585   0.590734, 0.152563, 0.400290,
    586   0.596940, 0.154848, 0.398125,
    587   0.603139, 0.157151, 0.395891,
    588   0.609330, 0.159474, 0.393589,
    589   0.615513, 0.161817, 0.391219,
    590   0.621685, 0.164184, 0.388781,
    591   0.627847, 0.166575, 0.386276,
    592   0.633998, 0.168992, 0.383704,
    593   0.640135, 0.171438, 0.381065,
    594   0.646260, 0.173914, 0.378359,
    595   0.652369, 0.176421, 0.375586,
    596   0.658463, 0.178962, 0.372748,
    597   0.664540, 0.181539, 0.369846,
    598   0.670599, 0.184153, 0.366879,
    599   0.676638, 0.186807, 0.363849,
    600   0.682656, 0.189501, 0.360757,
    601   0.688653, 0.192239, 0.357603,
    602   0.694627, 0.195021, 0.354388,
    603   0.700576, 0.197851, 0.351113,
    604   0.706500, 0.200728, 0.347777,
    605   0.712396, 0.203656, 0.344383,
    606   0.718264, 0.206636, 0.340931,
    607   0.724103, 0.209670, 0.337424,
    608   0.729909, 0.212759, 0.333861,
    609   0.735683, 0.215906, 0.330245,
    610   0.741423, 0.219112, 0.326576,
    611   0.747127, 0.222378, 0.322856,
    612   0.752794, 0.225706, 0.319085,
    613   0.758422, 0.229097, 0.315266,
    614   0.764010, 0.232554, 0.311399,
    615   0.769556, 0.236077, 0.307485,
    616   0.775059, 0.239667, 0.303526,
    617   0.780517, 0.243327, 0.299523,
    618   0.785929, 0.247056, 0.295477,
    619   0.791293, 0.250856, 0.291390,
    620   0.796607, 0.254728, 0.287264,
    621   0.801871, 0.258674, 0.283099,
    622   0.807082, 0.262692, 0.278898,
    623   0.812239, 0.266786, 0.274661,
    624   0.817341, 0.270954, 0.270390,
    625   0.822386, 0.275197, 0.266085,
    626   0.827372, 0.279517, 0.261750,
    627   0.832299, 0.283913, 0.257383,
    628   0.837165, 0.288385, 0.252988,
    629   0.841969, 0.292933, 0.248564,
    630   0.846709, 0.297559, 0.244113,
    631   0.851384, 0.302260, 0.239636,
    632   0.855992, 0.307038, 0.235133,
    633   0.860533, 0.311892, 0.230606,
    634   0.865006, 0.316822, 0.226055,
    635   0.869409, 0.321827, 0.221482,
    636   0.873741, 0.326906, 0.216886,
    637   0.878001, 0.332060, 0.212268,
    638   0.882188, 0.337287, 0.207628,
    639   0.886302, 0.342586, 0.202968,
    640   0.890341, 0.347957, 0.198286,
    641   0.894305, 0.353399, 0.193584,
    642   0.898192, 0.358911, 0.188860,
    643   0.902003, 0.364492, 0.184116,
    644   0.905735, 0.370140, 0.179350,
    645   0.909390, 0.375856, 0.174563,
    646   0.912966, 0.381636, 0.169755,
    647   0.916462, 0.387481, 0.164924,
    648   0.919879, 0.393389, 0.160070,
    649   0.923215, 0.399359, 0.155193,
    650   0.926470, 0.405389, 0.150292,
    651   0.929644, 0.411479, 0.145367,
    652   0.932737, 0.417627, 0.140417,
    653   0.935747, 0.423831, 0.135440,
    654   0.938675, 0.430091, 0.130438,
    655   0.941521, 0.436405, 0.125409,
    656   0.944285, 0.442772, 0.120354,
    657   0.946965, 0.449191, 0.115272,
    658   0.949562, 0.455660, 0.110164,
    659   0.952075, 0.462178, 0.105031,
    660   0.954506, 0.468744, 0.099874,
    661   0.956852, 0.475356, 0.094695,
    662   0.959114, 0.482014, 0.089499,
    663   0.961293, 0.488716, 0.084289,
    664   0.963387, 0.495462, 0.079073,
    665   0.965397, 0.502249, 0.073859,
    666   0.967322, 0.509078, 0.068659,
    667   0.969163, 0.515946, 0.063488,
    668   0.970919, 0.522853, 0.058367,
    669   0.972590, 0.529798, 0.053324,
    670   0.974176, 0.536780, 0.048392,
    671   0.975677, 0.543798, 0.043618,
    672   0.977092, 0.550850, 0.039050,
    673   0.978422, 0.557937, 0.034931,
    674   0.979666, 0.565057, 0.031409,
    675   0.980824, 0.572209, 0.028508,
    676   0.981895, 0.579392, 0.026250,
    677   0.982881, 0.586606, 0.024661,
    678   0.983779, 0.593849, 0.023770,
    679   0.984591, 0.601122, 0.023606,
    680   0.985315, 0.608422, 0.024202,
    681   0.985952, 0.615750, 0.025592,
    682   0.986502, 0.623105, 0.027814,
    683   0.986964, 0.630485, 0.030908,
    684   0.987337, 0.637890, 0.034916,
    685   0.987622, 0.645320, 0.039886,
    686   0.987819, 0.652773, 0.045581,
    687   0.987926, 0.660250, 0.051750,
    688   0.987945, 0.667748, 0.058329,
    689   0.987874, 0.675267, 0.065257,
    690   0.987714, 0.682807, 0.072489,
    691   0.987464, 0.690366, 0.079990,
    692   0.987124, 0.697944, 0.087731,
    693   0.986694, 0.705540, 0.095694,
    694   0.986175, 0.713153, 0.103863,
    695   0.985566, 0.720782, 0.112229,
    696   0.984865, 0.728427, 0.120785,
    697   0.984075, 0.736087, 0.129527,
    698   0.983196, 0.743758, 0.138453,
    699   0.982228, 0.751442, 0.147565,
    700   0.981173, 0.759135, 0.156863,
    701   0.980032, 0.766837, 0.166353,
    702   0.978806, 0.774545, 0.176037,
    703   0.977497, 0.782258, 0.185923,
    704   0.976108, 0.789974, 0.196018,
    705   0.974638, 0.797692, 0.206332,
    706   0.973088, 0.805409, 0.216877,
    707   0.971468, 0.813122, 0.227658,
    708   0.969783, 0.820825, 0.238686,
    709   0.968041, 0.828515, 0.249972,
    710   0.966243, 0.836191, 0.261534,
    711   0.964394, 0.843848, 0.273391,
    712   0.962517, 0.851476, 0.285546,
    713   0.960626, 0.859069, 0.298010,
    714   0.958720, 0.866624, 0.310820,
    715   0.956834, 0.874129, 0.323974,
    716   0.954997, 0.881569, 0.337475,
    717   0.953215, 0.888942, 0.351369,
    718   0.951546, 0.896226, 0.365627,
    719   0.950018, 0.903409, 0.380271,
    720   0.948683, 0.910473, 0.395289,
    721   0.947594, 0.917399, 0.410665,
    722   0.946809, 0.924168, 0.426373,
    723   0.946392, 0.930761, 0.442367,
    724   0.946403, 0.937159, 0.458592,
    725   0.946903, 0.943348, 0.474970,
    726   0.947937, 0.949318, 0.491426,
    727   0.949545, 0.955063, 0.507860,
    728   0.951740, 0.960587, 0.524203,
    729   0.954529, 0.965896, 0.540361,
    730   0.957896, 0.971003, 0.556275,
    731   0.961812, 0.975924, 0.571925,
    732   0.966249, 0.980678, 0.587206,
    733   0.971162, 0.985282, 0.602154,
    734   0.976511, 0.989753, 0.616760,
    735   0.982257, 0.994109, 0.631017,
    736   0.988362, 0.998364, 0.644924
    737 };
    738 
    739 static const double jet[] = {
    740   0.0, 0.0, 0.5,
    741   0.0, 0.0, 1.0,
    742   0.0, 0.5, 1.0,
    743   0.0, 1.0, 1.0,
    744   0.5, 1.0, 0.5,
    745   1.0, 1.0, 0.0,
    746   1.0, 0.5, 0.0,
    747   1.0, 0.0, 0.0,
    748   0.5, 0.0, 0.0,
    749 };
    750 
    751 static const double magma[] = {
    752   0.001462, 0.000466, 0.013866,
    753   0.002258, 0.001295, 0.018331,
    754   0.003279, 0.002305, 0.023708,
    755   0.004512, 0.003490, 0.029965,
    756   0.005950, 0.004843, 0.037130,
    757   0.007588, 0.006356, 0.044973,
    758   0.009426, 0.008022, 0.052844,
    759   0.011465, 0.009828, 0.060750,
    760   0.013708, 0.011771, 0.068667,
    761   0.016156, 0.013840, 0.076603,
    762   0.018815, 0.016026, 0.084584,
    763   0.021692, 0.018320, 0.092610,
    764   0.024792, 0.020715, 0.100676,
    765   0.028123, 0.023201, 0.108787,
    766   0.031696, 0.025765, 0.116965,
    767   0.035520, 0.028397, 0.125209,
    768   0.039608, 0.031090, 0.133515,
    769   0.043830, 0.033830, 0.141886,
    770   0.048062, 0.036607, 0.150327,
    771   0.052320, 0.039407, 0.158841,
    772   0.056615, 0.042160, 0.167446,
    773   0.060949, 0.044794, 0.176129,
    774   0.065330, 0.047318, 0.184892,
    775   0.069764, 0.049726, 0.193735,
    776   0.074257, 0.052017, 0.202660,
    777   0.078815, 0.054184, 0.211667,
    778   0.083446, 0.056225, 0.220755,
    779   0.088155, 0.058133, 0.229922,
    780   0.092949, 0.059904, 0.239164,
    781   0.097833, 0.061531, 0.248477,
    782   0.102815, 0.063010, 0.257854,
    783   0.107899, 0.064335, 0.267289,
    784   0.113094, 0.065492, 0.276784,
    785   0.118405, 0.066479, 0.286321,
    786   0.123833, 0.067295, 0.295879,
    787   0.129380, 0.067935, 0.305443,
    788   0.135053, 0.068391, 0.315000,
    789   0.140858, 0.068654, 0.324538,
    790   0.146785, 0.068738, 0.334011,
    791   0.152839, 0.068637, 0.343404,
    792   0.159018, 0.068354, 0.352688,
    793   0.165308, 0.067911, 0.361816,
    794   0.171713, 0.067305, 0.370771,
    795   0.178212, 0.066576, 0.379497,
    796   0.184801, 0.065732, 0.387973,
    797   0.191460, 0.064818, 0.396152,
    798   0.198177, 0.063862, 0.404009,
    799   0.204935, 0.062907, 0.411514,
    800   0.211718, 0.061992, 0.418647,
    801   0.218512, 0.061158, 0.425392,
    802   0.225302, 0.060445, 0.431742,
    803   0.232077, 0.059889, 0.437695,
    804   0.238826, 0.059517, 0.443256,
    805   0.245543, 0.059352, 0.448436,
    806   0.252220, 0.059415, 0.453248,
    807   0.258857, 0.059706, 0.457710,
    808   0.265447, 0.060237, 0.461840,
    809   0.271994, 0.060994, 0.465660,
    810   0.278493, 0.061978, 0.469190,
    811   0.284951, 0.063168, 0.472451,
    812   0.291366, 0.064553, 0.475462,
    813   0.297740, 0.066117, 0.478243,
    814   0.304081, 0.067835, 0.480812,
    815   0.310382, 0.069702, 0.483186,
    816   0.316654, 0.071690, 0.485380,
    817   0.322899, 0.073782, 0.487408,
    818   0.329114, 0.075972, 0.489287,
    819   0.335308, 0.078236, 0.491024,
    820   0.341482, 0.080564, 0.492631,
    821   0.347636, 0.082946, 0.494121,
    822   0.353773, 0.085373, 0.495501,
    823   0.359898, 0.087831, 0.496778,
    824   0.366012, 0.090314, 0.497960,
    825   0.372116, 0.092816, 0.499053,
    826   0.378211, 0.095332, 0.500067,
    827   0.384299, 0.097855, 0.501002,
    828   0.390384, 0.100379, 0.501864,
    829   0.396467, 0.102902, 0.502658,
    830   0.402548, 0.105420, 0.503386,
    831   0.408629, 0.107930, 0.504052,
    832   0.414709, 0.110431, 0.504662,
    833   0.420791, 0.112920, 0.505215,
    834   0.426877, 0.115395, 0.505714,
    835   0.432967, 0.117855, 0.506160,
    836   0.439062, 0.120298, 0.506555,
    837   0.445163, 0.122724, 0.506901,
    838   0.451271, 0.125132, 0.507198,
    839   0.457386, 0.127522, 0.507448,
    840   0.463508, 0.129893, 0.507652,
    841   0.469640, 0.132245, 0.507809,
    842   0.475780, 0.134577, 0.507921,
    843   0.481929, 0.136891, 0.507989,
    844   0.488088, 0.139186, 0.508011,
    845   0.494258, 0.141462, 0.507988,
    846   0.500438, 0.143719, 0.507920,
    847   0.506629, 0.145958, 0.507806,
    848   0.512831, 0.148179, 0.507648,
    849   0.519045, 0.150383, 0.507443,
    850   0.525270, 0.152569, 0.507192,
    851   0.531507, 0.154739, 0.506895,
    852   0.537755, 0.156894, 0.506551,
    853   0.544015, 0.159033, 0.506159,
    854   0.550287, 0.161158, 0.505719,
    855   0.556571, 0.163269, 0.505230,
    856   0.562866, 0.165368, 0.504692,
    857   0.569172, 0.167454, 0.504105,
    858   0.575490, 0.169530, 0.503466,
    859   0.581819, 0.171596, 0.502777,
    860   0.588158, 0.173652, 0.502035,
    861   0.594508, 0.175701, 0.501241,
    862   0.600868, 0.177743, 0.500394,
    863   0.607238, 0.179779, 0.499492,
    864   0.613617, 0.181811, 0.498536,
    865   0.620005, 0.183840, 0.497524,
    866   0.626401, 0.185867, 0.496456,
    867   0.632805, 0.187893, 0.495332,
    868   0.639216, 0.189921, 0.494150,
    869   0.645633, 0.191952, 0.492910,
    870   0.652056, 0.193986, 0.491611,
    871   0.658483, 0.196027, 0.490253,
    872   0.664915, 0.198075, 0.488836,
    873   0.671349, 0.200133, 0.487358,
    874   0.677786, 0.202203, 0.485819,
    875   0.684224, 0.204286, 0.484219,
    876   0.690661, 0.206384, 0.482558,
    877   0.697098, 0.208501, 0.480835,
    878   0.703532, 0.210638, 0.479049,
    879   0.709962, 0.212797, 0.477201,
    880   0.716387, 0.214982, 0.475290,
    881   0.722805, 0.217194, 0.473316,
    882   0.729216, 0.219437, 0.471279,
    883   0.735616, 0.221713, 0.469180,
    884   0.742004, 0.224025, 0.467018,
    885   0.748378, 0.226377, 0.464794,
    886   0.754737, 0.228772, 0.462509,
    887   0.761077, 0.231214, 0.460162,
    888   0.767398, 0.233705, 0.457755,
    889   0.773695, 0.236249, 0.455289,
    890   0.779968, 0.238851, 0.452765,
    891   0.786212, 0.241514, 0.450184,
    892   0.792427, 0.244242, 0.447543,
    893   0.798608, 0.247040, 0.444848,
    894   0.804752, 0.249911, 0.442102,
    895   0.810855, 0.252861, 0.439305,
    896   0.816914, 0.255895, 0.436461,
    897   0.822926, 0.259016, 0.433573,
    898   0.828886, 0.262229, 0.430644,
    899   0.834791, 0.265540, 0.427671,
    900   0.840636, 0.268953, 0.424666,
    901   0.846416, 0.272473, 0.421631,
    902   0.852126, 0.276106, 0.418573,
    903   0.857763, 0.279857, 0.415496,
    904   0.863320, 0.283729, 0.412403,
    905   0.868793, 0.287728, 0.409303,
    906   0.874176, 0.291859, 0.406205,
    907   0.879464, 0.296125, 0.403118,
    908   0.884651, 0.300530, 0.400047,
    909   0.889731, 0.305079, 0.397002,
    910   0.894700, 0.309773, 0.393995,
    911   0.899552, 0.314616, 0.391037,
    912   0.904281, 0.319610, 0.388137,
    913   0.908884, 0.324755, 0.385308,
    914   0.913354, 0.330052, 0.382563,
    915   0.917689, 0.335500, 0.379915,
    916   0.921884, 0.341098, 0.377376,
    917   0.925937, 0.346844, 0.374959,
    918   0.929845, 0.352734, 0.372677,
    919   0.933606, 0.358764, 0.370541,
    920   0.937221, 0.364929, 0.368567,
    921   0.940687, 0.371224, 0.366762,
    922   0.944006, 0.377643, 0.365136,
    923   0.947180, 0.384178, 0.363701,
    924   0.950210, 0.390820, 0.362468,
    925   0.953099, 0.397563, 0.361438,
    926   0.955849, 0.404400, 0.360619,
    927   0.958464, 0.411324, 0.360014,
    928   0.960949, 0.418323, 0.359630,
    929   0.963310, 0.425390, 0.359469,
    930   0.965549, 0.432519, 0.359529,
    931   0.967671, 0.439703, 0.359810,
    932   0.969680, 0.446936, 0.360311,
    933   0.971582, 0.454210, 0.361030,
    934   0.973381, 0.461520, 0.361965,
    935   0.975082, 0.468861, 0.363111,
    936   0.976690, 0.476226, 0.364466,
    937   0.978210, 0.483612, 0.366025,
    938   0.979645, 0.491014, 0.367783,
    939   0.981000, 0.498428, 0.369734,
    940   0.982279, 0.505851, 0.371874,
    941   0.983485, 0.513280, 0.374198,
    942   0.984622, 0.520713, 0.376698,
    943   0.985693, 0.528148, 0.379371,
    944   0.986700, 0.535582, 0.382210,
    945   0.987646, 0.543015, 0.385210,
    946   0.988533, 0.550446, 0.388365,
    947   0.989363, 0.557873, 0.391671,
    948   0.990138, 0.565296, 0.395122,
    949   0.990871, 0.572706, 0.398714,
    950   0.991558, 0.580107, 0.402441,
    951   0.992196, 0.587502, 0.406299,
    952   0.992785, 0.594891, 0.410283,
    953   0.993326, 0.602275, 0.414390,
    954   0.993834, 0.609644, 0.418613,
    955   0.994309, 0.616999, 0.422950,
    956   0.994738, 0.624350, 0.427397,
    957   0.995122, 0.631696, 0.431951,
    958   0.995480, 0.639027, 0.436607,
    959   0.995810, 0.646344, 0.441361,
    960   0.996096, 0.653659, 0.446213,
    961   0.996341, 0.660969, 0.451160,
    962   0.996580, 0.668256, 0.456192,
    963   0.996775, 0.675541, 0.461314,
    964   0.996925, 0.682828, 0.466526,
    965   0.997077, 0.690088, 0.471811,
    966   0.997186, 0.697349, 0.477182,
    967   0.997254, 0.704611, 0.482635,
    968   0.997325, 0.711848, 0.488154,
    969   0.997351, 0.719089, 0.493755,
    970   0.997351, 0.726324, 0.499428,
    971   0.997341, 0.733545, 0.505167,
    972   0.997285, 0.740772, 0.510983,
    973   0.997228, 0.747981, 0.516859,
    974   0.997138, 0.755190, 0.522806,
    975   0.997019, 0.762398, 0.528821,
    976   0.996898, 0.769591, 0.534892,
    977   0.996727, 0.776795, 0.541039,
    978   0.996571, 0.783977, 0.547233,
    979   0.996369, 0.791167, 0.553499,
    980   0.996162, 0.798348, 0.559820,
    981   0.995932, 0.805527, 0.566202,
    982   0.995680, 0.812706, 0.572645,
    983   0.995424, 0.819875, 0.579140,
    984   0.995131, 0.827052, 0.585701,
    985   0.994851, 0.834213, 0.592307,
    986   0.994524, 0.841387, 0.598983,
    987   0.994222, 0.848540, 0.605696,
    988   0.993866, 0.855711, 0.612482,
    989   0.993545, 0.862859, 0.619299,
    990   0.993170, 0.870024, 0.626189,
    991   0.992831, 0.877168, 0.633109,
    992   0.992440, 0.884330, 0.640099,
    993   0.992089, 0.891470, 0.647116,
    994   0.991688, 0.898627, 0.654202,
    995   0.991332, 0.905763, 0.661309,
    996   0.990930, 0.912915, 0.668481,
    997   0.990570, 0.920049, 0.675675,
    998   0.990175, 0.927196, 0.682926,
    999   0.989815, 0.934329, 0.690198,
   1000   0.989434, 0.941470, 0.697519,
   1001   0.989077, 0.948604, 0.704863,
   1002   0.988717, 0.955742, 0.712242,
   1003   0.988367, 0.962878, 0.719649,
   1004   0.988033, 0.970012, 0.727077,
   1005   0.987691, 0.977154, 0.734536,
   1006   0.987387, 0.984288, 0.742002,
   1007   0.987053, 0.991438, 0.749504,
   1008 };
   1009 
   1010 static const double moreland[] = {
   1011   0.2298057, 0.298717966, 0.753683153,
   1012   0.234299935, 0.305559204, 0.759874796,
   1013   0.238810063, 0.312388385, 0.766005866,
   1014   0.243336663, 0.319205292, 0.772075394,
   1015   0.247880265, 0.326009656, 0.778082421,
   1016   0.25244136, 0.332801165, 0.784026001,
   1017   0.257020396, 0.339579464, 0.789905199,
   1018   0.261617779, 0.346344164, 0.79571909,
   1019   0.26623388, 0.353094838, 0.801466763,
   1020   0.270869029, 0.359831032, 0.807147315,
   1021   0.275523523, 0.36655226, 0.812759858,
   1022   0.28019762, 0.373258014, 0.818303516,
   1023   0.284891546, 0.379947761, 0.823777422,
   1024   0.289605495, 0.386620945, 0.829180725,
   1025   0.294339624, 0.393276993, 0.834512584,
   1026   0.299094064, 0.399915313, 0.839772171,
   1027   0.30386891, 0.406535296, 0.84495867,
   1028   0.308664231, 0.413136319, 0.850071279,
   1029   0.313480065, 0.419717745, 0.855109207,
   1030   0.318316422, 0.426278924, 0.860071679,
   1031   0.323173283, 0.432819194, 0.864957929,
   1032   0.328050603, 0.439337884, 0.869767207,
   1033   0.332948312, 0.445834313, 0.874498775,
   1034   0.337866311, 0.45230779, 0.87915191,
   1035   0.342804478, 0.458757618, 0.883725899,
   1036   0.347762667, 0.465183092, 0.888220047,
   1037   0.352740705, 0.471583499, 0.892633669,
   1038   0.357738399, 0.477958123, 0.896966095,
   1039   0.362755532, 0.484306241, 0.90121667,
   1040   0.367791863, 0.490627125, 0.905384751,
   1041   0.372847134, 0.496920043, 0.909469711,
   1042   0.37792106, 0.503184261, 0.913470934,
   1043   0.38301334, 0.50941904, 0.917387822,
   1044   0.38812365, 0.515623638, 0.921219788,
   1045   0.39325165, 0.521797312, 0.924966262,
   1046   0.398396976, 0.527939316, 0.928626686,
   1047   0.40355925, 0.534048902, 0.932200518,
   1048   0.408738074, 0.540125323, 0.93568723,
   1049   0.413933033, 0.546167829, 0.939086309,
   1050   0.419143694, 0.552175668, 0.942397257,
   1051   0.424369608, 0.558148092, 0.945619588,
   1052   0.429610311, 0.564084349, 0.948752835,
   1053   0.434865321, 0.56998369, 0.951796543,
   1054   0.440134144, 0.575845364, 0.954750272,
   1055   0.445416268, 0.581668623, 0.957613599,
   1056   0.450711169, 0.587452719, 0.960386113,
   1057   0.456018308, 0.593196905, 0.96306742,
   1058   0.461337134, 0.598900436, 0.96565714,
   1059   0.46666708, 0.604562568, 0.968154911,
   1060   0.472007569, 0.61018256, 0.970560381,
   1061   0.477358011, 0.615759672, 0.972873218,
   1062   0.482717804, 0.621293167, 0.975093102,
   1063   0.488086336, 0.626782311, 0.97721973,
   1064   0.493462982, 0.632226371, 0.979252813,
   1065   0.498847107, 0.637624618, 0.981192078,
   1066   0.504238066, 0.642976326, 0.983037268,
   1067   0.509635204, 0.648280772, 0.98478814,
   1068   0.515037856, 0.653537236, 0.986444467,
   1069   0.520445349, 0.658745003, 0.988006036,
   1070   0.525857, 0.66390336, 0.989472652,
   1071   0.531272118, 0.669011598, 0.990844132,
   1072   0.536690004, 0.674069012, 0.99212031,
   1073   0.542109949, 0.679074903, 0.993301037,
   1074   0.54753124, 0.684028574, 0.994386177,
   1075   0.552953156, 0.688929332, 0.995375608,
   1076   0.558374965, 0.693776492, 0.996269227,
   1077   0.563795935, 0.698569369, 0.997066945,
   1078   0.569215322, 0.703307287, 0.997768685,
   1079   0.574632379, 0.707989572, 0.99837439,
   1080   0.580046354, 0.712615557, 0.998884016,
   1081   0.585456486, 0.717184578, 0.999297533,
   1082   0.590862011, 0.721695979, 0.999614929,
   1083   0.596262162, 0.726149107, 0.999836203,
   1084   0.601656165, 0.730543315, 0.999961374,
   1085   0.607043242, 0.734877964, 0.999990472,
   1086   0.61242261, 0.739152418, 0.999923544,
   1087   0.617793485, 0.743366047, 0.999760652,
   1088   0.623155076, 0.747518228, 0.999501871,
   1089   0.628506592, 0.751608345, 0.999147293,
   1090   0.633847237, 0.755635786, 0.998697024,
   1091   0.639176211, 0.759599947, 0.998151185,
   1092   0.644492714, 0.763500228, 0.99750991,
   1093   0.649795942, 0.767336039, 0.996773351,
   1094   0.655085089, 0.771106793, 0.995941671,
   1095   0.660359348, 0.774811913, 0.995015049,
   1096   0.665617908, 0.778450826, 0.993993679,
   1097   0.670859959, 0.782022968, 0.992877768,
   1098   0.676084688, 0.78552778, 0.991667539,
   1099   0.681291281, 0.788964712, 0.990363227,
   1100   0.686478925, 0.792333219, 0.988965083,
   1101   0.691646803, 0.795632765, 0.987473371,
   1102   0.696794099, 0.798862821, 0.985888369,
   1103   0.701919999, 0.802022864, 0.984210369,
   1104   0.707023684, 0.805112381, 0.982439677,
   1105   0.712104339, 0.808130864, 0.980576612,
   1106   0.717161148, 0.811077814, 0.978621507,
   1107   0.722193294, 0.813952739, 0.976574709,
   1108   0.727199962, 0.816755156, 0.974436577,
   1109   0.732180337, 0.81948459, 0.972207484,
   1110   0.737133606, 0.82214057, 0.969887816,
   1111   0.742058956, 0.824722639, 0.967477972,
   1112   0.746955574, 0.827230344, 0.964978364,
   1113   0.751822652, 0.829663241, 0.962389418,
   1114   0.756659379, 0.832020895, 0.959711569,
   1115   0.761464949, 0.834302879, 0.956945269,
   1116   0.766238556, 0.836508774, 0.95409098,
   1117   0.770979397, 0.838638169, 0.951149176,
   1118   0.775686671, 0.840690662, 0.948120345,
   1119   0.780359577, 0.842665861, 0.945004985,
   1120   0.78499732, 0.84456338, 0.941803607,
   1121   0.789599105, 0.846382843, 0.938516733,
   1122   0.79416414, 0.848123884, 0.935144898,
   1123   0.798691636, 0.849786142, 0.931688648,
   1124   0.803180808, 0.85136927, 0.928148539,
   1125   0.807630872, 0.852872925, 0.92452514,
   1126   0.812041048, 0.854296776, 0.92081903,
   1127   0.81641056, 0.855640499, 0.917030798,
   1128   0.820738635, 0.856903782, 0.913161047,
   1129   0.825024503, 0.85808632, 0.909210387,
   1130   0.829267397, 0.859187816, 0.90517944,
   1131   0.833466556, 0.860207984, 0.901068838,
   1132   0.837621221, 0.861146547, 0.896879224,
   1133   0.841730637, 0.862003236, 0.892611249,
   1134   0.845794055, 0.862777795, 0.888265576,
   1135   0.849810727, 0.863469972, 0.883842876,
   1136   0.853779913, 0.864079527, 0.87934383,
   1137   0.857700874, 0.864606232, 0.874769128,
   1138   0.861572878, 0.865049863, 0.870119469,
   1139   0.865395197, 0.86541021, 0.865395561,
   1140   0.86977749, 0.863633958, 0.859948576,
   1141   0.874064226, 0.861776352, 0.854466231,
   1142   0.878255583, 0.859837644, 0.848949435,
   1143   0.882351728, 0.857818097, 0.843399101,
   1144   0.886352818, 0.85571798, 0.837816138,
   1145   0.890259, 0.853537573, 0.832201453,
   1146   0.89407041, 0.851277164, 0.826555954,
   1147   0.897787179, 0.848937047, 0.820880546,
   1148   0.901409427, 0.846517528, 0.815176131,
   1149   0.904937269, 0.844018919, 0.809443611,
   1150   0.908370816, 0.841441541, 0.803683885,
   1151   0.911710171, 0.838785722, 0.79789785,
   1152   0.914955433, 0.836051799, 0.792086401,
   1153   0.918106696, 0.833240115, 0.786250429,
   1154   0.921164054, 0.830351023, 0.780390824,
   1155   0.924127593, 0.827384882, 0.774508472,
   1156   0.926997401, 0.824342058, 0.768604257,
   1157   0.929773562, 0.821222926, 0.76267906,
   1158   0.932456159, 0.818027865, 0.756733758,
   1159   0.935045272, 0.814757264, 0.750769226,
   1160   0.937540984, 0.811411517, 0.744786333,
   1161   0.939943375, 0.807991025, 0.738785947,
   1162   0.942252526, 0.804496196, 0.732768931,
   1163   0.944468518, 0.800927443, 0.726736146,
   1164   0.946591434, 0.797285187, 0.720688446,
   1165   0.948621357, 0.793569853, 0.714626683,
   1166   0.950558373, 0.789781872, 0.708551706,
   1167   0.952402567, 0.785921682, 0.702464356,
   1168   0.954154029, 0.781989725, 0.696365473,
   1169   0.955812849, 0.777986449, 0.690255891,
   1170   0.957379123, 0.773912305, 0.68413644,
   1171   0.958852946, 0.769767752, 0.678007945,
   1172   0.960234418, 0.765553251, 0.671871226,
   1173   0.961523642, 0.761269267, 0.665727098,
   1174   0.962720725, 0.756916272, 0.659576372,
   1175   0.963825777, 0.752494738, 0.653419853,
   1176   0.964838913, 0.748005143, 0.647258341,
   1177   0.965760251, 0.743447967, 0.64109263,
   1178   0.966589914, 0.738823693, 0.634923509,
   1179   0.96732803, 0.734132809, 0.628751763,
   1180   0.967974729, 0.729375802, 0.62257817,
   1181   0.96853015, 0.724553162, 0.616403502,
   1182   0.968994435, 0.719665383, 0.610228525,
   1183   0.969367729, 0.714712956, 0.604054002,
   1184   0.969650186, 0.709696378, 0.597880686,
   1185   0.969841963, 0.704616143, 0.591709328,
   1186   0.969943224, 0.699472746, 0.585540669,
   1187   0.969954137, 0.694266682, 0.579375448,
   1188   0.969874878, 0.688998447, 0.573214394,
   1189   0.969705626, 0.683668532, 0.567058232,
   1190   0.96944657, 0.678277431, 0.560907681,
   1191   0.969097901, 0.672825633, 0.554763452,
   1192   0.968659818, 0.667313624, 0.54862625,
   1193   0.968132528, 0.661741889, 0.542496774,
   1194   0.967516241, 0.656110908, 0.536375716,
   1195   0.966811177, 0.650421156, 0.530263762,
   1196   0.966017559, 0.644673104, 0.524161591,
   1197   0.965135621, 0.638867216, 0.518069875,
   1198   0.964165599, 0.63300395, 0.511989279,
   1199   0.963107739, 0.627083758, 0.505920462,
   1200   0.961962293, 0.621107082, 0.499864075,
   1201   0.960729521, 0.615074355, 0.493820764,
   1202   0.959409687, 0.608986, 0.487791167,
   1203   0.958003065, 0.602842431, 0.481775914,
   1204   0.956509936, 0.596644046, 0.475775629,
   1205   0.954930586, 0.590391232, 0.46979093,
   1206   0.95326531, 0.584084361, 0.463822426,
   1207   0.951514411, 0.57772379, 0.457870719,
   1208   0.949678196, 0.571309856, 0.451936407,
   1209   0.947756983, 0.564842879, 0.446020077,
   1210   0.945751096, 0.558323158, 0.440122312,
   1211   0.943660866, 0.551750968, 0.434243684,
   1212   0.941486631, 0.545126562, 0.428384763,
   1213   0.939228739, 0.538450165, 0.422546107,
   1214   0.936887543, 0.531721972, 0.41672827,
   1215   0.934463404, 0.524942147, 0.410931798,
   1216   0.931956691, 0.518110821, 0.40515723,
   1217   0.929367782, 0.511228087, 0.399405096,
   1218   0.92669706, 0.504293997, 0.393675922,
   1219   0.923944917, 0.49730856, 0.387970225,
   1220   0.921111753, 0.490271735, 0.382288516,
   1221   0.918197974, 0.483183431, 0.376631297,
   1222   0.915203996, 0.476043498, 0.370999065,
   1223   0.912130241, 0.468851724, 0.36539231,
   1224   0.908977139, 0.461607831, 0.359811513,
   1225   0.905745128, 0.454311462, 0.354257151,
   1226   0.902434654, 0.446962183, 0.348729691,
   1227   0.89904617, 0.439559467, 0.343229596,
   1228   0.895580136, 0.43210269, 0.33775732,
   1229   0.892037022, 0.424591118, 0.332313313,
   1230   0.888417303, 0.417023898, 0.326898016,
   1231   0.884721464, 0.409400045, 0.321511863,
   1232   0.880949996, 0.401718425, 0.316155284,
   1233   0.877103399, 0.393977745, 0.310828702,
   1234   0.873182178, 0.386176527, 0.305532531,
   1235   0.869186849, 0.378313092, 0.300267182,
   1236   0.865117934, 0.370385535, 0.295033059,
   1237   0.860975962, 0.362391695, 0.289830559,
   1238   0.85676147, 0.354329127, 0.284660075,
   1239   0.852475004, 0.346195061, 0.279521991,
   1240   0.848117114, 0.337986361, 0.27441669,
   1241   0.843688361, 0.329699471, 0.269344545,
   1242   0.839189312, 0.32133036, 0.264305927,
   1243   0.834620542, 0.312874446, 0.259301199,
   1244   0.829982631, 0.304326513, 0.254330723,
   1245   0.82527617, 0.295680611, 0.249394851,
   1246   0.820501754, 0.286929926, 0.244493934,
   1247   0.815659988, 0.278066636, 0.239628318,
   1248   0.810751482, 0.269081721, 0.234798343,
   1249   0.805776855, 0.259964733, 0.230004348,
   1250   0.800736732, 0.250703507, 0.225246666,
   1251   0.795631745, 0.24128379, 0.220525627,
   1252   0.790462533, 0.231688768, 0.215841558,
   1253   0.785229744, 0.221898442, 0.211194782,
   1254   0.779934029, 0.211888813, 0.20658562,
   1255   0.774576051, 0.201630762, 0.202014392,
   1256   0.769156474, 0.191088518, 0.197481414,
   1257   0.763675975, 0.180217488, 0.192987001,
   1258   0.758135232, 0.168961101, 0.188531467,
   1259   0.752534934, 0.157246067, 0.184115123,
   1260   0.746875773, 0.144974956, 0.179738284,
   1261   0.741158452, 0.132014017, 0.175401259,
   1262   0.735383675, 0.1181719, 0.171104363,
   1263   0.729552157, 0.103159409, 0.166847907,
   1264   0.723664618, 0.086504694, 0.162632207,
   1265   0.717721782, 0.067344036, 0.158457578,
   1266   0.711724383, 0.043755173, 0.154324339,
   1267   0.705673158, 0.01555616, 0.150232812
   1268 };
   1269 
   1270 static const double plasma[] = {
   1271   0.050383, 0.029803, 0.527975,
   1272   0.063536, 0.028426, 0.533124,
   1273   0.075353, 0.027206, 0.538007,
   1274   0.086222, 0.026125, 0.542658,
   1275   0.096379, 0.025165, 0.547103,
   1276   0.105980, 0.024309, 0.551368,
   1277   0.115124, 0.023556, 0.555468,
   1278   0.123903, 0.022878, 0.559423,
   1279   0.132381, 0.022258, 0.563250,
   1280   0.140603, 0.021687, 0.566959,
   1281   0.148607, 0.021154, 0.570562,
   1282   0.156421, 0.020651, 0.574065,
   1283   0.164070, 0.020171, 0.577478,
   1284   0.171574, 0.019706, 0.580806,
   1285   0.178950, 0.019252, 0.584054,
   1286   0.186213, 0.018803, 0.587228,
   1287   0.193374, 0.018354, 0.590330,
   1288   0.200445, 0.017902, 0.593364,
   1289   0.207435, 0.017442, 0.596333,
   1290   0.214350, 0.016973, 0.599239,
   1291   0.221197, 0.016497, 0.602083,
   1292   0.227983, 0.016007, 0.604867,
   1293   0.234715, 0.015502, 0.607592,
   1294   0.241396, 0.014979, 0.610259,
   1295   0.248032, 0.014439, 0.612868,
   1296   0.254627, 0.013882, 0.615419,
   1297   0.261183, 0.013308, 0.617911,
   1298   0.267703, 0.012716, 0.620346,
   1299   0.274191, 0.012109, 0.622722,
   1300   0.280648, 0.011488, 0.625038,
   1301   0.287076, 0.010855, 0.627295,
   1302   0.293478, 0.010213, 0.629490,
   1303   0.299855, 0.009561, 0.631624,
   1304   0.306210, 0.008902, 0.633694,
   1305   0.312543, 0.008239, 0.635700,
   1306   0.318856, 0.007576, 0.637640,
   1307   0.325150, 0.006915, 0.639512,
   1308   0.331426, 0.006261, 0.641316,
   1309   0.337683, 0.005618, 0.643049,
   1310   0.343925, 0.004991, 0.644710,
   1311   0.350150, 0.004382, 0.646298,
   1312   0.356359, 0.003798, 0.647810,
   1313   0.362553, 0.003243, 0.649245,
   1314   0.368733, 0.002724, 0.650601,
   1315   0.374897, 0.002245, 0.651876,
   1316   0.381047, 0.001814, 0.653068,
   1317   0.387183, 0.001434, 0.654177,
   1318   0.393304, 0.001114, 0.655199,
   1319   0.399411, 0.000859, 0.656133,
   1320   0.405503, 0.000678, 0.656977,
   1321   0.411580, 0.000577, 0.657730,
   1322   0.417642, 0.000564, 0.658390,
   1323   0.423689, 0.000646, 0.658956,
   1324   0.429719, 0.000831, 0.659425,
   1325   0.435734, 0.001127, 0.659797,
   1326   0.441732, 0.001540, 0.660069,
   1327   0.447714, 0.002080, 0.660240,
   1328   0.453677, 0.002755, 0.660310,
   1329   0.459623, 0.003574, 0.660277,
   1330   0.465550, 0.004545, 0.660139,
   1331   0.471457, 0.005678, 0.659897,
   1332   0.477344, 0.006980, 0.659549,
   1333   0.483210, 0.008460, 0.659095,
   1334   0.489055, 0.010127, 0.658534,
   1335   0.494877, 0.011990, 0.657865,
   1336   0.500678, 0.014055, 0.657088,
   1337   0.506454, 0.016333, 0.656202,
   1338   0.512206, 0.018833, 0.655209,
   1339   0.517933, 0.021563, 0.654109,
   1340   0.523633, 0.024532, 0.652901,
   1341   0.529306, 0.027747, 0.651586,
   1342   0.534952, 0.031217, 0.650165,
   1343   0.540570, 0.034950, 0.648640,
   1344   0.546157, 0.038954, 0.647010,
   1345   0.551715, 0.043136, 0.645277,
   1346   0.557243, 0.047331, 0.643443,
   1347   0.562738, 0.051545, 0.641509,
   1348   0.568201, 0.055778, 0.639477,
   1349   0.573632, 0.060028, 0.637349,
   1350   0.579029, 0.064296, 0.635126,
   1351   0.584391, 0.068579, 0.632812,
   1352   0.589719, 0.072878, 0.630408,
   1353   0.595011, 0.077190, 0.627917,
   1354   0.600266, 0.081516, 0.625342,
   1355   0.605485, 0.085854, 0.622686,
   1356   0.610667, 0.090204, 0.619951,
   1357   0.615812, 0.094564, 0.617140,
   1358   0.620919, 0.098934, 0.614257,
   1359   0.625987, 0.103312, 0.611305,
   1360   0.631017, 0.107699, 0.608287,
   1361   0.636008, 0.112092, 0.605205,
   1362   0.640959, 0.116492, 0.602065,
   1363   0.645872, 0.120898, 0.598867,
   1364   0.650746, 0.125309, 0.595617,
   1365   0.655580, 0.129725, 0.592317,
   1366   0.660374, 0.134144, 0.588971,
   1367   0.665129, 0.138566, 0.585582,
   1368   0.669845, 0.142992, 0.582154,
   1369   0.674522, 0.147419, 0.578688,
   1370   0.679160, 0.151848, 0.575189,
   1371   0.683758, 0.156278, 0.571660,
   1372   0.688318, 0.160709, 0.568103,
   1373   0.692840, 0.165141, 0.564522,
   1374   0.697324, 0.169573, 0.560919,
   1375   0.701769, 0.174005, 0.557296,
   1376   0.706178, 0.178437, 0.553657,
   1377   0.710549, 0.182868, 0.550004,
   1378   0.714883, 0.187299, 0.546338,
   1379   0.719181, 0.191729, 0.542663,
   1380   0.723444, 0.196158, 0.538981,
   1381   0.727670, 0.200586, 0.535293,
   1382   0.731862, 0.205013, 0.531601,
   1383   0.736019, 0.209439, 0.527908,
   1384   0.740143, 0.213864, 0.524216,
   1385   0.744232, 0.218288, 0.520524,
   1386   0.748289, 0.222711, 0.516834,
   1387   0.752312, 0.227133, 0.513149,
   1388   0.756304, 0.231555, 0.509468,
   1389   0.760264, 0.235976, 0.505794,
   1390   0.764193, 0.240396, 0.502126,
   1391   0.768090, 0.244817, 0.498465,
   1392   0.771958, 0.249237, 0.494813,
   1393   0.775796, 0.253658, 0.491171,
   1394   0.779604, 0.258078, 0.487539,
   1395   0.783383, 0.262500, 0.483918,
   1396   0.787133, 0.266922, 0.480307,
   1397   0.790855, 0.271345, 0.476706,
   1398   0.794549, 0.275770, 0.473117,
   1399   0.798216, 0.280197, 0.469538,
   1400   0.801855, 0.284626, 0.465971,
   1401   0.805467, 0.289057, 0.462415,
   1402   0.809052, 0.293491, 0.458870,
   1403   0.812612, 0.297928, 0.455338,
   1404   0.816144, 0.302368, 0.451816,
   1405   0.819651, 0.306812, 0.448306,
   1406   0.823132, 0.311261, 0.444806,
   1407   0.826588, 0.315714, 0.441316,
   1408   0.830018, 0.320172, 0.437836,
   1409   0.833422, 0.324635, 0.434366,
   1410   0.836801, 0.329105, 0.430905,
   1411   0.840155, 0.333580, 0.427455,
   1412   0.843484, 0.338062, 0.424013,
   1413   0.846788, 0.342551, 0.420579,
   1414   0.850066, 0.347048, 0.417153,
   1415   0.853319, 0.351553, 0.413734,
   1416   0.856547, 0.356066, 0.410322,
   1417   0.859750, 0.360588, 0.406917,
   1418   0.862927, 0.365119, 0.403519,
   1419   0.866078, 0.369660, 0.400126,
   1420   0.869203, 0.374212, 0.396738,
   1421   0.872303, 0.378774, 0.393355,
   1422   0.875376, 0.383347, 0.389976,
   1423   0.878423, 0.387932, 0.386600,
   1424   0.881443, 0.392529, 0.383229,
   1425   0.884436, 0.397139, 0.379860,
   1426   0.887402, 0.401762, 0.376494,
   1427   0.890340, 0.406398, 0.373130,
   1428   0.893250, 0.411048, 0.369768,
   1429   0.896131, 0.415712, 0.366407,
   1430   0.898984, 0.420392, 0.363047,
   1431   0.901807, 0.425087, 0.359688,
   1432   0.904601, 0.429797, 0.356329,
   1433   0.907365, 0.434524, 0.352970,
   1434   0.910098, 0.439268, 0.349610,
   1435   0.912800, 0.444029, 0.346251,
   1436   0.915471, 0.448807, 0.342890,
   1437   0.918109, 0.453603, 0.339529,
   1438   0.920714, 0.458417, 0.336166,
   1439   0.923287, 0.463251, 0.332801,
   1440   0.925825, 0.468103, 0.329435,
   1441   0.928329, 0.472975, 0.326067,
   1442   0.930798, 0.477867, 0.322697,
   1443   0.933232, 0.482780, 0.319325,
   1444   0.935630, 0.487712, 0.315952,
   1445   0.937990, 0.492667, 0.312575,
   1446   0.940313, 0.497642, 0.309197,
   1447   0.942598, 0.502639, 0.305816,
   1448   0.944844, 0.507658, 0.302433,
   1449   0.947051, 0.512699, 0.299049,
   1450   0.949217, 0.517763, 0.295662,
   1451   0.951344, 0.522850, 0.292275,
   1452   0.953428, 0.527960, 0.288883,
   1453   0.955470, 0.533093, 0.285490,
   1454   0.957469, 0.538250, 0.282096,
   1455   0.959424, 0.543431, 0.278701,
   1456   0.961336, 0.548636, 0.275305,
   1457   0.963203, 0.553865, 0.271909,
   1458   0.965024, 0.559118, 0.268513,
   1459   0.966798, 0.564396, 0.265118,
   1460   0.968526, 0.569700, 0.261721,
   1461   0.970205, 0.575028, 0.258325,
   1462   0.971835, 0.580382, 0.254931,
   1463   0.973416, 0.585761, 0.251540,
   1464   0.974947, 0.591165, 0.248151,
   1465   0.976428, 0.596595, 0.244767,
   1466   0.977856, 0.602051, 0.241387,
   1467   0.979233, 0.607532, 0.238013,
   1468   0.980556, 0.613039, 0.234646,
   1469   0.981826, 0.618572, 0.231287,
   1470   0.983041, 0.624131, 0.227937,
   1471   0.984199, 0.629718, 0.224595,
   1472   0.985301, 0.635330, 0.221265,
   1473   0.986345, 0.640969, 0.217948,
   1474   0.987332, 0.646633, 0.214648,
   1475   0.988260, 0.652325, 0.211364,
   1476   0.989128, 0.658043, 0.208100,
   1477   0.989935, 0.663787, 0.204859,
   1478   0.990681, 0.669558, 0.201642,
   1479   0.991365, 0.675355, 0.198453,
   1480   0.991985, 0.681179, 0.195295,
   1481   0.992541, 0.687030, 0.192170,
   1482   0.993032, 0.692907, 0.189084,
   1483   0.993456, 0.698810, 0.186041,
   1484   0.993814, 0.704741, 0.183043,
   1485   0.994103, 0.710698, 0.180097,
   1486   0.994324, 0.716681, 0.177208,
   1487   0.994474, 0.722691, 0.174381,
   1488   0.994553, 0.728728, 0.171622,
   1489   0.994561, 0.734791, 0.168938,
   1490   0.994495, 0.740880, 0.166335,
   1491   0.994355, 0.746995, 0.163821,
   1492   0.994141, 0.753137, 0.161404,
   1493   0.993851, 0.759304, 0.159092,
   1494   0.993482, 0.765499, 0.156891,
   1495   0.993033, 0.771720, 0.154808,
   1496   0.992505, 0.777967, 0.152855,
   1497   0.991897, 0.784239, 0.151042,
   1498   0.991209, 0.790537, 0.149377,
   1499   0.990439, 0.796859, 0.147870,
   1500   0.989587, 0.803205, 0.146529,
   1501   0.988648, 0.809579, 0.145357,
   1502   0.987621, 0.815978, 0.144363,
   1503   0.986509, 0.822401, 0.143557,
   1504   0.985314, 0.828846, 0.142945,
   1505   0.984031, 0.835315, 0.142528,
   1506   0.982653, 0.841812, 0.142303,
   1507   0.981190, 0.848329, 0.142279,
   1508   0.979644, 0.854866, 0.142453,
   1509   0.977995, 0.861432, 0.142808,
   1510   0.976265, 0.868016, 0.143351,
   1511   0.974443, 0.874622, 0.144061,
   1512   0.972530, 0.881250, 0.144923,
   1513   0.970533, 0.887896, 0.145919,
   1514   0.968443, 0.894564, 0.147014,
   1515   0.966271, 0.901249, 0.148180,
   1516   0.964021, 0.907950, 0.149370,
   1517   0.961681, 0.914672, 0.150520,
   1518   0.959276, 0.921407, 0.151566,
   1519   0.956808, 0.928152, 0.152409,
   1520   0.954287, 0.934908, 0.152921,
   1521   0.951726, 0.941671, 0.152925,
   1522   0.949151, 0.948435, 0.152178,
   1523   0.946602, 0.955190, 0.150328,
   1524   0.944152, 0.961916, 0.146861,
   1525   0.941896, 0.968590, 0.140956,
   1526   0.940015, 0.975158, 0.131326,
   1527 };
   1528 
   1529 static const double viridis[] = {
   1530   0.267004, 0.004874, 0.329415,
   1531   0.268510, 0.009605, 0.335427,
   1532   0.269944, 0.014625, 0.341379,
   1533   0.271305, 0.019942, 0.347269,
   1534   0.272594, 0.025563, 0.353093,
   1535   0.273809, 0.031497, 0.358853,
   1536   0.274952, 0.037752, 0.364543,
   1537   0.276022, 0.044167, 0.370164,
   1538   0.277018, 0.050344, 0.375715,
   1539   0.277941, 0.056324, 0.381191,
   1540   0.278791, 0.062145, 0.386592,
   1541   0.279566, 0.067836, 0.391917,
   1542   0.280267, 0.073417, 0.397163,
   1543   0.280894, 0.078907, 0.402329,
   1544   0.281446, 0.084320, 0.407414,
   1545   0.281924, 0.089666, 0.412415,
   1546   0.282327, 0.094955, 0.417331,
   1547   0.282656, 0.100196, 0.422160,
   1548   0.282910, 0.105393, 0.426902,
   1549   0.283091, 0.110553, 0.431554,
   1550   0.283197, 0.115680, 0.436115,
   1551   0.283229, 0.120777, 0.440584,
   1552   0.283187, 0.125848, 0.444960,
   1553   0.283072, 0.130895, 0.449241,
   1554   0.282884, 0.135920, 0.453427,
   1555   0.282623, 0.140926, 0.457517,
   1556   0.282290, 0.145912, 0.461510,
   1557   0.281887, 0.150881, 0.465405,
   1558   0.281412, 0.155834, 0.469201,
   1559   0.280868, 0.160771, 0.472899,
   1560   0.280255, 0.165693, 0.476498,
   1561   0.279574, 0.170599, 0.479997,
   1562   0.278826, 0.175490, 0.483397,
   1563   0.278012, 0.180367, 0.486697,
   1564   0.277134, 0.185228, 0.489898,
   1565   0.276194, 0.190074, 0.493001,
   1566   0.275191, 0.194905, 0.496005,
   1567   0.274128, 0.199721, 0.498911,
   1568   0.273006, 0.204520, 0.501721,
   1569   0.271828, 0.209303, 0.504434,
   1570   0.270595, 0.214069, 0.507052,
   1571   0.269308, 0.218818, 0.509577,
   1572   0.267968, 0.223549, 0.512008,
   1573   0.266580, 0.228262, 0.514349,
   1574   0.265145, 0.232956, 0.516599,
   1575   0.263663, 0.237631, 0.518762,
   1576   0.262138, 0.242286, 0.520837,
   1577   0.260571, 0.246922, 0.522828,
   1578   0.258965, 0.251537, 0.524736,
   1579   0.257322, 0.256130, 0.526563,
   1580   0.255645, 0.260703, 0.528312,
   1581   0.253935, 0.265254, 0.529983,
   1582   0.252194, 0.269783, 0.531579,
   1583   0.250425, 0.274290, 0.533103,
   1584   0.248629, 0.278775, 0.534556,
   1585   0.246811, 0.283237, 0.535941,
   1586   0.244972, 0.287675, 0.537260,
   1587   0.243113, 0.292092, 0.538516,
   1588   0.241237, 0.296485, 0.539709,
   1589   0.239346, 0.300855, 0.540844,
   1590   0.237441, 0.305202, 0.541921,
   1591   0.235526, 0.309527, 0.542944,
   1592   0.233603, 0.313828, 0.543914,
   1593   0.231674, 0.318106, 0.544834,
   1594   0.229739, 0.322361, 0.545706,
   1595   0.227802, 0.326594, 0.546532,
   1596   0.225863, 0.330805, 0.547314,
   1597   0.223925, 0.334994, 0.548053,
   1598   0.221989, 0.339161, 0.548752,
   1599   0.220057, 0.343307, 0.549413,
   1600   0.218130, 0.347432, 0.550038,
   1601   0.216210, 0.351535, 0.550627,
   1602   0.214298, 0.355619, 0.551184,
   1603   0.212395, 0.359683, 0.551710,
   1604   0.210503, 0.363727, 0.552206,
   1605   0.208623, 0.367752, 0.552675,
   1606   0.206756, 0.371758, 0.553117,
   1607   0.204903, 0.375746, 0.553533,
   1608   0.203063, 0.379716, 0.553925,
   1609   0.201239, 0.383670, 0.554294,
   1610   0.199430, 0.387607, 0.554642,
   1611   0.197636, 0.391528, 0.554969,
   1612   0.195860, 0.395433, 0.555276,
   1613   0.194100, 0.399323, 0.555565,
   1614   0.192357, 0.403199, 0.555836,
   1615   0.190631, 0.407061, 0.556089,
   1616   0.188923, 0.410910, 0.556326,
   1617   0.187231, 0.414746, 0.556547,
   1618   0.185556, 0.418570, 0.556753,
   1619   0.183898, 0.422383, 0.556944,
   1620   0.182256, 0.426184, 0.557120,
   1621   0.180629, 0.429975, 0.557282,
   1622   0.179019, 0.433756, 0.557430,
   1623   0.177423, 0.437527, 0.557565,
   1624   0.175841, 0.441290, 0.557685,
   1625   0.174274, 0.445044, 0.557792,
   1626   0.172719, 0.448791, 0.557885,
   1627   0.171176, 0.452530, 0.557965,
   1628   0.169646, 0.456262, 0.558030,
   1629   0.168126, 0.459988, 0.558082,
   1630   0.166617, 0.463708, 0.558119,
   1631   0.165117, 0.467423, 0.558141,
   1632   0.163625, 0.471133, 0.558148,
   1633   0.162142, 0.474838, 0.558140,
   1634   0.160665, 0.478540, 0.558115,
   1635   0.159194, 0.482237, 0.558073,
   1636   0.157729, 0.485932, 0.558013,
   1637   0.156270, 0.489624, 0.557936,
   1638   0.154815, 0.493313, 0.557840,
   1639   0.153364, 0.497000, 0.557724,
   1640   0.151918, 0.500685, 0.557587,
   1641   0.150476, 0.504369, 0.557430,
   1642   0.149039, 0.508051, 0.557250,
   1643   0.147607, 0.511733, 0.557049,
   1644   0.146180, 0.515413, 0.556823,
   1645   0.144759, 0.519093, 0.556572,
   1646   0.143343, 0.522773, 0.556295,
   1647   0.141935, 0.526453, 0.555991,
   1648   0.140536, 0.530132, 0.555659,
   1649   0.139147, 0.533812, 0.555298,
   1650   0.137770, 0.537492, 0.554906,
   1651   0.136408, 0.541173, 0.554483,
   1652   0.135066, 0.544853, 0.554029,
   1653   0.133743, 0.548535, 0.553541,
   1654   0.132444, 0.552216, 0.553018,
   1655   0.131172, 0.555899, 0.552459,
   1656   0.129933, 0.559582, 0.551864,
   1657   0.128729, 0.563265, 0.551229,
   1658   0.127568, 0.566949, 0.550556,
   1659   0.126453, 0.570633, 0.549841,
   1660   0.125394, 0.574318, 0.549086,
   1661   0.124395, 0.578002, 0.548287,
   1662   0.123463, 0.581687, 0.547445,
   1663   0.122606, 0.585371, 0.546557,
   1664   0.121831, 0.589055, 0.545623,
   1665   0.121148, 0.592739, 0.544641,
   1666   0.120565, 0.596422, 0.543611,
   1667   0.120092, 0.600104, 0.542530,
   1668   0.119738, 0.603785, 0.541400,
   1669   0.119512, 0.607464, 0.540218,
   1670   0.119423, 0.611141, 0.538982,
   1671   0.119483, 0.614817, 0.537692,
   1672   0.119699, 0.618490, 0.536347,
   1673   0.120081, 0.622161, 0.534946,
   1674   0.120638, 0.625828, 0.533488,
   1675   0.121380, 0.629492, 0.531973,
   1676   0.122312, 0.633153, 0.530398,
   1677   0.123444, 0.636809, 0.528763,
   1678   0.124780, 0.640461, 0.527068,
   1679   0.126326, 0.644107, 0.525311,
   1680   0.128087, 0.647749, 0.523491,
   1681   0.130067, 0.651384, 0.521608,
   1682   0.132268, 0.655014, 0.519661,
   1683   0.134692, 0.658636, 0.517649,
   1684   0.137339, 0.662252, 0.515571,
   1685   0.140210, 0.665859, 0.513427,
   1686   0.143303, 0.669459, 0.511215,
   1687   0.146616, 0.673050, 0.508936,
   1688   0.150148, 0.676631, 0.506589,
   1689   0.153894, 0.680203, 0.504172,
   1690   0.157851, 0.683765, 0.501686,
   1691   0.162016, 0.687316, 0.499129,
   1692   0.166383, 0.690856, 0.496502,
   1693   0.170948, 0.694384, 0.493803,
   1694   0.175707, 0.697900, 0.491033,
   1695   0.180653, 0.701402, 0.488189,
   1696   0.185783, 0.704891, 0.485273,
   1697   0.191090, 0.708366, 0.482284,
   1698   0.196571, 0.711827, 0.479221,
   1699   0.202219, 0.715272, 0.476084,
   1700   0.208030, 0.718701, 0.472873,
   1701   0.214000, 0.722114, 0.469588,
   1702   0.220124, 0.725509, 0.466226,
   1703   0.226397, 0.728888, 0.462789,
   1704   0.232815, 0.732247, 0.459277,
   1705   0.239374, 0.735588, 0.455688,
   1706   0.246070, 0.738910, 0.452024,
   1707   0.252899, 0.742211, 0.448284,
   1708   0.259857, 0.745492, 0.444467,
   1709   0.266941, 0.748751, 0.440573,
   1710   0.274149, 0.751988, 0.436601,
   1711   0.281477, 0.755203, 0.432552,
   1712   0.288921, 0.758394, 0.428426,
   1713   0.296479, 0.761561, 0.424223,
   1714   0.304148, 0.764704, 0.419943,
   1715   0.311925, 0.767822, 0.415586,
   1716   0.319809, 0.770914, 0.411152,
   1717   0.327796, 0.773980, 0.406640,
   1718   0.335885, 0.777018, 0.402049,
   1719   0.344074, 0.780029, 0.397381,
   1720   0.352360, 0.783011, 0.392636,
   1721   0.360741, 0.785964, 0.387814,
   1722   0.369214, 0.788888, 0.382914,
   1723   0.377779, 0.791781, 0.377939,
   1724   0.386433, 0.794644, 0.372886,
   1725   0.395174, 0.797475, 0.367757,
   1726   0.404001, 0.800275, 0.362552,
   1727   0.412913, 0.803041, 0.357269,
   1728   0.421908, 0.805774, 0.351910,
   1729   0.430983, 0.808473, 0.346476,
   1730   0.440137, 0.811138, 0.340967,
   1731   0.449368, 0.813768, 0.335384,
   1732   0.458674, 0.816363, 0.329727,
   1733   0.468053, 0.818921, 0.323998,
   1734   0.477504, 0.821444, 0.318195,
   1735   0.487026, 0.823929, 0.312321,
   1736   0.496615, 0.826376, 0.306377,
   1737   0.506271, 0.828786, 0.300362,
   1738   0.515992, 0.831158, 0.294279,
   1739   0.525776, 0.833491, 0.288127,
   1740   0.535621, 0.835785, 0.281908,
   1741   0.545524, 0.838039, 0.275626,
   1742   0.555484, 0.840254, 0.269281,
   1743   0.565498, 0.842430, 0.262877,
   1744   0.575563, 0.844566, 0.256415,
   1745   0.585678, 0.846661, 0.249897,
   1746   0.595839, 0.848717, 0.243329,
   1747   0.606045, 0.850733, 0.236712,
   1748   0.616293, 0.852709, 0.230052,
   1749   0.626579, 0.854645, 0.223353,
   1750   0.636902, 0.856542, 0.216620,
   1751   0.647257, 0.858400, 0.209861,
   1752   0.657642, 0.860219, 0.203082,
   1753   0.668054, 0.861999, 0.196293,
   1754   0.678489, 0.863742, 0.189503,
   1755   0.688944, 0.865448, 0.182725,
   1756   0.699415, 0.867117, 0.175971,
   1757   0.709898, 0.868751, 0.169257,
   1758   0.720391, 0.870350, 0.162603,
   1759   0.730889, 0.871916, 0.156029,
   1760   0.741388, 0.873449, 0.149561,
   1761   0.751884, 0.874951, 0.143228,
   1762   0.762373, 0.876424, 0.137064,
   1763   0.772852, 0.877868, 0.131109,
   1764   0.783315, 0.879285, 0.125405,
   1765   0.793760, 0.880678, 0.120005,
   1766   0.804182, 0.882046, 0.114965,
   1767   0.814576, 0.883393, 0.110347,
   1768   0.824940, 0.884720, 0.106217,
   1769   0.835270, 0.886029, 0.102646,
   1770   0.845561, 0.887322, 0.099702,
   1771   0.855810, 0.888601, 0.097452,
   1772   0.866013, 0.889868, 0.095953,
   1773   0.876168, 0.891125, 0.095250,
   1774   0.886271, 0.892374, 0.095374,
   1775   0.896320, 0.893616, 0.096335,
   1776   0.906311, 0.894855, 0.098125,
   1777   0.916242, 0.896091, 0.100717,
   1778   0.926106, 0.897330, 0.104071,
   1779   0.935904, 0.898570, 0.108131,
   1780   0.945636, 0.899815, 0.112838,
   1781   0.955300, 0.901065, 0.118128,
   1782   0.964894, 0.902323, 0.123941,
   1783   0.974417, 0.903590, 0.130215,
   1784   0.983868, 0.904867, 0.136897,
   1785   0.993248, 0.906157, 0.143936,
   1786 };
   1787 
   1788 static INLINE void
   1789 get_color(const size_t icolor, double color[3], void* context)
   1790 {
   1791   const double* colors = context;
   1792   ASSERT(color && context);
   1793   color[0] = colors[icolor*3+0];
   1794   color[1] = colors[icolor*3+1];
   1795   color[2] = colors[icolor*3+2];
   1796 }
   1797 
   1798 #define DEFINE_PALETTE(Palette)                                                \
   1799   const struct scmap_palette scmap_palette_##Palette = {                       \
   1800     get_color, sizeof(Palette)/(sizeof(double)*3), (void*)Palette              \
   1801   }
   1802 DEFINE_PALETTE(accent);
   1803 DEFINE_PALETTE(blues);
   1804 DEFINE_PALETTE(brbg);
   1805 DEFINE_PALETTE(bugn);
   1806 DEFINE_PALETTE(bupu);
   1807 DEFINE_PALETTE(chromajs);
   1808 DEFINE_PALETTE(dark2);
   1809 DEFINE_PALETTE(gnbu);
   1810 DEFINE_PALETTE(gnpu);
   1811 DEFINE_PALETTE(greens);
   1812 DEFINE_PALETTE(greys);
   1813 DEFINE_PALETTE(inferno);
   1814 DEFINE_PALETTE(jet);
   1815 DEFINE_PALETTE(magma);
   1816 DEFINE_PALETTE(moreland);
   1817 DEFINE_PALETTE(oranges);
   1818 DEFINE_PALETTE(orrd);
   1819 DEFINE_PALETTE(paired);
   1820 DEFINE_PALETTE(parula);
   1821 DEFINE_PALETTE(pastel1);
   1822 DEFINE_PALETTE(pastel2);
   1823 DEFINE_PALETTE(piyg);
   1824 DEFINE_PALETTE(plasma);
   1825 DEFINE_PALETTE(prgn);
   1826 DEFINE_PALETTE(pubu);
   1827 DEFINE_PALETTE(pubugn);
   1828 DEFINE_PALETTE(puor);
   1829 DEFINE_PALETTE(purd);
   1830 DEFINE_PALETTE(purples);
   1831 DEFINE_PALETTE(rdbu);
   1832 DEFINE_PALETTE(rdgy);
   1833 DEFINE_PALETTE(rdpu);
   1834 DEFINE_PALETTE(rdylbu);
   1835 DEFINE_PALETTE(rdylgn);
   1836 DEFINE_PALETTE(reds);
   1837 DEFINE_PALETTE(sand);
   1838 DEFINE_PALETTE(set1);
   1839 DEFINE_PALETTE(set2);
   1840 DEFINE_PALETTE(set3);
   1841 DEFINE_PALETTE(spectral);
   1842 DEFINE_PALETTE(viridis);
   1843 DEFINE_PALETTE(whgnbu);
   1844 DEFINE_PALETTE(whylrd);
   1845 DEFINE_PALETTE(ylgn);
   1846 DEFINE_PALETTE(ylgnbu);
   1847 DEFINE_PALETTE(ylorbr);
   1848 DEFINE_PALETTE(ylorrd);
   1849 DEFINE_PALETTE(ylrd);
   1850 #undef DEFINE_PALETTE
   1851 
   1852 const struct scmap_palette*
   1853 scmap_get_builtin_palette(const char* name)
   1854 {
   1855   const struct scmap_palette* palette = NULL;
   1856   if(!strcmp(name, "accent")) palette = &scmap_palette_accent;
   1857   else if(!strcmp(name, "blues")) palette = &scmap_palette_blues;
   1858   else if(!strcmp(name, "brbg")) palette = &scmap_palette_brbg;
   1859   else if(!strcmp(name, "bugn")) palette = &scmap_palette_bugn;
   1860   else if(!strcmp(name, "bupu")) palette = &scmap_palette_bupu;
   1861   else if(!strcmp(name, "chromajs")) palette = &scmap_palette_chromajs;
   1862   else if(!strcmp(name, "dark2")) palette = &scmap_palette_dark2;
   1863   else if(!strcmp(name, "gnbu")) palette = &scmap_palette_gnbu;
   1864   else if(!strcmp(name, "gnpu")) palette = &scmap_palette_gnpu;
   1865   else if(!strcmp(name, "greens")) palette = &scmap_palette_greens;
   1866   else if(!strcmp(name, "greys")) palette = &scmap_palette_greys;
   1867   else if(!strcmp(name, "inferno")) palette = &scmap_palette_inferno;
   1868   else if(!strcmp(name, "jet")) palette = &scmap_palette_jet;
   1869   else if(!strcmp(name, "magma")) palette = &scmap_palette_magma;
   1870   else if(!strcmp(name, "moreland")) palette = &scmap_palette_moreland;
   1871   else if(!strcmp(name, "oranges")) palette = &scmap_palette_oranges;
   1872   else if(!strcmp(name, "orrd")) palette = &scmap_palette_orrd;
   1873   else if(!strcmp(name, "paired")) palette = &scmap_palette_paired;
   1874   else if(!strcmp(name, "parula")) palette = &scmap_palette_parula;
   1875   else if(!strcmp(name, "pastel1")) palette = &scmap_palette_pastel1;
   1876   else if(!strcmp(name, "pastel2")) palette = &scmap_palette_pastel2;
   1877   else if(!strcmp(name, "piyg")) palette = &scmap_palette_piyg;
   1878   else if(!strcmp(name, "plasma")) palette = &scmap_palette_plasma;
   1879   else if(!strcmp(name, "prgn")) palette = &scmap_palette_prgn;
   1880   else if(!strcmp(name, "pubu")) palette = &scmap_palette_pubu;
   1881   else if(!strcmp(name, "pubugn")) palette = &scmap_palette_pubugn;
   1882   else if(!strcmp(name, "puor")) palette = &scmap_palette_puor;
   1883   else if(!strcmp(name, "purd")) palette = &scmap_palette_purd;
   1884   else if(!strcmp(name, "purples")) palette = &scmap_palette_purples;
   1885   else if(!strcmp(name, "rdbu")) palette = &scmap_palette_rdbu;
   1886   else if(!strcmp(name, "rdgy")) palette = &scmap_palette_rdgy;
   1887   else if(!strcmp(name, "rdpu")) palette = &scmap_palette_rdpu;
   1888   else if(!strcmp(name, "rdylbu")) palette = &scmap_palette_rdylbu;
   1889   else if(!strcmp(name, "rdylgn")) palette = &scmap_palette_rdylgn;
   1890   else if(!strcmp(name, "reds")) palette = &scmap_palette_reds;
   1891   else if(!strcmp(name, "sand")) palette = &scmap_palette_sand;
   1892   else if(!strcmp(name, "set1")) palette = &scmap_palette_set1;
   1893   else if(!strcmp(name, "set2")) palette = &scmap_palette_set2;
   1894   else if(!strcmp(name, "set3")) palette = &scmap_palette_set3;
   1895   else if(!strcmp(name, "spectral")) palette = &scmap_palette_spectral;
   1896   else if(!strcmp(name, "viridis")) palette = &scmap_palette_viridis;
   1897   else if(!strcmp(name, "whgnbu")) palette = &scmap_palette_whgnbu;
   1898   else if(!strcmp(name, "whylrd")) palette = &scmap_palette_whylrd;
   1899   else if(!strcmp(name, "ylgn")) palette = &scmap_palette_ylgn;
   1900   else if(!strcmp(name, "ylgnbu")) palette = &scmap_palette_ylgnbu;
   1901   else if(!strcmp(name, "ylorbr")) palette = &scmap_palette_ylorbr;
   1902   else if(!strcmp(name, "ylorrd")) palette = &scmap_palette_ylorrd;
   1903   else if(!strcmp(name, "ylrd")) palette = &scmap_palette_ylrd;
   1904   return palette;
   1905 }
   1906