485 lines
11 KiB
PHP
485 lines
11 KiB
PHP
MAX_SPHERES = 4 ;8
|
|
MAX_LIGHTS = 6
|
|
|
|
main_loop:
|
|
call normalize_surface
|
|
xor eax,eax ; y
|
|
.next_line:
|
|
xor ebx,ebx ; x
|
|
|
|
@@:
|
|
push eax
|
|
push ebx
|
|
|
|
|
|
|
|
call find_intersection
|
|
pop ebx
|
|
pop eax
|
|
add ebx,1
|
|
cmp ebx,XRES
|
|
jnz @b
|
|
|
|
add eax,1
|
|
cmp eax,YRES
|
|
jnz .next_line
|
|
|
|
ret
|
|
|
|
normalize_surface:
|
|
movss xmm0,[surface+8]
|
|
movlhps xmm0,xmm0
|
|
movlps xmm0,[surface]
|
|
movaps xmm1,xmm0
|
|
mulps xmm0,xmm0
|
|
haddps xmm0,xmm0
|
|
haddps xmm0,xmm0
|
|
sqrtss xmm0,xmm0
|
|
shufps xmm0,xmm0,0
|
|
divps xmm1,xmm0
|
|
movaps xmm2,xmm1
|
|
pslldq xmm2,4
|
|
psrldq xmm2,4
|
|
movaps [surface_without_d],xmm2
|
|
movlps [surface],xmm1
|
|
movhlps xmm1,xmm1
|
|
movss [surface+8],xmm1
|
|
ret
|
|
|
|
find_intersection:
|
|
; eax - y
|
|
; ebx - x
|
|
push ebp
|
|
mov ebp,esp
|
|
sub esp,256
|
|
and ebp,0xfffffff0
|
|
.dz equ dword[ebp-8]
|
|
.dy equ dword[ebp-12]
|
|
.dx equ [ebp-16]
|
|
.a equ dword[ebp-20]
|
|
.b equ dword[ebp-24]
|
|
.c equ dword[ebp-28]
|
|
.delta equ dword[ebp-32]
|
|
.iy equ dword[ebp-36]
|
|
.ix equ [ebp-40]
|
|
.t1 equ [ebp-44]
|
|
.t2 equ [ebp-48]
|
|
.n equ [ebp-64]
|
|
.sph_xyz equ [ebp-80]
|
|
.one_write equ byte [ebp-81] ;tells if sth written in 'nearest' data
|
|
.sph_counter equ dword [ebp-85]
|
|
.dx_sh equ [ebp-96]
|
|
.a_sh equ [ebp-100]
|
|
.b_sh equ [ebp-104]
|
|
.c_sh equ [ebp-108]
|
|
.delta_sh equ [ebp-112]
|
|
.sph_counter_sh equ dword[ebp-116]
|
|
.shadow_mark equ dword [ebp-120]
|
|
.nearest_surface equ [ebp-144]
|
|
.vd equ [ebp-148] ; denominator in plane inters. equation
|
|
|
|
mov .iy,eax
|
|
mov .ix,ebx
|
|
|
|
mov .one_write,0
|
|
mov .shadow_mark,0
|
|
xorps xmm0,xmm0
|
|
cvtpi2ps xmm0,.ix
|
|
mov ecx,XRES
|
|
cvtsi2ss xmm2,ecx
|
|
shufps xmm2,xmm2,0
|
|
divps xmm0,xmm2
|
|
subps xmm0,[camera]
|
|
movaps .dx,xmm0
|
|
movaps xmm1,xmm0
|
|
mulps xmm1,xmm0
|
|
haddps xmm1,xmm1
|
|
haddps xmm1,xmm1
|
|
movss .a,xmm1
|
|
mov .sph_counter,0
|
|
.next_sph: ; intersection with sphere
|
|
movaps xmm5,[camera]
|
|
mov edx,.sph_counter
|
|
shl edx,4
|
|
add edx,sphere
|
|
subps xmm5,[edx]
|
|
mulps xmm5,[float2]
|
|
mulps xmm5,.dx
|
|
haddps xmm5,xmm5
|
|
haddps xmm5,xmm5
|
|
movss .b,xmm5
|
|
|
|
|
|
movaps xmm4,[edx]
|
|
mulps xmm4,xmm4
|
|
movaps xmm5,[camera]
|
|
mulps xmm5,xmm5
|
|
addps xmm4,xmm5
|
|
haddps xmm4,xmm4
|
|
haddps xmm4,xmm4
|
|
movaps xmm5,[edx] ;;[sphere] ;; [edx]
|
|
mulps xmm5,[camera]
|
|
haddps xmm5,xmm5
|
|
haddps xmm5,xmm5
|
|
mulss xmm5,[float2]
|
|
subss xmm4,xmm5
|
|
mov ebx,.sph_counter
|
|
shl ebx,2
|
|
add ebx,sph_radius
|
|
movss xmm5,[ebx] ;[R] ; [ebx]
|
|
mulss xmm5,xmm5
|
|
subss xmm4,xmm5
|
|
movss .c,xmm4
|
|
|
|
movss xmm5,.b
|
|
mulss xmm5,xmm5
|
|
mulss xmm4,.a
|
|
mulss xmm4,[float4]
|
|
subss xmm5,xmm4
|
|
movss .delta,xmm5
|
|
xorps xmm6,xmm6
|
|
cmpnltss xmm5,xmm6
|
|
movd ecx,xmm5 ; ecx = -1 => greater than 0.0
|
|
cmp ecx,0
|
|
jnz @f
|
|
jmp .next_s ; no intersection
|
|
|
|
|
|
@@:
|
|
movss xmm5,.delta
|
|
sqrtss xmm5,xmm5
|
|
movss xmm4,xmm5
|
|
subss xmm6,.b
|
|
movss xmm7,xmm6
|
|
subss xmm6,xmm5
|
|
divss xmm6,[float2]
|
|
divss xmm6,.a
|
|
movss .t1,xmm6
|
|
addss xmm4,xmm7
|
|
divss xmm4,[float2]
|
|
divss xmm4,.a
|
|
movss .t2,xmm4
|
|
maxss xmm6,xmm4
|
|
cmp .one_write,0 ; test if sth in 'nearest' data is written
|
|
jz @f
|
|
movss xmm4,xmm6 ;5
|
|
cmpnltss xmm4,[smalest_t]
|
|
movd ecx,xmm4
|
|
or ecx,ecx
|
|
jz .next_s
|
|
@@:
|
|
movss [smalest_t],xmm6 ;5
|
|
movaps xmm0,[edx]
|
|
movaps [nearest_sphere],xmm0
|
|
push dword[ebx]
|
|
pop dword[nearest_radius]
|
|
mov .one_write,1 ; one_write - object index -> 1 = sphere
|
|
.next_s:
|
|
add .sph_counter,1
|
|
cmp .sph_counter,MAX_SPHERES
|
|
jnz .next_sph
|
|
|
|
if 1
|
|
movaps xmm0,[surface_without_d] ; find with plane intersection
|
|
mulps xmm0,[camera] ; only one surface is computed
|
|
haddps xmm0,xmm0
|
|
haddps xmm0,xmm0
|
|
addss xmm0,[surface+12]
|
|
movaps xmm1,[surface_without_d]
|
|
mulps xmm1,.dx
|
|
haddps xmm1,xmm1
|
|
haddps xmm1,xmm1
|
|
xorps xmm2,xmm2
|
|
cmpnless xmm2,xmm1
|
|
movd ecx,xmm2
|
|
cmp ecx,0 ;-1
|
|
je .put_pixel ; denominator equal 'zero' - no intersection
|
|
xorps xmm2,xmm2 ; denominator > 0 -> inters. not in screen area
|
|
movss .vd,xmm1 ; write to memory this denom.
|
|
divss xmm0,xmm1
|
|
subss xmm2,xmm0
|
|
cmp .one_write,0
|
|
jz @f
|
|
movss xmm0,xmm2
|
|
cmpnltss xmm2,[smalest_t]
|
|
movd ecx,xmm2
|
|
cmp ecx,0
|
|
je .put_pixel
|
|
@@:
|
|
movss [smalest_t],xmm0
|
|
; test [smalest_t],0x80000000
|
|
; jz @f
|
|
; and [smalest_t],0x7fffffff
|
|
; @@:
|
|
movaps xmm2,[surface]
|
|
movaps .nearest_surface,xmm2
|
|
mov .one_write,2 ; nearest object -> 2 = flat plane
|
|
|
|
end if
|
|
|
|
.put_pixel:
|
|
cmp .one_write,0 ; end if no intersection
|
|
je .end
|
|
|
|
movss xmm5,[smalest_t]
|
|
shufps xmm5,xmm5,0 ; calc and put pixel
|
|
movaps xmm6,.dx
|
|
mulps xmm6,xmm5
|
|
movaps xmm4,[camera]
|
|
addps xmm4,xmm6 ; xmm4 - x,y,z on the sphere or on surface
|
|
movaps .sph_xyz,xmm4
|
|
if 1
|
|
cmp .one_write,2
|
|
jne .shadow
|
|
; movaps xmm4,.sph_xyz
|
|
movaps xmm7,xmm4
|
|
lea ebx,.nearest_surface
|
|
movss xmm4,[ebx+8]
|
|
movlhps xmm4,xmm4
|
|
movlps xmm4,[ebx] ; xmm4 - normal to surface vector
|
|
test dword .vd,0x80000000
|
|
jz @f
|
|
andps xmm4,[positive_mask] ;0x7fffffff
|
|
@@:
|
|
|
|
jmp .calc_pix
|
|
|
|
end if
|
|
.shadow:
|
|
if 1
|
|
; to find shadow intersect:
|
|
; P0 - point on sphere
|
|
; P1 - light
|
|
; with every other sphere in scene if any intersection occured -
|
|
; - point is in shadow
|
|
; next_sph_shad:
|
|
mov ecx,MAX_LIGHTS
|
|
.next_light_sh:
|
|
push ecx
|
|
shl ecx,4
|
|
movaps xmm0,[ecx+light] ;xmm4 - point on nearest sphere
|
|
subps xmm0,xmm4
|
|
movaps .dx_sh,xmm0
|
|
|
|
mulps xmm0,xmm0
|
|
haddps xmm0,xmm0
|
|
haddps xmm0,xmm0
|
|
movss .a_sh,xmm0
|
|
|
|
mov .sph_counter_sh,0
|
|
.next_sph_sh: ; be sure you not intersect nearest sphere with itself
|
|
movaps xmm5,.sph_xyz ;[light]
|
|
mov edx,.sph_counter_sh
|
|
shl edx,4
|
|
add edx,sphere
|
|
movaps xmm7,[edx]
|
|
cmpeqps xmm7,[nearest_sphere]
|
|
movmskps ecx,xmm7
|
|
and ecx,0111b
|
|
cmp ecx,0
|
|
jne .next_s_sh
|
|
|
|
subps xmm5,[edx] ; [edx] - cur sph
|
|
mulps xmm5,[float2]
|
|
mulps xmm5,.dx_sh
|
|
haddps xmm5,xmm5
|
|
haddps xmm5,xmm5
|
|
movss .b_sh,xmm5
|
|
|
|
movaps xmm4,[edx]
|
|
mulps xmm4,xmm4
|
|
movaps xmm5,.sph_xyz
|
|
mulps xmm5,xmm5
|
|
addps xmm4,xmm5
|
|
haddps xmm4,xmm4
|
|
haddps xmm4,xmm4
|
|
movaps xmm5,.sph_xyz
|
|
mulps xmm5,[edx]
|
|
haddps xmm5,xmm5
|
|
haddps xmm5,xmm5
|
|
mulss xmm5,[float2]
|
|
subss xmm4,xmm5
|
|
mov ebx,.sph_counter_sh
|
|
shl ebx,2
|
|
add ebx,sph_radius
|
|
movss xmm5,[ebx]
|
|
mulss xmm5,xmm5
|
|
subss xmm4,xmm5
|
|
movss .c_sh,xmm4
|
|
|
|
movss xmm5,.b_sh
|
|
mulss xmm5,xmm5
|
|
mulss xmm4,.a_sh
|
|
mulss xmm4,[float4]
|
|
subss xmm5,xmm4
|
|
movss .delta_sh,xmm5
|
|
xorps xmm6,xmm6
|
|
cmpnltss xmm5,xmm6
|
|
movd ecx,xmm5 ; ecx = -1 greater than 0.0
|
|
cmp ecx,0
|
|
jnz @f
|
|
jmp .next_s_sh ; no intersection
|
|
@@:
|
|
add .shadow_mark,1 ; mark ->point in shadow
|
|
pop ecx
|
|
sub ecx,1
|
|
jnz .next_light_sh
|
|
jmp .put_pix
|
|
.next_s_sh:
|
|
add .sph_counter_sh,1
|
|
cmp .sph_counter_sh,MAX_SPHERES
|
|
jnz .next_sph_sh
|
|
pop ecx
|
|
sub ecx,1
|
|
jnz .next_light_sh
|
|
|
|
end if
|
|
|
|
|
|
|
|
|
|
|
|
|
|
.put_pix:
|
|
movaps xmm4,.sph_xyz
|
|
movaps xmm7,xmm4
|
|
subps xmm4,[nearest_sphere]
|
|
|
|
movss xmm0,[nearest_radius]
|
|
shufps xmm0,xmm0,0
|
|
divps xmm4,xmm0 ; xmm4 - normal to surface vector
|
|
.calc_pix: ; normal computed
|
|
movaps xmm1,xmm4 ; copy of normal in xmm1
|
|
xor eax,eax
|
|
xorps xmm3,xmm3
|
|
; movss xmm2,[light_factor] ; other model of lighting
|
|
; shufps xmm2,xmm2,0
|
|
|
|
.next_light:
|
|
mov ebx,eax
|
|
shl ebx,4
|
|
movaps xmm5,[light+ebx]
|
|
subps xmm5,xmm7 ; calc light unit vector
|
|
movaps xmm6,xmm5
|
|
mulps xmm5,xmm5
|
|
haddps xmm5,xmm5
|
|
haddps xmm5,xmm5
|
|
sqrtss xmm5,xmm5
|
|
shufps xmm5,xmm5,0
|
|
divps xmm6,xmm5 ; xmm6 - normalized light vector
|
|
; dot_product
|
|
movaps xmm4,xmm1 ; xmm4 - normal to surface
|
|
mulps xmm4,xmm6
|
|
haddps xmm4,xmm4
|
|
haddps xmm4,xmm4
|
|
shufps xmm4,xmm4,0
|
|
mulps xmm4,[lights_color+ebx] ; xmm4 - computed col. light vector dep.
|
|
|
|
; mulps xmm4,xmm2 ; other model of lighting
|
|
; addps xmm3,xmm4
|
|
|
|
maxps xmm3,xmm4 ; will be this better ?
|
|
|
|
add eax,1
|
|
cmp eax,MAX_LIGHTS
|
|
jnz .next_light
|
|
|
|
minps xmm3,[float255]
|
|
cmp .shadow_mark,0
|
|
je @f
|
|
cvtsi2ss xmm2,.shadow_mark
|
|
shufps xmm2,xmm2,0
|
|
mulps xmm2,[shadow_factor]
|
|
subps xmm3,xmm2
|
|
xorps xmm0,xmm0
|
|
maxps xmm3,xmm0
|
|
@@:
|
|
|
|
cvtps2dq xmm3,xmm3
|
|
packssdw xmm3,xmm3
|
|
packuswb xmm3,xmm3
|
|
paddusb xmm3,[ambient_col]
|
|
|
|
if 1
|
|
cmp .one_write,2
|
|
jne .perspective
|
|
movaps xmm0,xmm3 ; calc texture on plane
|
|
movaps xmm1,.sph_xyz
|
|
mov ecx,XRES
|
|
cvtsi2ss xmm4,ecx
|
|
shufps xmm4,xmm4,0
|
|
mulps xmm1,xmm4
|
|
cvtps2dq xmm1,xmm1
|
|
movd ecx,xmm1
|
|
test ecx,0x8
|
|
jz @f
|
|
mov ecx,0xffffffff
|
|
jmp .next_tex_test
|
|
@@:
|
|
xor ecx,ecx
|
|
.next_tex_test:
|
|
psrldq xmm1,8
|
|
movd ebx,xmm1
|
|
test ebx,0x20
|
|
jz @f
|
|
mov ebx,0xffffffff
|
|
jmp .set_tex
|
|
@@:
|
|
xor ebx,ebx
|
|
.set_tex:
|
|
xor ebx,ecx
|
|
shr ebx,28
|
|
mov bh,bl
|
|
movd xmm7,ebx
|
|
paddusb xmm3,xmm7
|
|
|
|
end if
|
|
|
|
|
|
.perspective:
|
|
movaps xmm0,.sph_xyz ; perspective correction
|
|
subps xmm0,[camera]
|
|
movss xmm1,[camera]
|
|
movss xmm2,xmm0
|
|
movaps xmm4,xmm0
|
|
shufps xmm4,xmm4,00000010b
|
|
divss xmm2,xmm4
|
|
mulss xmm2,[camera+8]
|
|
subss xmm1,xmm2 ; xmm1 - x
|
|
|
|
movaps xmm2,xmm0
|
|
shufps xmm2,xmm2,00000001b
|
|
movaps xmm4,xmm0
|
|
shufps xmm4,xmm4,00000010b
|
|
divss xmm2,xmm4
|
|
mulss xmm2,[camera+8]
|
|
movss xmm4,[camera+4]
|
|
subss xmm4,xmm2 ; xmm4 - y
|
|
|
|
mov ebx,XRES
|
|
cvtsi2ss xmm2,ebx
|
|
mulss xmm1,xmm2
|
|
mulss xmm4,xmm2
|
|
cvtss2si ecx,xmm1
|
|
mov .ix,ecx
|
|
cvtss2si edx,xmm4
|
|
mov .iy,edx
|
|
|
|
mov edi,screen
|
|
mov ecx,XRES
|
|
imul ecx,.iy
|
|
add ecx,.ix
|
|
lea ecx,[ecx*3]
|
|
add edi,ecx
|
|
movd [edi],xmm3
|
|
.end:
|
|
add esp,256
|
|
pop ebp
|
|
|
|
ret
|
|
|
|
|
|
|
|
|
|
|