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