"AIscm" array JIT

[Jan Wedekind <jan@wedesoft.de>]

Hi,
   I am working on a compact library [1] for JIT compilation of array 
operations. It only runs on AMD64 processors. Currently it supports array 
operations using booleans, integers, and integer RGB and integer complex 
numbers.
   There are still important things missing: floating point numbers, 
compiling calls to C methods (e.g. sin, cos, ...), tensor operations, 
convolutions, ... I would like to eventually do numerical processing 
similar to Python's NumPy (but more generic), Theano (but more compact 
syntax as facilitated by macros), and OpenCV.
   Here is an example adding an integer to each element of a 2D array and 
returning the result:

     scheme@(guile-user)> (use-modules (oop goops) (aiscm jit) (aiscm int)
                                       (aiscm pointer) (aiscm sequence))
     scheme@(guile-user)> (+ (arr (2 3 5) (7 11 13)) 3)
     $1 = #<sequence<sequence<int<8,unsigned>>>>:
     ((5 6 8)
      (10 14 16))

   The fallback method for the GOOPS generic "+" adds a JIT compiled plus 
operation for the specific array types to the generic and then calls "+" 
again.
   The corresponding machine code to produce the unsigned byte array is 
shown below:

      0:   4c 89 64 24 f0          mov    QWORD PTR [rsp-0x10],r12
      5:   48 89 6c 24 e8          mov    QWORD PTR [rsp-0x18],rbp
      a:   4c 89 7c 24 e0          mov    QWORD PTR [rsp-0x20],r15
      f:   4c 89 74 24 d8          mov    QWORD PTR [rsp-0x28],r14
     14:   4c 89 6c 24 d0          mov    QWORD PTR [rsp-0x30],r13
     19:   48 89 5c 24 c8          mov    QWORD PTR [rsp-0x38],rbx
     1e:   48 89 7c 24 f8          mov    QWORD PTR [rsp-0x8],rdi
     23:   4c 8b 64 24 08          mov    r12,QWORD PTR [rsp+0x8]
     28:   48 8b 7c 24 18          mov    rdi,QWORD PTR [rsp+0x18]
     2d:   48 8b 6c 24 20          mov    rbp,QWORD PTR [rsp+0x20]
     32:   8a 44 24 28             mov    al,BYTE PTR [rsp+0x28]
     36:   48 6b de 01             imul   rbx,rsi,0x1
     3a:   49 8b f0                mov    rsi,r8
     3d:   4d 6b cc 01             imul   r9,r12,0x1
     41:   4c 8b fd                mov    r15,rbp
     44:   49 be 00 00 00 00 00    movabs r14,0x0
     4b:   00 00 00
     4e:   4c 8b 44 24 f8          mov    r8,QWORD PTR [rsp-0x8]
     53:   4d 3b f0                cmp    r14,r8
     56:   74 3e                   je     0x96
     58:   49 ff c6                inc    r14
     5b:   4c 6b d9 01             imul   r11,rcx,0x1
     5f:   4c 8b ee                mov    r13,rsi
     62:   4c 6b d7 01             imul   r10,rdi,0x1
     66:   4d 8b e7                mov    r12,r15
     69:   48 bd 00 00 00 00 00    movabs rbp,0x0
     70:   00 00 00
     73:   48 3b ea                cmp    rbp,rdx
     76:   74 16                   je     0x8e
     78:   48 ff c5                inc    rbp
     7b:   45 8a 04 24             mov    r8b,BYTE PTR [r12]
     7f:   44 02 c0                add    r8b,al
     82:   45 88 45 00             mov    BYTE PTR [r13+0x0],r8b
     86:   4d 03 eb                add    r13,r11
     89:   4d 03 e2                add    r12,r10
     8c:   eb e5                   jmp    0x73
     8e:   48 03 f3                add    rsi,rbx
     91:   4d 03 f9                add    r15,r9
     94:   eb b8                   jmp    0x4e
     96:   4c 8b 64 24 f0          mov    r12,QWORD PTR [rsp-0x10]
     9b:   48 8b 6c 24 e8          mov    rbp,QWORD PTR [rsp-0x18]
     a0:   4c 8b 7c 24 e0          mov    r15,QWORD PTR [rsp-0x20]
     a5:   4c 8b 74 24 d8          mov    r14,QWORD PTR [rsp-0x28]
     aa:   4c 8b 6c 24 d0          mov    r13,QWORD PTR [rsp-0x30]
     af:   48 8b 5c 24 c8          mov    rbx,QWORD PTR [rsp-0x38]
     b4:   c3                      ret



   Any comments,suggestions, and feedback are welcome!

Regards
Jan

[1] https://github.com/wedesoft/aiscm