Annotation of parser3/src/classes/math.C, revision 1.41
1.1 parser 1: /** @file
2: Parser: @b math parser class.
3:
1.40 paf 4: Copyright(c) 2001-2005 ArtLebedev Group(http://www.artlebedev.com)
1.20 paf 5: Author: Alexandr Petrosian <paf@design.ru>(http://paf.design.ru)
1.28 paf 6:
7: portions from gen_uuid.c,
8: Copyright (C) 1996, 1997, 1998, 1999 Theodore Ts'o.
1.22 paf 9: */
1.1 parser 10:
1.41 ! paf 11: static const char * const IDENT_MATH_C="$Date: 2005/08/09 08:14:48 $";
1.1 parser 12:
1.34 paf 13: #include "pa_vmethod_frame.h"
1.1 parser 14: #include "pa_common.h"
15: #include "pa_vint.h"
1.2 parser 16: #include "pa_vmath.h"
1.1 parser 17: #include "pa_request.h"
1.19 paf 18: #include "pa_md5.h"
1.28 paf 19: #include "pa_threads.h"
1.1 parser 20:
1.5 parser 21: #ifdef WIN32
1.28 paf 22: # define _WIN32_WINNT 0x400
1.5 parser 23: # include <windows.h>
1.28 paf 24: # include <wincrypt.h>
1.5 parser 25: #endif
26:
1.34 paf 27: #ifdef HAVE_CRYPT
28: # ifdef HAVE_CRYPT_H
29: # include <crypt.h>
30: # endif
31: #else
32: extern char *crypt(const char* , const char* );
1.20 paf 33: #endif
34:
1.1 parser 35: // defines
36:
1.20 paf 37: #define MAX_SALT 8
1.1 parser 38:
39: // class
40:
1.34 paf 41: class MMath: public Methoded {
1.1 parser 42: public:
1.34 paf 43: MMath();
1.5 parser 44:
1.1 parser 45: public: // Methoded
1.15 paf 46: bool used_directly() { return false; }
1.1 parser 47: };
48:
1.34 paf 49: // global variables
50:
51: DECLARE_CLASS_VAR(math, 0 /*fictive*/, new MMath);
52:
1.28 paf 53: #ifdef WIN32
54: class Random_provider {
55: HCRYPTPROV fhProv;
56:
57: void acquire() {
58: SYNCHRONIZED;
59:
60: if(fhProv)
61: return;
62:
63: if(!CryptAcquireContext(&fhProv, NULL, NULL, PROV_RSA_FULL, CRYPT_VERIFYCONTEXT))
64: throw Exception(0,
65: 0,
66: "CryptAcquireContext failed");
67: }
68: void release() {
69: if(fhProv)
70: CryptReleaseContext(fhProv, 0);
71: }
72:
73: public:
74: Random_provider(): fhProv(0) {}
75: ~Random_provider() { release(); }
76: void generate(void *buffer, size_t size) {
77: acquire();
78:
79: if(!CryptGenRandom(fhProv, size, (BYTE*)buffer))
80: throw Exception(0,
81: 0,
82: "CryptGenRandom failed");
83: }
84: }
85: random_provider;
86:
87: #else
88:
89: /// from gen_uuid.c
90: static int get_random_fd(void)
91: {
92: struct timeval tv;
93: static int fd = -2;
94: int i;
95:
96: if (fd == -2) {
97: gettimeofday(&tv, 0);
98: fd = open("/dev/urandom", O_RDONLY);
99: if (fd == -1)
100: fd = open("/dev/random", O_RDONLY | O_NONBLOCK);
101: srand((getpid() << 16) ^ getuid() ^ tv.tv_sec ^ tv.tv_usec);
102: }
103: /* Crank the random number generator a few times */
104: gettimeofday(&tv, 0);
105: for (i = (tv.tv_sec ^ tv.tv_usec) & 0x1F; i > 0; i--)
106: rand();
107: return fd;
108: }
109:
110:
111: /*
112: * Generate a series of random bytes. Use /dev/urandom if possible,
113: * and if not, use srandom/random.
114: */
115: static void get_random_bytes(void *buf, int nbytes)
116: {
117: int i, fd = get_random_fd();
118: int lose_counter = 0;
119: char *cp = (char *) buf;
120:
121: if (fd >= 0) {
122: while (nbytes > 0) {
123: i = read(fd, cp, nbytes);
124: if (i <= 0) {
125: if (lose_counter++ > 16)
126: break;
127: continue;
128: }
129: nbytes -= i;
130: cp += i;
131: lose_counter = 0;
132: }
133: }
134:
135: /* XXX put something better here if no /dev/random! */
136: for (i = 0; i < nbytes; i++)
137: *cp++ = rand() & 0xFF;
138: return;
139: }
140:
141:
142: #endif
143:
144:
145: // helpers
146:
147: static void random(void *buffer, size_t size) {
148: #ifdef WIN32
149: random_provider.generate(buffer, size);
150: #else
151: get_random_bytes(buffer, size);
152: #endif
153: }
154:
1.34 paf 155:
1.1 parser 156: // methods
1.34 paf 157:
158: #define MAX_UINT 0xFFFFFFFFu
159:
1.20 paf 160: static inline int _random(uint top) {
1.28 paf 161: uint raw;
162: random(&raw, sizeof(raw));
1.34 paf 163: return int(double(raw) / MAX_UINT * top );
1.20 paf 164: }
1.1 parser 165:
1.34 paf 166: static void _random(Request& r, MethodParams& params) {
1.41 ! paf 167: double top=params.as_double(0, "range must be expression", r);
1.34 paf 168: if(top<=0 || top>MAX_UINT)
1.16 paf 169: throw Exception("parser.runtime",
1.34 paf 170: 0,
171: "top(%g) must be [1..%u]", top, MAX_UINT);
1.1 parser 172:
1.34 paf 173: r.write_no_lang(*new VInt(_random(uint(top))));
1.1 parser 174: }
175:
176:
1.20 paf 177: typedef double(*math1_func_ptr)(double);
1.3 parser 178: static double frac(double param) { return param-trunc(param); }
179: static double degrees(double param) { return param /PI *180; }
180: static double radians(double param) { return param /180 *PI; }
1.1 parser 181:
1.34 paf 182: static void math1(Request& r, MethodParams& params, math1_func_ptr func) {
1.41 ! paf 183: double param=params.as_double(0, "parameter must be expression", r);
! 184: double result=func(param);
1.34 paf 185: r.write_no_lang(*new VDouble(result));
1.1 parser 186: }
187:
188: #define MATH1(name) \
1.34 paf 189: static void _##name(Request& r, MethodParams& params) {\
190: math1(r, params, &name);\
1.1 parser 191: }
192: #define MATH1P(name_parser, name_c) \
1.34 paf 193: static void _##name_parser(Request& r, MethodParams& params) {\
194: math1(r, params, &name_c);\
1.1 parser 195: }
196: MATH1(round); MATH1(floor); MATH1P(ceiling, ceil);
1.3 parser 197: MATH1(trunc); MATH1(frac);
1.1 parser 198: MATH1P(abs, fabs); MATH1(sign);
199: MATH1(exp); MATH1(log);
200: MATH1(sin); MATH1(asin);
201: MATH1(cos); MATH1(acos);
202: MATH1(tan); MATH1(atan);
1.3 parser 203: MATH1(degrees); MATH1(radians);
1.1 parser 204: MATH1(sqrt);
205:
206:
207: typedef double (*math2_func_ptr)(double, double);
1.34 paf 208: static void math2(Request& r, MethodParams& params, math2_func_ptr func) {
1.41 ! paf 209: double a=params.as_double(0, "parameter must be expression", r);
! 210: double b=params.as_double(1, "parameter must be expression", r);
! 211: double result=func(a, b);
1.34 paf 212: r.write_no_lang(*new VDouble(result));
1.1 parser 213: }
214:
215: #define MATH2(name) \
1.34 paf 216: static void _##name(Request& r, MethodParams& params) {\
217: math2(r, params, &name);\
1.1 parser 218: }
219: MATH2(pow);
220:
1.20 paf 221: inline bool is_salt_body_char(int c) {
222: return isalnum(c) || c == '.' || c=='/';
223: }
1.34 paf 224: static size_t calc_prefix_size(const char* salt) {
1.37 paf 225: if(strlen(salt)) {
1.39 paf 226: if(!is_salt_body_char((unsigned char)salt[0])) { // $... {...
1.34 paf 227: const char* cur=salt+1; // skip
1.39 paf 228: while(is_salt_body_char((unsigned char)*cur++)) // ...$ ...}
1.20 paf 229: ;
230: return cur-salt;
231: } else
232: return 0;
233: } else
234: return 0;
235: }
1.34 paf 236: static void _crypt(Request& r, MethodParams& params) {
237: const char* password=params.as_string(0, "password must be string").cstr();
238: const char* maybe_bodyless_salt=params.as_string(1, "salt must be string").cstr();
1.20 paf 239:
240: size_t prefix_size=calc_prefix_size(maybe_bodyless_salt);
1.34 paf 241: const char* normal_salt;
1.20 paf 242: char normalize_buf[MAX_STRING];
243: if(prefix_size==strlen(maybe_bodyless_salt)) { // bodyless?
244: strncpy(normalize_buf, maybe_bodyless_salt, MAX_STRING-MAX_SALT-1);
245: char *cur=normalize_buf+strlen(normalize_buf);
246: // sould add up MAX_SALT random chars
247: static unsigned char itoa64[] = /* 0 ... 63 => ASCII - 64 */
248: "./0123456789ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz";
249: for(int i=0; i<MAX_SALT; i++)
250: *cur++=itoa64[_random(64)];
251: *cur=0;
252: normal_salt=normalize_buf;
253: } else
254: normal_salt=maybe_bodyless_salt;
1.19 paf 255:
256: /* FreeBSD style MD5 string
257: */
1.20 paf 258: if(strncmp(normal_salt, PA_MD5PW_ID, PA_MD5PW_IDLEN) == 0) {
1.19 paf 259: const size_t sample_size=120;
1.34 paf 260: char *sample_buf=new(PointerFreeGC) char[sample_size];
261: pa_MD5Encode((const unsigned char *)password,
262: (const unsigned char *)normal_salt, sample_buf, sample_size);
263: String sample(sample_buf);
264: r.write_pass_lang(sample);
1.20 paf 265: } else {
266: #ifdef HAVE_CRYPT
1.34 paf 267: const char* static_sample_buf=crypt(password, normal_salt);
268: if(!static_sample_buf // nothing generated
269: || !static_sample_buf[0] // generated nothing
270: || strncmp(static_sample_buf, normal_salt, prefix_size)!=0) // salt prefix not preserved
1.20 paf 271: throw Exception("parser.runtime",
1.34 paf 272: 0,
273: "crypt on this platform does not support '%.*s' salt prefix", prefix_size, normal_salt);
1.20 paf 274:
1.34 paf 275: r.write_pass_lang(String(pa_strdup(static_sample_buf)));
1.20 paf 276: #else
1.19 paf 277: throw Exception("parser.runtime",
1.34 paf 278: 0,
1.19 paf 279: "salt must start with '" PA_MD5PW_ID "'");
1.20 paf 280: #endif
281: }
1.19 paf 282: }
283:
1.34 paf 284: static const char* hex_string(unsigned char* bytes, size_t size, bool upcase) {
285: char *bytes_hex=new(PointerFreeGC) char [size*2/*byte->hh*/+1/*for zero-teminator*/];
1.30 paf 286: unsigned char *src=bytes;
287: unsigned char *end=bytes+size;
288: char *dest=bytes_hex;
1.31 paf 289:
290: static const char *hex=upcase?"0123456789ABCDEF":"0123456789abcdef";
291:
292: for(; src<end; src++) {
293: *dest++=hex[*src/0x10];
294: *dest++=hex[*src%0x10];
295: }
296: *dest=0;
1.30 paf 297:
298: return bytes_hex;
299: }
300:
1.34 paf 301: static void _md5(Request& r, MethodParams& params) {
302: const char *string=params.as_string(0, "parameter must be string").cstr();
1.27 paf 303:
304:
305: PA_MD5_CTX context;
306: unsigned char digest[16];
1.34 paf 307: pa_MD5Init(&context);
308: pa_MD5Update(&context, (const unsigned char*)string, strlen(string));
309: pa_MD5Final(digest, &context);
1.27 paf 310:
1.34 paf 311: r.write_pass_lang(*new String(hex_string(digest, sizeof(digest), false)));
1.26 paf 312: }
313:
1.28 paf 314: /// to hell with extra bytes on 64bit platforms
315: struct uuid {
316: unsigned int time_low;
317: unsigned short time_mid;
318: unsigned short time_hi_and_version;
319: unsigned short clock_seq;
320: unsigned char node[6];
321: };
1.34 paf 322: static void _uuid(Request& r, MethodParams& /*params*/) {
1.28 paf 323:
324: // random
325: struct uuid uuid;
326: random(&uuid, sizeof(uuid));
327:
328: // http://www.opengroup.org/onlinepubs/9629399/apdxa.htm#tagtcjh_35
329: // ~
330: // version = DCE Security version, with embedded POSIX UIDs.
331: // variant = DCE
332: //
333: // DCE=Distributed Computing Environment
334: // http://www.opengroup.org/dce/
335: //
336: // they say this influences comparison&such,
337: // but could not figure out how, hence structure layout specified strictly
338: // anyhow, uuidgen on Win32 yield those values
339: //
340: // xxxxxxxx-xxxx-4xxx-{8,9,A,B}xxx-xxxxxxxxxxxx
341: uuid.clock_seq = (uuid.clock_seq & 0x3FFF) | 0x8000;
342: uuid.time_hi_and_version = (uuid.time_hi_and_version & 0x0FFF) | 0x4000;
343:
344: // format
1.34 paf 345: const int uuid_cstr_bufsize=36+1/*for zero-teminator*/+1/*for faulty snprintfs*/;
346: char *uuid_cstr=new(PointerFreeGC) char[uuid_cstr_bufsize];
1.29 paf 347: snprintf(uuid_cstr, uuid_cstr_bufsize,
1.28 paf 348: "%08X-%04X-%04X-%02X%02X-%02X%02X%02X%02X%02X%02X",
349: uuid.time_low, uuid.time_mid, uuid.time_hi_and_version,
350: uuid.clock_seq >> 8, uuid.clock_seq & 0xFF,
351: uuid.node[0], uuid.node[1], uuid.node[2],
352: uuid.node[3], uuid.node[4], uuid.node[5]);
353:
1.34 paf 354: r.write_pass_lang(*new String(uuid_cstr));
1.28 paf 355: }
356:
1.34 paf 357: static void _uid64(Request& r, MethodParams& /*params*/) {
1.30 paf 358:
359: unsigned char id[64/8];
360: random(&id, sizeof(id));
361:
1.34 paf 362: r.write_pass_lang(*new String(hex_string(id, sizeof(id), true)));
1.30 paf 363: }
364:
1.1 parser 365: // constructor
366:
1.34 paf 367: MMath::MMath(): Methoded("math") {
1.1 parser 368: // ^FUNC(expr)
1.28 paf 369: #define ADDX(name, X) \
370: add_native_method(#name, Method::CT_STATIC, _##name, X, X)
371: #define ADD0(name) ADDX(name, 0)
372: #define ADD1(name) ADDX(name, 1)
373: #define ADD2(name) ADDX(name, 2)
1.1 parser 374:
375: ADD1(round); ADD1(floor); ADD1(ceiling);
1.3 parser 376: ADD1(trunc); ADD1(frac);
1.1 parser 377: ADD1(abs); ADD1(sign);
378: ADD1(exp); ADD1(log);
379: ADD1(sin); ADD1(asin);
380: ADD1(cos); ADD1(acos);
381: ADD1(tan); ADD1(atan);
1.3 parser 382: ADD1(degrees); ADD1(radians);
1.1 parser 383: ADD1(sqrt);
1.3 parser 384: ADD1(random);
1.1 parser 385:
386: // ^pow(x;y)
387: ADD2(pow);
388:
1.19 paf 389: // ^crypt[password;salt]
390: ADD2(crypt);
1.26 paf 391:
392: // ^md5[string]
393: ADD1(md5);
1.5 parser 394:
1.29 paf 395: // ^uuid[]
396: ADD0(uuid);
397:
1.30 paf 398: // ^uid64[]
399: ADD0(uid64);
1.1 parser 400: }
E-mail:
![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to math.C CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [parser3project] / parser3 / src / classes