00001
00002 #include <basix/double.hpp>
00003 #include <basix/int.hpp>
00004 #include <basix/vector.hpp>
00005 #include <basix/port.hpp>
00006 #include <basix/literal.hpp>
00007 #include <numerix/integer.hpp>
00008 #include <numerix/modular.hpp>
00009 #include <numerix/modular_integer.hpp>
00010 #include <numerix/rational.hpp>
00011 #include <numerix/complex.hpp>
00012 #include <numerix/complex_double.hpp>
00013 #include <basix/glue.hpp>
00014
00015 #define double_literal(x) as_double (as_string (x))
00016 #define int_literal(x) as_int (as_string (x))
00017 #define is_generic_literal is<literal>
00018 #define gen_literal_apply(f,v) gen (as<generic> (f), v)
00019 #define gen_literal_access(f,v) access (as<generic> (f), v)
00020
00021 namespace mmx {
00022 static complex<double>
00023 GLUE_1 (const double &arg_1) {
00024 return complex<double > (arg_1);
00025 }
00026
00027 static complex<double>
00028 GLUE_2 (const double &arg_1, const double &arg_2) {
00029 return complex<double > (arg_1, arg_2);
00030 }
00031
00032 static complex<double>
00033 GLUE_3 (const double &arg_1) {
00034 return complex<double > (arg_1);
00035 }
00036
00037 static double
00038 GLUE_4 (const complex<double> &arg_1) {
00039 return Re (arg_1);
00040 }
00041
00042 static double
00043 GLUE_5 (const complex<double> &arg_1) {
00044 return Im (arg_1);
00045 }
00046
00047 static complex<double>
00048 GLUE_6 (const complex<double> &arg_1) {
00049 return conj (arg_1);
00050 }
00051
00052 static complex<double>
00053 GLUE_7 (const complex<double> &arg_1) {
00054 return times_i (arg_1);
00055 }
00056
00057 static complex<double>
00058 GLUE_8 (const complex<double> &arg_1) {
00059 return over_i (arg_1);
00060 }
00061
00062 static complex<double>
00063 GLUE_9 (const complex<double> &arg_1) {
00064 return -arg_1;
00065 }
00066
00067 static complex<double>
00068 GLUE_10 (const complex<double> &arg_1) {
00069 return square (arg_1);
00070 }
00071
00072 static complex<double>
00073 GLUE_11 (const complex<double> &arg_1, const complex<double> &arg_2) {
00074 return arg_1 + arg_2;
00075 }
00076
00077 static complex<double>
00078 GLUE_12 (const complex<double> &arg_1, const complex<double> &arg_2) {
00079 return arg_1 - arg_2;
00080 }
00081
00082 static complex<double>
00083 GLUE_13 (const complex<double> &arg_1, const complex<double> &arg_2) {
00084 return arg_1 * arg_2;
00085 }
00086
00087 static complex<double>
00088 GLUE_14 (const double &arg_1, const complex<double> &arg_2) {
00089 return arg_1 + arg_2;
00090 }
00091
00092 static complex<double>
00093 GLUE_15 (const complex<double> &arg_1, const double &arg_2) {
00094 return arg_1 + arg_2;
00095 }
00096
00097 static complex<double>
00098 GLUE_16 (const double &arg_1, const complex<double> &arg_2) {
00099 return arg_1 - arg_2;
00100 }
00101
00102 static complex<double>
00103 GLUE_17 (const complex<double> &arg_1, const double &arg_2) {
00104 return arg_1 - arg_2;
00105 }
00106
00107 static complex<double>
00108 GLUE_18 (const double &arg_1, const complex<double> &arg_2) {
00109 return arg_1 * arg_2;
00110 }
00111
00112 static complex<double>
00113 GLUE_19 (const complex<double> &arg_1, const double &arg_2) {
00114 return arg_1 * arg_2;
00115 }
00116
00117 static bool
00118 GLUE_20 (const complex<double> &arg_1, const complex<double> &arg_2) {
00119 return arg_1 == arg_2;
00120 }
00121
00122 static bool
00123 GLUE_21 (const complex<double> &arg_1, const complex<double> &arg_2) {
00124 return arg_1 != arg_2;
00125 }
00126
00127 static bool
00128 GLUE_22 (const complex<double> &arg_1, const double &arg_2) {
00129 return arg_1 == arg_2;
00130 }
00131
00132 static bool
00133 GLUE_23 (const complex<double> &arg_1, const double &arg_2) {
00134 return arg_1 != arg_2;
00135 }
00136
00137 static bool
00138 GLUE_24 (const double &arg_1, const complex<double> &arg_2) {
00139 return arg_1 == arg_2;
00140 }
00141
00142 static bool
00143 GLUE_25 (const double &arg_1, const complex<double> &arg_2) {
00144 return arg_1 != arg_2;
00145 }
00146
00147 static complex<double>
00148 GLUE_26 (const complex<double> &arg_1) {
00149 return invert (arg_1);
00150 }
00151
00152 static complex<double>
00153 GLUE_27 (const complex<double> &arg_1, const complex<double> &arg_2) {
00154 return arg_1 / arg_2;
00155 }
00156
00157 static complex<double>
00158 GLUE_28 (const double &arg_1, const complex<double> &arg_2) {
00159 return arg_1 / arg_2;
00160 }
00161
00162 static complex<double>
00163 GLUE_29 (const complex<double> &arg_1, const double &arg_2) {
00164 return arg_1 / arg_2;
00165 }
00166
00167 static double
00168 GLUE_30 (const complex<double> &arg_1) {
00169 return abs (arg_1);
00170 }
00171
00172 static double
00173 GLUE_31 (const complex<double> &arg_1) {
00174 return arg (arg_1);
00175 }
00176
00177 static complex<double>
00178 GLUE_32 (const complex<double> &arg_1) {
00179 return sqrt (arg_1);
00180 }
00181
00182 static complex<double>
00183 GLUE_33 (const complex<double> &arg_1) {
00184 return exp (arg_1);
00185 }
00186
00187 static complex<double>
00188 GLUE_34 (const complex<double> &arg_1) {
00189 return log (arg_1);
00190 }
00191
00192 static complex<double>
00193 GLUE_35 (const complex<double> &arg_1) {
00194 return cos (arg_1);
00195 }
00196
00197 static complex<double>
00198 GLUE_36 (const complex<double> &arg_1) {
00199 return sin (arg_1);
00200 }
00201
00202 static complex<double>
00203 GLUE_37 (const complex<double> &arg_1) {
00204 return tan (arg_1);
00205 }
00206
00207 static complex<double>
00208 GLUE_38 (const complex<double> &arg_1) {
00209 return acos (arg_1);
00210 }
00211
00212 static complex<double>
00213 GLUE_39 (const complex<double> &arg_1) {
00214 return asin (arg_1);
00215 }
00216
00217 static complex<double>
00218 GLUE_40 (const complex<double> &arg_1) {
00219 return atan (arg_1);
00220 }
00221
00222 static complex<double>
00223 GLUE_41 (const complex<double> &arg_1, const complex<double> &arg_2) {
00224 return pow (arg_1, arg_2);
00225 }
00226
00227 static complex<double>
00228 GLUE_42 (const complex<double> &arg_1, const double &arg_2) {
00229 return pow (arg_1, arg_2);
00230 }
00231
00232 static complex<double>
00233 GLUE_43 (const double &arg_1, const complex<double> &arg_2) {
00234 return pow (arg_1, arg_2);
00235 }
00236
00237 static complex<double>
00238 GLUE_44 (const complex<double> &arg_1, const int &arg_2) {
00239 return sympow (arg_1, arg_2);
00240 }
00241
00242 static bool
00243 GLUE_45 (const complex<double> &arg_1) {
00244 return is_finite (arg_1);
00245 }
00246
00247 static bool
00248 GLUE_46 (const complex<double> &arg_1) {
00249 return is_infinite (arg_1);
00250 }
00251
00252 static bool
00253 GLUE_47 (const complex<double> &arg_1) {
00254 return is_nan (arg_1);
00255 }
00256
00257 static complex<double>
00258 GLUE_48 (const complex<double> &arg_1) {
00259 return times_infinity (arg_1);
00260 }
00261
00262 static int
00263 GLUE_49 (const complex<double> &arg_1) {
00264 return precision (arg_1);
00265 }
00266
00267 static complex<double>
00268 GLUE_50 (const complex<double> &arg_1, const int &arg_2) {
00269 return change_precision (arg_1, arg_2);
00270 }
00271
00272 static int
00273 GLUE_51 (const complex<double> &arg_1) {
00274 return exponent (arg_1);
00275 }
00276
00277 static double
00278 GLUE_52 (const complex<double> &arg_1) {
00279 return magnitude (arg_1);
00280 }
00281
00282 static complex<double>
00283 GLUE_53 (const complex<double> &arg_1, const int &arg_2) {
00284 return incexp2 (arg_1, arg_2);
00285 }
00286
00287 static complex<double>
00288 GLUE_54 (const complex<double> &arg_1, const int &arg_2) {
00289 return decexp2 (arg_1, arg_2);
00290 }
00291
00292 static double
00293 GLUE_55 (const complex<double> &arg_1) {
00294 return rounding_error (arg_1);
00295 }
00296
00297 static double
00298 GLUE_56 (const complex<double> &arg_1) {
00299 return additive_error (arg_1);
00300 }
00301
00302 static double
00303 GLUE_57 (const complex<double> &arg_1) {
00304 return multiplicative_error (arg_1);
00305 }
00306
00307 static double
00308 GLUE_58 (const complex<double> &arg_1) {
00309 return elementary_error (arg_1);
00310 }
00311
00312 static complex<generic>
00313 GLUE_59 (const complex<double> &arg_1) {
00314 return as<complex<generic> > (arg_1);
00315 }
00316
00317 static complex<double>
00318 GLUE_60 (const complex<rational> &arg_1) {
00319 return as<complex<double> > (arg_1);
00320 }
00321
00322 void
00323 glue_complex_double () {
00324 static bool done = false;
00325 if (done) return;
00326 done = true;
00327 call_glue (string ("glue_double"));
00328 call_glue (string ("glue_complex_rational"));
00329 define_type<complex<double> > (gen (lit ("Complex"), lit ("Double")));
00330 define ("complex", GLUE_1);
00331 define ("complex", GLUE_2);
00332 define_converter ("upgrade", GLUE_3, PENALTY_HOMOMORPHISM);
00333 define ("Re", GLUE_4);
00334 define ("Im", GLUE_5);
00335 define ("conj", GLUE_6);
00336 define ("times_i", GLUE_7);
00337 define ("over_i", GLUE_8);
00338 define ("-", GLUE_9);
00339 define ("square", GLUE_10);
00340 define ("+", GLUE_11);
00341 define ("-", GLUE_12);
00342 define ("*", GLUE_13);
00343 define ("+", GLUE_14);
00344 define ("+", GLUE_15);
00345 define ("-", GLUE_16);
00346 define ("-", GLUE_17);
00347 define ("*", GLUE_18);
00348 define ("*", GLUE_19);
00349 define ("=", GLUE_20);
00350 define ("!=", GLUE_21);
00351 define ("=", GLUE_22);
00352 define ("!=", GLUE_23);
00353 define ("=", GLUE_24);
00354 define ("!=", GLUE_25);
00355 define ("invert", GLUE_26);
00356 define ("/", GLUE_27);
00357 define ("/", GLUE_28);
00358 define ("/", GLUE_29);
00359 define ("abs", GLUE_30);
00360 define ("arg", GLUE_31);
00361 define ("sqrt", GLUE_32);
00362 define ("exp", GLUE_33);
00363 define ("log", GLUE_34);
00364 define ("cos", GLUE_35);
00365 define ("sin", GLUE_36);
00366 define ("tan", GLUE_37);
00367 define ("arccos", GLUE_38);
00368 define ("arcsin", GLUE_39);
00369 define ("arctan", GLUE_40);
00370 define ("^", GLUE_41);
00371 define ("^", GLUE_42);
00372 define ("^", GLUE_43);
00373 define ("sym_pow", GLUE_44);
00374 define ("finite?", GLUE_45);
00375 define ("infinite?", GLUE_46);
00376 define ("nan?", GLUE_47);
00377 define ("times_infinity", GLUE_48);
00378 define ("precision", GLUE_49);
00379 define ("change_precision", GLUE_50);
00380 define ("exponent", GLUE_51);
00381 define ("magnitude", GLUE_52);
00382 define ("increase_exponent", GLUE_53);
00383 define ("decrease_exponent", GLUE_54);
00384 define ("rounding_error", GLUE_55);
00385 define ("additive_error", GLUE_56);
00386 define ("multiplicative_error", GLUE_57);
00387 define ("elementary_error", GLUE_58);
00388 define_converter (":>", GLUE_59, PENALTY_PROMOTE_GENERIC);
00389 define_converter (":>", GLUE_60, PENALTY_INCLUSION);
00390 }
00391 }
