summaryrefslogtreecommitdiff
path: root/Build/source/utils/asymptote/runmath.cc
diff options
context:
space:
mode:
authorKarl Berry <karl@freefriends.org>2017-03-07 18:10:02 +0000
committerKarl Berry <karl@freefriends.org>2017-03-07 18:10:02 +0000
commita6dc131e8a2ff8ce660543a262d771380fdf1b4e (patch)
tree33ab70fd89ae0842cedcfce676aaa6d64a9eb7dc /Build/source/utils/asymptote/runmath.cc
parent78cc0e8d7d1977f9cf5c8c9ba3d9de81ed6e03a3 (diff)
asy 2.40 sources
git-svn-id: svn://tug.org/texlive/trunk@43422 c570f23f-e606-0410-a88d-b1316a301751
Diffstat (limited to 'Build/source/utils/asymptote/runmath.cc')
-rw-r--r--Build/source/utils/asymptote/runmath.cc240
1 files changed, 123 insertions, 117 deletions
diff --git a/Build/source/utils/asymptote/runmath.cc b/Build/source/utils/asymptote/runmath.cc
index e30d84676bd..14b4a980b57 100644
--- a/Build/source/utils/asymptote/runmath.cc
+++ b/Build/source/utils/asymptote/runmath.cc
@@ -63,9 +63,15 @@ function *realRealFunction();
#line 12 "runmath.in"
#include <inttypes.h>
+
#include "mathop.h"
#include "path.h"
+#ifdef __CYGWIN__
+extern "C" double yn(int, double);
+extern "C" double jn(int, double);
+#endif
+
using namespace camp;
typedef array realarray;
@@ -131,233 +137,233 @@ void Srand(Int seed)
#endif
namespace run {
-#line 75 "runmath.in"
+#line 81 "runmath.in"
// real ^(real x, Int y);
void gen_runmath0(stack *Stack)
{
Int y=vm::pop<Int>(Stack);
real x=vm::pop<real>(Stack);
-#line 76 "runmath.in"
+#line 82 "runmath.in"
{Stack->push<real>(pow(x,y)); return;}
}
-#line 80 "runmath.in"
+#line 86 "runmath.in"
// pair ^(pair z, Int y);
void gen_runmath1(stack *Stack)
{
Int y=vm::pop<Int>(Stack);
pair z=vm::pop<pair>(Stack);
-#line 81 "runmath.in"
+#line 87 "runmath.in"
{Stack->push<pair>(pow(z,y)); return;}
}
-#line 85 "runmath.in"
+#line 91 "runmath.in"
// Int quotient(Int x, Int y);
void gen_runmath2(stack *Stack)
{
Int y=vm::pop<Int>(Stack);
Int x=vm::pop<Int>(Stack);
-#line 86 "runmath.in"
+#line 92 "runmath.in"
{Stack->push<Int>(quotient<Int>()(x,y)); return;}
}
-#line 90 "runmath.in"
+#line 96 "runmath.in"
// Int abs(Int x);
void gen_runmath3(stack *Stack)
{
Int x=vm::pop<Int>(Stack);
-#line 91 "runmath.in"
+#line 97 "runmath.in"
{Stack->push<Int>(Abs(x)); return;}
}
-#line 95 "runmath.in"
+#line 101 "runmath.in"
// Int sgn(real x);
void gen_runmath4(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 96 "runmath.in"
+#line 102 "runmath.in"
{Stack->push<Int>(sgn(x)); return;}
}
-#line 100 "runmath.in"
+#line 106 "runmath.in"
// Int rand();
void gen_runmath5(stack *Stack)
{
-#line 101 "runmath.in"
+#line 107 "runmath.in"
if(initializeRandom)
Srand(1);
{Stack->push<Int>(random()); return;}
}
-#line 107 "runmath.in"
+#line 113 "runmath.in"
// void srand(Int seed);
void gen_runmath6(stack *Stack)
{
Int seed=vm::pop<Int>(Stack);
-#line 108 "runmath.in"
+#line 114 "runmath.in"
Srand(seed);
}
// a random number uniformly distributed in the interval [0,1]
-#line 113 "runmath.in"
+#line 119 "runmath.in"
// real unitrand();
void gen_runmath7(stack *Stack)
{
-#line 114 "runmath.in"
+#line 120 "runmath.in"
{Stack->push<real>(((real) random())/RANDOM_MAX); return;}
}
-#line 118 "runmath.in"
+#line 124 "runmath.in"
// Int ceil(real x);
void gen_runmath8(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 119 "runmath.in"
+#line 125 "runmath.in"
{Stack->push<Int>(Intcast(ceil(x))); return;}
}
-#line 123 "runmath.in"
+#line 129 "runmath.in"
// Int floor(real x);
void gen_runmath9(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 124 "runmath.in"
+#line 130 "runmath.in"
{Stack->push<Int>(Intcast(floor(x))); return;}
}
-#line 128 "runmath.in"
+#line 134 "runmath.in"
// Int round(real x);
void gen_runmath10(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 129 "runmath.in"
+#line 135 "runmath.in"
if(validInt(x)) {Stack->push<Int>(Round(x)); return;}
integeroverflow(0);
}
-#line 134 "runmath.in"
+#line 140 "runmath.in"
// Int Ceil(real x);
void gen_runmath11(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 135 "runmath.in"
+#line 141 "runmath.in"
{Stack->push<Int>(Ceil(x)); return;}
}
-#line 139 "runmath.in"
+#line 145 "runmath.in"
// Int Floor(real x);
void gen_runmath12(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 140 "runmath.in"
+#line 146 "runmath.in"
{Stack->push<Int>(Floor(x)); return;}
}
-#line 144 "runmath.in"
+#line 150 "runmath.in"
// Int Round(real x);
void gen_runmath13(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 145 "runmath.in"
+#line 151 "runmath.in"
{Stack->push<Int>(Round(Intcap(x))); return;}
}
-#line 149 "runmath.in"
+#line 155 "runmath.in"
// real fmod(real x, real y);
void gen_runmath14(stack *Stack)
{
real y=vm::pop<real>(Stack);
real x=vm::pop<real>(Stack);
-#line 150 "runmath.in"
+#line 156 "runmath.in"
if (y == 0.0) dividebyzero();
{Stack->push<real>(fmod(x,y)); return;}
}
-#line 155 "runmath.in"
+#line 161 "runmath.in"
// real atan2(real y, real x);
void gen_runmath15(stack *Stack)
{
real x=vm::pop<real>(Stack);
real y=vm::pop<real>(Stack);
-#line 156 "runmath.in"
+#line 162 "runmath.in"
{Stack->push<real>(atan2(y,x)); return;}
}
-#line 160 "runmath.in"
+#line 166 "runmath.in"
// real hypot(real x, real y);
void gen_runmath16(stack *Stack)
{
real y=vm::pop<real>(Stack);
real x=vm::pop<real>(Stack);
-#line 161 "runmath.in"
+#line 167 "runmath.in"
{Stack->push<real>(hypot(x,y)); return;}
}
-#line 165 "runmath.in"
+#line 171 "runmath.in"
// real remainder(real x, real y);
void gen_runmath17(stack *Stack)
{
real y=vm::pop<real>(Stack);
real x=vm::pop<real>(Stack);
-#line 166 "runmath.in"
+#line 172 "runmath.in"
{Stack->push<real>(remainder(x,y)); return;}
}
-#line 170 "runmath.in"
+#line 176 "runmath.in"
// real Jn(Int n, real x);
void gen_runmath18(stack *Stack)
{
real x=vm::pop<real>(Stack);
Int n=vm::pop<Int>(Stack);
-#line 171 "runmath.in"
+#line 177 "runmath.in"
{Stack->push<real>(jn(n,x)); return;}
}
-#line 175 "runmath.in"
+#line 181 "runmath.in"
// real Yn(Int n, real x);
void gen_runmath19(stack *Stack)
{
real x=vm::pop<real>(Stack);
Int n=vm::pop<Int>(Stack);
-#line 176 "runmath.in"
+#line 182 "runmath.in"
{Stack->push<real>(yn(n,x)); return;}
}
-#line 180 "runmath.in"
+#line 186 "runmath.in"
// real erf(real x);
void gen_runmath20(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 181 "runmath.in"
+#line 187 "runmath.in"
{Stack->push<real>(erf(x)); return;}
}
-#line 185 "runmath.in"
+#line 191 "runmath.in"
// real erfc(real x);
void gen_runmath21(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 186 "runmath.in"
+#line 192 "runmath.in"
{Stack->push<real>(erfc(x)); return;}
}
-#line 190 "runmath.in"
+#line 196 "runmath.in"
// Int factorial(Int n);
void gen_runmath22(stack *Stack)
{
Int n=vm::pop<Int>(Stack);
-#line 191 "runmath.in"
+#line 197 "runmath.in"
if(n < 0) error(invalidargument);
{Stack->push<Int>(factorial(n)); return;}
}
-#line 195 "runmath.in"
+#line 201 "runmath.in"
// Int choose(Int n, Int k);
void gen_runmath23(stack *Stack)
{
Int k=vm::pop<Int>(Stack);
Int n=vm::pop<Int>(Stack);
-#line 196 "runmath.in"
+#line 202 "runmath.in"
if(n < 0 || k < 0 || k > n) error(invalidargument);
Int f=1;
Int r=n-k;
@@ -368,12 +374,12 @@ void gen_runmath23(stack *Stack)
{Stack->push<Int>(f); return;}
}
-#line 206 "runmath.in"
+#line 212 "runmath.in"
// real gamma(real x);
void gen_runmath24(stack *Stack)
{
real x=vm::pop<real>(Stack);
-#line 207 "runmath.in"
+#line 213 "runmath.in"
#ifdef HAVE_TGAMMA
{Stack->push<real>(tgamma(x)); return;}
#else
@@ -382,14 +388,14 @@ void gen_runmath24(stack *Stack)
#endif
}
-#line 216 "runmath.in"
+#line 222 "runmath.in"
// realarray* quadraticroots(real a, real b, real c);
void gen_runmath25(stack *Stack)
{
real c=vm::pop<real>(Stack);
real b=vm::pop<real>(Stack);
real a=vm::pop<real>(Stack);
-#line 217 "runmath.in"
+#line 223 "runmath.in"
quadraticroots q(a,b,c);
array *roots=new array(q.roots);
if(q.roots >= 1) (*roots)[0]=q.t1;
@@ -397,14 +403,14 @@ void gen_runmath25(stack *Stack)
{Stack->push<realarray*>(roots); return;}
}
-#line 225 "runmath.in"
+#line 231 "runmath.in"
// pairarray* quadraticroots(explicit pair a, explicit pair b, explicit pair c);
void gen_runmath26(stack *Stack)
{
pair c=vm::pop<pair>(Stack);
pair b=vm::pop<pair>(Stack);
pair a=vm::pop<pair>(Stack);
-#line 226 "runmath.in"
+#line 232 "runmath.in"
Quadraticroots q(a,b,c);
array *roots=new array(q.roots);
if(q.roots >= 1) (*roots)[0]=q.z1;
@@ -412,7 +418,7 @@ void gen_runmath26(stack *Stack)
{Stack->push<pairarray*>(roots); return;}
}
-#line 234 "runmath.in"
+#line 240 "runmath.in"
// realarray* cubicroots(real a, real b, real c, real d);
void gen_runmath27(stack *Stack)
{
@@ -420,7 +426,7 @@ void gen_runmath27(stack *Stack)
real c=vm::pop<real>(Stack);
real b=vm::pop<real>(Stack);
real a=vm::pop<real>(Stack);
-#line 235 "runmath.in"
+#line 241 "runmath.in"
cubicroots q(a,b,c,d);
array *roots=new array(q.roots);
if(q.roots >= 1) (*roots)[0]=q.t1;
@@ -431,108 +437,108 @@ void gen_runmath27(stack *Stack)
// Logical operations
-#line 246 "runmath.in"
+#line 252 "runmath.in"
// bool !(bool b);
void gen_runmath28(stack *Stack)
{
bool b=vm::pop<bool>(Stack);
-#line 247 "runmath.in"
+#line 253 "runmath.in"
{Stack->push<bool>(!b); return;}
}
-#line 252 "runmath.in"
+#line 258 "runmath.in"
void boolMemEq(stack *Stack)
{
frame * b=vm::pop<frame *>(Stack);
frame * a=vm::pop<frame *>(Stack);
-#line 253 "runmath.in"
+#line 259 "runmath.in"
{Stack->push<bool>(a == b); return;}
}
-#line 257 "runmath.in"
+#line 263 "runmath.in"
void boolMemNeq(stack *Stack)
{
frame * b=vm::pop<frame *>(Stack);
frame * a=vm::pop<frame *>(Stack);
-#line 258 "runmath.in"
+#line 264 "runmath.in"
{Stack->push<bool>(a != b); return;}
}
-#line 262 "runmath.in"
+#line 268 "runmath.in"
void boolFuncEq(stack *Stack)
{
callable * b=vm::pop<callable *>(Stack);
callable * a=vm::pop<callable *>(Stack);
-#line 263 "runmath.in"
+#line 269 "runmath.in"
{Stack->push<bool>(a->compare(b)); return;}
}
-#line 267 "runmath.in"
+#line 273 "runmath.in"
void boolFuncNeq(stack *Stack)
{
callable * b=vm::pop<callable *>(Stack);
callable * a=vm::pop<callable *>(Stack);
-#line 268 "runmath.in"
+#line 274 "runmath.in"
{Stack->push<bool>(!(a->compare(b))); return;}
}
// Bit operations
-#line 274 "runmath.in"
+#line 280 "runmath.in"
// Int AND(Int a, Int b);
void gen_runmath33(stack *Stack)
{
Int b=vm::pop<Int>(Stack);
Int a=vm::pop<Int>(Stack);
-#line 275 "runmath.in"
+#line 281 "runmath.in"
{Stack->push<Int>(a & b); return;}
}
-#line 280 "runmath.in"
+#line 286 "runmath.in"
// Int OR(Int a, Int b);
void gen_runmath34(stack *Stack)
{
Int b=vm::pop<Int>(Stack);
Int a=vm::pop<Int>(Stack);
-#line 281 "runmath.in"
+#line 287 "runmath.in"
{Stack->push<Int>(a | b); return;}
}
-#line 285 "runmath.in"
+#line 291 "runmath.in"
// Int XOR(Int a, Int b);
void gen_runmath35(stack *Stack)
{
Int b=vm::pop<Int>(Stack);
Int a=vm::pop<Int>(Stack);
-#line 286 "runmath.in"
+#line 292 "runmath.in"
{Stack->push<Int>(a ^ b); return;}
}
-#line 290 "runmath.in"
+#line 296 "runmath.in"
// Int NOT(Int a);
void gen_runmath36(stack *Stack)
{
Int a=vm::pop<Int>(Stack);
-#line 291 "runmath.in"
+#line 297 "runmath.in"
{Stack->push<Int>(~a); return;}
}
-#line 295 "runmath.in"
+#line 301 "runmath.in"
// Int CLZ(Int a);
void gen_runmath37(stack *Stack)
{
Int a=vm::pop<Int>(Stack);
-#line 296 "runmath.in"
+#line 302 "runmath.in"
if((uint32_t) a > 0xFFFFFFFF) {Stack->push<Int>(-1); return;}
{Stack->push<Int>(CLZ((uint32_t) a)); return;}
}
-#line 301 "runmath.in"
+#line 307 "runmath.in"
// Int CTZ(Int a);
void gen_runmath38(stack *Stack)
{
Int a=vm::pop<Int>(Stack);
-#line 302 "runmath.in"
+#line 308 "runmath.in"
if((uint32_t) a > 0xFFFFFFFF) {Stack->push<Int>(-1); return;}
#if __GNUC__
{Stack->push<Int>(__builtin_ctz(a)); return;}
@@ -553,83 +559,83 @@ namespace trans {
void gen_runmath_venv(venv &ve)
{
-#line 75 "runmath.in"
+#line 81 "runmath.in"
addFunc(ve, run::gen_runmath0, primReal(), SYM_CARET, formal(primReal(), SYM(x), false, false), formal(primInt(), SYM(y), false, false));
-#line 80 "runmath.in"
+#line 86 "runmath.in"
addFunc(ve, run::gen_runmath1, primPair(), SYM_CARET, formal(primPair(), SYM(z), false, false), formal(primInt(), SYM(y), false, false));
-#line 85 "runmath.in"
+#line 91 "runmath.in"
addFunc(ve, run::gen_runmath2, primInt(), SYM(quotient), formal(primInt(), SYM(x), false, false), formal(primInt(), SYM(y), false, false));
-#line 90 "runmath.in"
+#line 96 "runmath.in"
addFunc(ve, run::gen_runmath3, primInt(), SYM(abs), formal(primInt(), SYM(x), false, false));
-#line 95 "runmath.in"
+#line 101 "runmath.in"
addFunc(ve, run::gen_runmath4, primInt(), SYM(sgn), formal(primReal(), SYM(x), false, false));
-#line 100 "runmath.in"
+#line 106 "runmath.in"
addFunc(ve, run::gen_runmath5, primInt(), SYM(rand));
-#line 107 "runmath.in"
+#line 113 "runmath.in"
addFunc(ve, run::gen_runmath6, primVoid(), SYM(srand), formal(primInt(), SYM(seed), false, false));
-#line 112 "runmath.in"
- addFunc(ve, run::gen_runmath7, primReal(), SYM(unitrand));
#line 118 "runmath.in"
+ addFunc(ve, run::gen_runmath7, primReal(), SYM(unitrand));
+#line 124 "runmath.in"
addFunc(ve, run::gen_runmath8, primInt(), SYM(ceil), formal(primReal(), SYM(x), false, false));
-#line 123 "runmath.in"
+#line 129 "runmath.in"
addFunc(ve, run::gen_runmath9, primInt(), SYM(floor), formal(primReal(), SYM(x), false, false));
-#line 128 "runmath.in"
- addFunc(ve, run::gen_runmath10, primInt(), SYM(round), formal(primReal(), SYM(x), false, false));
#line 134 "runmath.in"
+ addFunc(ve, run::gen_runmath10, primInt(), SYM(round), formal(primReal(), SYM(x), false, false));
+#line 140 "runmath.in"
addFunc(ve, run::gen_runmath11, primInt(), SYM(Ceil), formal(primReal(), SYM(x), false, false));
-#line 139 "runmath.in"
+#line 145 "runmath.in"
addFunc(ve, run::gen_runmath12, primInt(), SYM(Floor), formal(primReal(), SYM(x), false, false));
-#line 144 "runmath.in"
+#line 150 "runmath.in"
addFunc(ve, run::gen_runmath13, primInt(), SYM(Round), formal(primReal(), SYM(x), false, false));
-#line 149 "runmath.in"
- addFunc(ve, run::gen_runmath14, primReal(), SYM(fmod), formal(primReal(), SYM(x), false, false), formal(primReal(), SYM(y), false, false));
#line 155 "runmath.in"
+ addFunc(ve, run::gen_runmath14, primReal(), SYM(fmod), formal(primReal(), SYM(x), false, false), formal(primReal(), SYM(y), false, false));
+#line 161 "runmath.in"
addFunc(ve, run::gen_runmath15, primReal(), SYM(atan2), formal(primReal(), SYM(y), false, false), formal(primReal(), SYM(x), false, false));
-#line 160 "runmath.in"
+#line 166 "runmath.in"
addFunc(ve, run::gen_runmath16, primReal(), SYM(hypot), formal(primReal(), SYM(x), false, false), formal(primReal(), SYM(y), false, false));
-#line 165 "runmath.in"
+#line 171 "runmath.in"
addFunc(ve, run::gen_runmath17, primReal(), SYM(remainder), formal(primReal(), SYM(x), false, false), formal(primReal(), SYM(y), false, false));
-#line 170 "runmath.in"
+#line 176 "runmath.in"
addFunc(ve, run::gen_runmath18, primReal(), SYM(Jn), formal(primInt(), SYM(n), false, false), formal(primReal(), SYM(x), false, false));
-#line 175 "runmath.in"
+#line 181 "runmath.in"
addFunc(ve, run::gen_runmath19, primReal(), SYM(Yn), formal(primInt(), SYM(n), false, false), formal(primReal(), SYM(x), false, false));
-#line 180 "runmath.in"
+#line 186 "runmath.in"
addFunc(ve, run::gen_runmath20, primReal(), SYM(erf), formal(primReal(), SYM(x), false, false));
-#line 185 "runmath.in"
+#line 191 "runmath.in"
addFunc(ve, run::gen_runmath21, primReal(), SYM(erfc), formal(primReal(), SYM(x), false, false));
-#line 190 "runmath.in"
+#line 196 "runmath.in"
addFunc(ve, run::gen_runmath22, primInt(), SYM(factorial), formal(primInt(), SYM(n), false, false));
-#line 195 "runmath.in"
+#line 201 "runmath.in"
addFunc(ve, run::gen_runmath23, primInt(), SYM(choose), formal(primInt(), SYM(n), false, false), formal(primInt(), SYM(k), false, false));
-#line 206 "runmath.in"
+#line 212 "runmath.in"
addFunc(ve, run::gen_runmath24, primReal(), SYM(gamma), formal(primReal(), SYM(x), false, false));
-#line 216 "runmath.in"
+#line 222 "runmath.in"
addFunc(ve, run::gen_runmath25, realArray(), SYM(quadraticroots), formal(primReal(), SYM(a), false, false), formal(primReal(), SYM(b), false, false), formal(primReal(), SYM(c), false, false));
-#line 225 "runmath.in"
+#line 231 "runmath.in"
addFunc(ve, run::gen_runmath26, pairArray(), SYM(quadraticroots), formal(primPair(), SYM(a), false, true), formal(primPair(), SYM(b), false, true), formal(primPair(), SYM(c), false, true));
-#line 234 "runmath.in"
+#line 240 "runmath.in"
addFunc(ve, run::gen_runmath27, realArray(), SYM(cubicroots), formal(primReal(), SYM(a), false, false), formal(primReal(), SYM(b), false, false), formal(primReal(), SYM(c), false, false), formal(primReal(), SYM(d), false, false));
-#line 244 "runmath.in"
+#line 250 "runmath.in"
addFunc(ve, run::gen_runmath28, primBoolean(), SYM_LOGNOT, formal(primBoolean(), SYM(b), false, false));
-#line 252 "runmath.in"
+#line 258 "runmath.in"
REGISTER_BLTIN(run::boolMemEq,"boolMemEq");
-#line 257 "runmath.in"
+#line 263 "runmath.in"
REGISTER_BLTIN(run::boolMemNeq,"boolMemNeq");
-#line 262 "runmath.in"
+#line 268 "runmath.in"
REGISTER_BLTIN(run::boolFuncEq,"boolFuncEq");
-#line 267 "runmath.in"
+#line 273 "runmath.in"
REGISTER_BLTIN(run::boolFuncNeq,"boolFuncNeq");
-#line 272 "runmath.in"
+#line 278 "runmath.in"
addFunc(ve, run::gen_runmath33, primInt(), SYM(AND), formal(primInt(), SYM(a), false, false), formal(primInt(), SYM(b), false, false));
-#line 280 "runmath.in"
+#line 286 "runmath.in"
addFunc(ve, run::gen_runmath34, primInt(), SYM(OR), formal(primInt(), SYM(a), false, false), formal(primInt(), SYM(b), false, false));
-#line 285 "runmath.in"
+#line 291 "runmath.in"
addFunc(ve, run::gen_runmath35, primInt(), SYM(XOR), formal(primInt(), SYM(a), false, false), formal(primInt(), SYM(b), false, false));
-#line 290 "runmath.in"
+#line 296 "runmath.in"
addFunc(ve, run::gen_runmath36, primInt(), SYM(NOT), formal(primInt(), SYM(a), false, false));
-#line 295 "runmath.in"
- addFunc(ve, run::gen_runmath37, primInt(), SYM(CLZ), formal(primInt(), SYM(a), false, false));
#line 301 "runmath.in"
+ addFunc(ve, run::gen_runmath37, primInt(), SYM(CLZ), formal(primInt(), SYM(a), false, false));
+#line 307 "runmath.in"
addFunc(ve, run::gen_runmath38, primInt(), SYM(CTZ), formal(primInt(), SYM(a), false, false));
}