用内联汇编写的RT核心

__asm {

; 152 : ray_dir = ray->dir;

; 153 :

; 154 : dot_nd = - ( normal.x * ray_dir.x + normal.y * ray_dir.y + normal.z * ray_dir.z );

xorps xmm3, xmm3 ;

mov ebx, dword ptr ray ;

mov edi, dword ptr state ;

mov esi, dword ptr this ;

mov ecx, dword ptr [ebx]e_Ray.dir.x ;

mov edx, dword ptr [ebx]e_Ray.dir.y ;

mov eax, dword ptr [ebx]e_Ray.dir.z ;

movss xmm0, dword ptr [esi]this.normal.z ;

movss xmm1, dword ptr [esi]this.normal.y ;

prefetchnta byte ptr [esi+TYPE this] ;

mov dword ptr ray_dir.x, ecx ;

mov dword ptr ray_dir.y, edx ;

mulss xmm1, dword ptr ray_dir.y ;

mov dword ptr ray_dir.z, eax ;

mulss xmm0, dword ptr ray_dir.z ;

addss xmm0, xmm1 ;

movss xmm1, dword ptr [esi]this.normal.x ;

mulss xmm1, dword ptr ray_dir.x ;

; 156 : ray_src = ray->src;

mov edx, dword ptr [ebx]e_Ray.src.x ;

mov eax, dword ptr [ebx]e_Ray.src.y ;

mov ecx, dword ptr [ebx]e_Ray.src.z ;

addss xmm0, xmm1 ;

movaps xmm1, xmm3 ;

subss xmm1, xmm0 ;

; 158 : if( dot_nd <= 0.0f )

comiss xmm3, xmm1 ;

movss dword ptr dot_nd, xmm1 ;

mov dword ptr ray_src.x, edx ;

mov dword ptr ray_src.y, eax ;

mov dword ptr ray_src.z, ecx ;

; 159 : return false;

jae RETURN_FALSE ;

; 161 : t = ray_src.x * normal.x + ray_src.y * normal.y + ray_src.z * normal.z + normal.w;

movss xmm7, dword ptr [esi]this.normal.z ;

mulss xmm7, dword ptr ray_src.z ;

movss xmm0, dword ptr [esi]this.normal.y ;

mulss xmm0, dword ptr ray_src.y ;

addss xmm7, xmm0 ;

movss xmm0, dword ptr [esi]this.normal.x ;

mulss xmm0, dword ptr ray_src.x ;

addss xmm7, xmm0 ;

addss xmm7, dword ptr [esi]this.normal.w ;

; 163 : t_near = ray->t_near;

movss xmm0, dword ptr [ebx]e_Ray.t_near ;

; 165 : if( t <= t_near * dot_nd )

mulss xmm0, xmm1 ;

comiss xmm0, xmm7 ;

movss dword ptr t, xmm7 ;

; 166 : return false;

jae RETURN_FALSE ;

; 168 : t_far = MIN( ray->t_far, state->t );

movss xmm0, dword ptr [edi]e_RayState.t ;

comiss xmm0, dword ptr [ebx]e_Ray.t_far ;

jbe CASE_BELOW ;

movss xmm0, dword ptr [ebx]e_Ray.t_far ;

//jmp FAR_TEST ;

CASE_BELOW:

//movss xmm0, dword ptr [edi]e_RayState.t ;

//FAR_TEST:

; 169 :

; 170 : if( t >= t_far * dot_nd )

mulss xmm0, xmm1 ;

comiss xmm7, xmm0 ;

; 171 : return false;

jae RETURN_FALSE ;

; 177 : hit.arr[ projX ] = ray_src.arr[ projX ] * dot_nd + ray_dir.arr[ projX ] * t;

movzx eax, byte ptr [esi]this.projX ;

add eax, eax ;

movss xmm0, dword ptr [ray_dir+eax] ;

movss xmm2, dword ptr [ray_src+eax] ;

mulss xmm0, xmm7 ;

mulss xmm2, xmm1 ;

addss xmm0, xmm2 ;

movss dword ptr [hit+eax], xmm0 ;

; 178 : hit.arr[ projY ] = ray_src.arr[ projY ] * dot_nd + ray_dir.arr[ projY ] * t;

movzx ecx, byte ptr [esi]this.projY ;

add ecx, ecx ;

; 182 : if( inited != TRUE )

//cmp byte ptr [esi]this.inited, 1 ;

movss xmm0, dword ptr [ray_dir+ecx] ;

movss xmm2, dword ptr [ray_src+ecx] ;

mulss xmm0, xmm7 ;

mulss xmm2, xmm1 ;

addss xmm0, xmm2 ;

movss dword ptr [hit+ecx], xmm0 ;

//je BARY_TEST ;

; 183 : prepare_intersect();

//call prepare_intersect ;

//BARY_TEST:

; 187 : bary.x = hit.arr[ projX ] * la.x + hit.arr[ projY ] * la.y + la.z * dot_nd;

movss xmm6, dword ptr [esi]this.la.y ;

mulss xmm6, dword ptr [hit+ecx] ;

movss xmm0, dword ptr [esi]this.la.z ;

lea ecx, dword ptr [hit+ecx] ;

mulss xmm0, xmm1 ;

addss xmm6, xmm0 ;

movss xmm0, dword ptr [esi]this.la.x ;

mulss xmm0, dword ptr [hit+eax] ;

lea eax, dword ptr [hit+eax] ;

addss xmm6, xmm0 ;

; 189 : if( bary.x < 0.0f || bary.x > dot_nd )

comiss xmm3, xmm6 ;

ja RETURN_FALSE ;

comiss xmm6, xmm1 ;

ja RETURN_FALSE ;

; 192 : bary.y = hit.arr[ projX ] * lb.x + hit.arr[ projY ] * lb.y + lb.z * dot_nd;

movss xmm5, dword ptr [esi]this.lb.y ;

mulss xmm5, dword ptr [ecx] ;

movss xmm0, dword ptr [esi]this.lb.z ;

mulss xmm0, xmm1 ;

addss xmm5, xmm0 ;

movss xmm0, dword ptr [esi]this.lb.x ;

mulss xmm0, dword ptr [eax] ;

addss xmm5, xmm0 ;

; 194 : if( bary.y < 0.0f || bary.y > dot_nd )

comiss xmm3, xmm5 ;

ja RETURN_FALSE ;

comiss xmm5, xmm1 ;

ja RETURN_FALSE ;

; 197 : bary.z = dot_nd - bary.x - bary.y;

movaps xmm4, xmm1 ;

subss xmm4, xmm6 ;

subss xmm4, xmm5 ;

; 199 : if( bary.z < 0.0f || bary.z > dot_nd )

comiss xmm3, xmm4 ;

ja RETURN_FALSE ;

comiss xmm4, xmm1 ;

ja RETURN_FALSE ;

movss dword ptr bary.x, xmm6 ;

movss dword ptr bary.y, xmm5 ;

movss dword ptr bary.z, xmm4 ;

} // __asm