1 ! Copyright (C) 2009 Slava Pestov, Joe Groff.
2 ! See http://factorcode.org/license.txt for BSD license.
3 USING: accessors alien alien.c-types byte-arrays fry
4 classes.algebra cpu.architecture kernel math sequences
5 math.vectors math.vectors.simd.intrinsics
6 macros generalizations combinators combinators.short-circuit
7 arrays locals compiler.tree.propagation.info
8 compiler.cfg.builder.blocks
9 compiler.cfg.comparisons
10 compiler.cfg.stacks compiler.cfg.stacks.local compiler.cfg.hats
11 compiler.cfg.instructions compiler.cfg.registers
12 compiler.cfg.intrinsics
13 compiler.cfg.intrinsics.alien
14 compiler.cfg.intrinsics.simd.backend
16 FROM: alien.c-types => heap-size char short int longlong float double ;
17 SPECIALIZED-ARRAYS: char uchar short ushort int uint longlong ulonglong float double ;
18 IN: compiler.cfg.intrinsics.simd
22 : sign-bit-mask ( rep -- byte-array )
24 { char-16-rep [ uchar-array{
25 HEX: 80 HEX: 80 HEX: 80 HEX: 80
26 HEX: 80 HEX: 80 HEX: 80 HEX: 80
27 HEX: 80 HEX: 80 HEX: 80 HEX: 80
28 HEX: 80 HEX: 80 HEX: 80 HEX: 80
30 { short-8-rep [ ushort-array{
31 HEX: 8000 HEX: 8000 HEX: 8000 HEX: 8000
32 HEX: 8000 HEX: 8000 HEX: 8000 HEX: 8000
34 { int-4-rep [ uint-array{
35 HEX: 8000,0000 HEX: 8000,0000
36 HEX: 8000,0000 HEX: 8000,0000
38 { longlong-2-rep [ ulonglong-array{
39 HEX: 8000,0000,0000,0000
40 HEX: 8000,0000,0000,0000
44 : ^load-neg-zero-vector ( rep -- dst )
46 { float-4-rep [ float-array{ -0.0 -0.0 -0.0 -0.0 } underlying>> ^^load-constant ] }
47 { double-2-rep [ double-array{ -0.0 -0.0 } underlying>> ^^load-constant ] }
50 : ^load-add-sub-vector ( rep -- dst )
52 { float-4-rep [ float-array{ -0.0 0.0 -0.0 0.0 } underlying>> ^^load-constant ] }
53 { double-2-rep [ double-array{ -0.0 0.0 } underlying>> ^^load-constant ] }
54 { char-16-rep [ char-array{ -1 0 -1 0 -1 0 -1 0 -1 0 -1 0 -1 0 -1 0 } underlying>> ^^load-constant ] }
55 { short-8-rep [ short-array{ -1 0 -1 0 -1 0 -1 0 } underlying>> ^^load-constant ] }
56 { int-4-rep [ int-array{ -1 0 -1 0 } underlying>> ^^load-constant ] }
57 { longlong-2-rep [ longlong-array{ -1 0 } underlying>> ^^load-constant ] }
60 : >variable-shuffle ( shuffle rep -- shuffle' )
61 rep-component-type heap-size
62 [ dup <repetition> >byte-array ]
63 [ iota >byte-array ] bi
64 '[ _ n*v _ v+ ] map concat ;
66 : ^load-immediate-shuffle ( shuffle rep -- dst )
67 >variable-shuffle ^^load-constant ;
69 :: ^blend-vector ( mask true false rep -- dst )
70 true mask rep ^^and-vector
71 mask false rep ^^andn-vector
74 : ^not-vector ( src rep -- dst )
77 [ [ ^^fill-vector ] [ ^^xor-vector ] bi ]
80 :: ^((compare-vector)) ( src1 src2 rep {cc,swap} -- dst )
81 {cc,swap} first2 :> ( cc swap? )
83 [ src2 src1 rep cc ^^compare-vector ]
84 [ src1 src2 rep cc ^^compare-vector ] if ;
86 :: ^(compare-vector) ( src1 src2 rep orig-cc -- dst )
87 rep orig-cc %compare-vector-ccs :> ( ccs not? )
90 [ rep not? [ ^^fill-vector ] [ ^^zero-vector ] if ]
92 ccs unclip :> ( rest-ccs first-cc )
93 src1 src2 rep first-cc ^((compare-vector)) :> first-dst
96 [ [ src1 src2 rep ] dip ^((compare-vector)) rep ^^or-vector ]
99 not? [ rep ^not-vector ] when
102 :: ^minmax-compare-vector ( src1 src2 rep cc -- dst )
104 { cc< [ src1 src2 rep ^^max-vector src1 rep cc/= ^(compare-vector) ] }
105 { cc<= [ src1 src2 rep ^^min-vector src1 rep cc= ^(compare-vector) ] }
106 { cc> [ src1 src2 rep ^^min-vector src1 rep cc/= ^(compare-vector) ] }
107 { cc>= [ src1 src2 rep ^^max-vector src1 rep cc= ^(compare-vector) ] }
110 : ^compare-vector ( src1 src2 rep cc -- dst )
112 [ ^(compare-vector) ]
113 [ ^minmax-compare-vector ]
114 { unsigned-int-vector-rep [| src1 src2 rep cc |
115 rep sign-bit-mask ^^load-constant :> sign-bits
116 src1 sign-bits rep ^^xor-vector
117 src2 sign-bits rep ^^xor-vector
118 rep unsign-rep cc ^(compare-vector)
122 : ^unpack-vector-head ( src rep -- dst )
124 [ ^^unpack-vector-head ]
125 { unsigned-int-vector-rep [ [ ^^zero-vector ] [ ^^merge-vector-head ] bi ] }
126 { signed-int-vector-rep [| src rep |
127 src src rep ^^merge-vector-head :> merged
128 rep rep-component-type heap-size 8 * :> bits
129 merged bits rep widen-vector-rep ^^shr-vector-imm
131 { signed-int-vector-rep [| src rep |
132 rep ^^zero-vector :> zero
133 zero src rep cc> ^compare-vector :> sign
134 src sign rep ^^merge-vector-head
138 : ^unpack-vector-tail ( src rep -- dst )
140 [ ^^unpack-vector-tail ]
141 [ [ ^^tail>head-vector ] [ ^^unpack-vector-head ] bi ]
142 { unsigned-int-vector-rep [ [ ^^zero-vector ] [ ^^merge-vector-tail ] bi ] }
143 { signed-int-vector-rep [| src rep |
144 src src rep ^^merge-vector-tail :> merged
145 rep rep-component-type heap-size 8 * :> bits
146 merged bits rep widen-vector-rep ^^shr-vector-imm
148 { signed-int-vector-rep [| src rep |
149 rep ^^zero-vector :> zero
150 zero src rep cc> ^compare-vector :> sign
151 src sign rep ^^merge-vector-tail
155 : ^(sum-vector-2) ( src rep -- dst )
157 [ dupd ^^horizontal-add-vector ]
159 src src rep ^^merge-vector-head :> head
160 src src rep ^^merge-vector-tail :> tail
161 head tail rep ^^add-vector
165 : ^(sum-vector-4) ( src rep -- dst )
168 [ dupd ^^horizontal-add-vector ]
169 [ dupd ^^horizontal-add-vector ] bi
172 src src rep ^^merge-vector-head :> head
173 src src rep ^^merge-vector-tail :> tail
174 head tail rep ^^add-vector :> src'
176 rep widen-vector-rep :> rep'
177 src' src' rep' ^^merge-vector-head :> head'
178 src' src' rep' ^^merge-vector-tail :> tail'
179 head' tail' rep ^^add-vector
183 : ^(sum-vector-8) ( src rep -- dst )
186 [ dupd ^^horizontal-add-vector ]
187 [ dupd ^^horizontal-add-vector ]
188 [ dupd ^^horizontal-add-vector ] tri
191 src src rep ^^merge-vector-head :> head
192 src src rep ^^merge-vector-tail :> tail
193 head tail rep ^^add-vector :> src'
195 rep widen-vector-rep :> rep'
196 src' src' rep' ^^merge-vector-head :> head'
197 src' src' rep' ^^merge-vector-tail :> tail'
198 head' tail' rep ^^add-vector :> src''
200 rep' widen-vector-rep :> rep''
201 src'' src'' rep'' ^^merge-vector-head :> head''
202 src'' src'' rep'' ^^merge-vector-tail :> tail''
203 head'' tail'' rep ^^add-vector
207 : ^(sum-vector-16) ( src rep -- dst )
211 [ dupd ^^horizontal-add-vector ]
212 [ dupd ^^horizontal-add-vector ]
213 [ dupd ^^horizontal-add-vector ]
214 [ dupd ^^horizontal-add-vector ]
218 src src rep ^^merge-vector-head :> head
219 src src rep ^^merge-vector-tail :> tail
220 head tail rep ^^add-vector :> src'
222 rep widen-vector-rep :> rep'
223 src' src' rep' ^^merge-vector-head :> head'
224 src' src' rep' ^^merge-vector-tail :> tail'
225 head' tail' rep ^^add-vector :> src''
227 rep' widen-vector-rep :> rep''
228 src'' src'' rep'' ^^merge-vector-head :> head''
229 src'' src'' rep'' ^^merge-vector-tail :> tail''
230 head'' tail'' rep ^^add-vector :> src'''
232 rep'' widen-vector-rep :> rep'''
233 src''' src''' rep''' ^^merge-vector-head :> head'''
234 src''' src''' rep''' ^^merge-vector-tail :> tail'''
235 head''' tail''' rep ^^add-vector
239 : ^(sum-vector) ( src rep -- dst )
242 { 2 [ ^(sum-vector-2) ] }
243 { 4 [ ^(sum-vector-4) ] }
244 { 8 [ ^(sum-vector-8) ] }
245 { 16 [ ^(sum-vector-16) ] }
247 ] [ ^^vector>scalar ] bi ;
249 : ^sum-vector ( src rep -- dst )
251 { float-vector-rep [ ^(sum-vector) ] }
252 { int-vector-rep [| src rep |
253 src rep ^unpack-vector-head :> head
254 src rep ^unpack-vector-tail :> tail
255 rep widen-vector-rep :> wide-rep
256 head tail wide-rep ^^add-vector wide-rep
261 : shuffle? ( obj -- ? ) { [ array? ] [ [ integer? ] all? ] } 1&& ;
263 : ^shuffle-vector-imm ( src1 shuffle rep -- dst )
264 [ rep-length 0 pad-tail ] keep {
265 [ ^^shuffle-vector-imm ]
266 [ [ ^load-immediate-shuffle ] [ ^^shuffle-vector ] bi ]
269 : ^broadcast-vector ( src n rep -- dst )
270 [ rep-length swap <array> ] keep
271 ^shuffle-vector-imm ;
273 : ^with-vector ( src rep -- dst )
274 [ ^^scalar>vector ] keep [ 0 ] dip ^broadcast-vector ;
276 : ^select-vector ( src n rep -- dst )
277 [ ^broadcast-vector ] keep ^^vector>scalar ;
281 : emit-simd-v+ ( node -- )
284 } emit-vv-vector-op ;
286 : emit-simd-v- ( node -- )
289 } emit-vv-vector-op ;
291 : emit-simd-vneg ( node -- )
293 { float-vector-rep [ [ ^load-neg-zero-vector ] [ ^^sub-vector ] bi ] }
294 { int-vector-rep [ [ ^^zero-vector ] [ ^^sub-vector ] bi ] }
297 : emit-simd-v+- ( node -- )
300 { float-vector-rep [| src1 src2 rep |
301 rep ^load-add-sub-vector :> signs
302 src2 signs rep ^^xor-vector :> src2'
303 src1 src2' rep ^^add-vector
305 { int-vector-rep [| src1 src2 rep |
306 rep ^load-add-sub-vector :> signs
307 src2 signs rep ^^xor-vector :> src2'
308 src2' signs rep ^^sub-vector :> src2''
309 src1 src2'' rep ^^add-vector
311 } emit-vv-vector-op ;
313 : emit-simd-vs+ ( node -- )
315 { float-vector-rep [ ^^add-vector ] }
316 { int-vector-rep [ ^^saturated-add-vector ] }
317 } emit-vv-vector-op ;
319 : emit-simd-vs- ( node -- )
321 { float-vector-rep [ ^^sub-vector ] }
322 { int-vector-rep [ ^^saturated-sub-vector ] }
323 } emit-vv-vector-op ;
325 : emit-simd-vs* ( node -- )
327 { float-vector-rep [ ^^mul-vector ] }
328 { int-vector-rep [ ^^saturated-mul-vector ] }
329 } emit-vv-vector-op ;
331 : emit-simd-v* ( node -- )
334 } emit-vv-vector-op ;
336 : emit-simd-v/ ( node -- )
339 } emit-vv-vector-op ;
341 : emit-simd-vmin ( node -- )
345 [ cc< ^compare-vector ]
346 [ ^blend-vector ] 3bi
348 } emit-vv-vector-op ;
350 : emit-simd-vmax ( node -- )
354 [ cc> ^compare-vector ]
355 [ ^blend-vector ] 3bi
357 } emit-vv-vector-op ;
359 : emit-simd-v. ( node -- )
362 { float-vector-rep [ [ ^^mul-vector ] [ ^sum-vector ] bi ] }
363 } emit-vv-vector-op ;
365 : emit-simd-vsqrt ( node -- )
370 : emit-simd-sum ( node -- )
375 : emit-simd-vabs ( node -- )
377 { unsigned-int-vector-rep [ drop ] }
379 { float-vector-rep [ [ ^load-neg-zero-vector ] [ swapd ^^andn-vector ] bi ] }
380 { int-vector-rep [| src rep |
381 rep ^^zero-vector :> zero
382 zero src rep ^^sub-vector :> -src
383 zero src rep cc> ^compare-vector :> sign
384 sign -src src rep ^blend-vector
388 : emit-simd-vand ( node -- )
391 } emit-vv-vector-op ;
393 : emit-simd-vandn ( node -- )
396 } emit-vv-vector-op ;
398 : emit-simd-vor ( node -- )
401 } emit-vv-vector-op ;
403 : emit-simd-vxor ( node -- )
406 } emit-vv-vector-op ;
408 : emit-simd-vnot ( node -- )
413 : emit-simd-vlshift ( node -- )
418 } [ integer? ] emit-vv-or-vl-vector-op ;
420 : emit-simd-vrshift ( node -- )
425 } [ integer? ] emit-vv-or-vl-vector-op ;
427 : emit-simd-hlshift ( node -- )
429 [ ^^horizontal-shl-vector-imm ]
430 } [ integer? ] emit-vl-vector-op ;
432 : emit-simd-hrshift ( node -- )
434 [ ^^horizontal-shr-vector-imm ]
435 } [ integer? ] emit-vl-vector-op ;
437 : emit-simd-vshuffle-elements ( node -- )
439 [ ^shuffle-vector-imm ]
440 } [ shuffle? ] emit-vl-vector-op ;
442 : emit-simd-vshuffle-bytes ( node -- )
445 } emit-vv-vector-op ;
447 : emit-simd-vmerge-head ( node -- )
449 [ ^^merge-vector-head ]
450 } emit-vv-vector-op ;
452 : emit-simd-vmerge-tail ( node -- )
454 [ ^^merge-vector-tail ]
455 } emit-vv-vector-op ;
457 : emit-simd-v<= ( node -- )
459 [ cc<= ^compare-vector ]
460 } emit-vv-vector-op ;
461 : emit-simd-v< ( node -- )
463 [ cc< ^compare-vector ]
464 } emit-vv-vector-op ;
465 : emit-simd-v= ( node -- )
467 [ cc= ^compare-vector ]
468 } emit-vv-vector-op ;
469 : emit-simd-v> ( node -- )
471 [ cc> ^compare-vector ]
472 } emit-vv-vector-op ;
473 : emit-simd-v>= ( node -- )
475 [ cc>= ^compare-vector ]
476 } emit-vv-vector-op ;
477 : emit-simd-vunordered? ( node -- )
479 [ cc/<>= ^compare-vector ]
480 } emit-vv-vector-op ;
482 : emit-simd-vany? ( node -- )
484 [ vcc-any ^^test-vector ]
486 : emit-simd-vall? ( node -- )
488 [ vcc-all ^^test-vector ]
490 : emit-simd-vnone? ( node -- )
492 [ vcc-none ^^test-vector ]
495 : emit-simd-v>float ( node -- )
497 { float-vector-rep [ drop ] }
498 { int-vector-rep [ ^^integer>float-vector ] }
501 : emit-simd-v>integer ( node -- )
503 { float-vector-rep [ ^^float>integer-vector ] }
504 { int-vector-rep [ drop ] }
507 : emit-simd-vpack-signed ( node -- )
509 [ ^^signed-pack-vector ]
510 } emit-vv-vector-op ;
512 : emit-simd-vpack-unsigned ( node -- )
514 [ ^^unsigned-pack-vector ]
515 } emit-vv-vector-op ;
517 : emit-simd-vunpack-head ( node -- )
519 [ ^unpack-vector-head ]
522 : emit-simd-vunpack-tail ( node -- )
524 [ ^unpack-vector-tail ]
527 : emit-simd-with ( node -- )
532 : emit-simd-gather-2 ( node -- )
534 [ ^^gather-vector-2 ]
535 } emit-vv-vector-op ;
537 : emit-simd-gather-4 ( node -- )
539 [ ^^gather-vector-4 ]
540 } emit-vvvv-vector-op ;
542 : emit-simd-select ( node -- )
545 } [ integer? ] emit-vl-vector-op ;
547 : emit-alien-vector ( node -- )
550 ds-drop prepare-alien-getter
551 _ ^^alien-vector ds-push
553 [ inline-alien-getter? ] inline-alien
554 ] with { [ %alien-vector-reps member? ] } if-literals-match ;
556 : emit-set-alien-vector ( node -- )
559 ds-drop prepare-alien-setter ds-pop
562 [ byte-array inline-alien-setter? ]
564 ] with { [ %alien-vector-reps member? ] } if-literals-match ;
568 { (simd-v+) [ emit-simd-v+ ] }
569 { (simd-v-) [ emit-simd-v- ] }
570 { (simd-vneg) [ emit-simd-vneg ] }
571 { (simd-v+-) [ emit-simd-v+- ] }
572 { (simd-vs+) [ emit-simd-vs+ ] }
573 { (simd-vs-) [ emit-simd-vs- ] }
574 { (simd-vs*) [ emit-simd-vs* ] }
575 { (simd-v*) [ emit-simd-v* ] }
576 { (simd-v/) [ emit-simd-v/ ] }
577 { (simd-vmin) [ emit-simd-vmin ] }
578 { (simd-vmax) [ emit-simd-vmax ] }
579 { (simd-v.) [ emit-simd-v. ] }
580 { (simd-vsqrt) [ emit-simd-vsqrt ] }
581 { (simd-sum) [ emit-simd-sum ] }
582 { (simd-vabs) [ emit-simd-vabs ] }
583 { (simd-vbitand) [ emit-simd-vand ] }
584 { (simd-vbitandn) [ emit-simd-vandn ] }
585 { (simd-vbitor) [ emit-simd-vor ] }
586 { (simd-vbitxor) [ emit-simd-vxor ] }
587 { (simd-vbitnot) [ emit-simd-vnot ] }
588 { (simd-vand) [ emit-simd-vand ] }
589 { (simd-vandn) [ emit-simd-vandn ] }
590 { (simd-vor) [ emit-simd-vor ] }
591 { (simd-vxor) [ emit-simd-vxor ] }
592 { (simd-vnot) [ emit-simd-vnot ] }
593 { (simd-vlshift) [ emit-simd-vlshift ] }
594 { (simd-vrshift) [ emit-simd-vrshift ] }
595 { (simd-hlshift) [ emit-simd-hlshift ] }
596 { (simd-hrshift) [ emit-simd-hrshift ] }
597 { (simd-vshuffle-elements) [ emit-simd-vshuffle-elements ] }
598 { (simd-vshuffle-bytes) [ emit-simd-vshuffle-bytes ] }
599 { (simd-vmerge-head) [ emit-simd-vmerge-head ] }
600 { (simd-vmerge-tail) [ emit-simd-vmerge-tail ] }
601 { (simd-v<=) [ emit-simd-v<= ] }
602 { (simd-v<) [ emit-simd-v< ] }
603 { (simd-v=) [ emit-simd-v= ] }
604 { (simd-v>) [ emit-simd-v> ] }
605 { (simd-v>=) [ emit-simd-v>= ] }
606 { (simd-vunordered?) [ emit-simd-vunordered? ] }
607 { (simd-vany?) [ emit-simd-vany? ] }
608 { (simd-vall?) [ emit-simd-vall? ] }
609 { (simd-vnone?) [ emit-simd-vnone? ] }
610 { (simd-v>float) [ emit-simd-v>float ] }
611 { (simd-v>integer) [ emit-simd-v>integer ] }
612 { (simd-vpack-signed) [ emit-simd-vpack-signed ] }
613 { (simd-vpack-unsigned) [ emit-simd-vpack-unsigned ] }
614 { (simd-vunpack-head) [ emit-simd-vunpack-head ] }
615 { (simd-vunpack-tail) [ emit-simd-vunpack-tail ] }
616 { (simd-with) [ emit-simd-with ] }
617 { (simd-gather-2) [ emit-simd-gather-2 ] }
618 { (simd-gather-4) [ emit-simd-gather-4 ] }
619 { (simd-select) [ emit-simd-select ] }
620 { alien-vector [ emit-alien-vector ] }
621 { set-alien-vector [ emit-set-alien-vector ] }
622 } enable-intrinsics ;