Skip to content

WIP: SSE Funhouse #29

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
sbryant wants to merge 5 commits into
base: master
Choose a base branch
from
Open

WIP: SSE Funhouse #29

Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
0e7eef9
SSE dot product anyone?
sbryant Jun 12, 2013
1331ecf
This is can be an aligned load!
sbryant Jun 12, 2013
ed51e49
Conditionally compile SSE f3.
sbryant Jun 12, 2013
c7e94f0
Just worry about the function body.
sbryant Jun 12, 2013
0d6b32a
it's actually __SSE4_1__
sbryant Jun 12, 2013
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
25 changes: 22 additions & 3 deletions src/vector.c
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,7 @@
#include "vector.h"
#if defined __SSE4_1__
#include <smmintrin.h>
#endif

float3 *clone_f3(float3 f) {
float3 *clone = malloc(sizeof(float3));
Expand Down Expand Up @@ -27,12 +30,28 @@ float3 *f3_scale(float3 *v, float c) {
return v;
}


float f3_dot(float3 v1, float3 v2) {
return (v1[0] * v2[0] +
v1[1] * v2[1] +
v1[2] * v2[2]);
#if defined __SSE4_1__
float result;
float arr[4] = { v1[0], v1[1], v1[2], 0.0 };
float arr2[4] = { v2[0], v2[1], v2[2], 0.0 };

__m128 a = _mm_load_ps(arr);
__m128 b = _mm_load_ps(arr2);
const int mask = 0xff;
__m128 sse_result = _mm_dp_ps(a, b, mask);

_mm_store_ss(&result, sse_result);
return result;
#else
return (v1[0] * v2[0] +
v1[1] * v2[1] +
v1[2] * v2[2]);
#endif
}


float3 *f3_cross(float3 *result, float3 v1, float3 v2) {
(*result)[0] = v1[1]*v2[2] - v1[2]*v2[1];
(*result)[1] = v1[2]*v2[0] - v1[0]*v2[2];
Expand Down