/**
* graphene_matrix_transform_vec3:
* @m: a #graphene_matrix_t
* @v: a #graphene_vec3_t
* @res: (out caller-allocates): return location for a #graphene_vec3_t
*
* Transforms the given #graphene_vec3_t using the matrix @m.
*
* Since: 1.0
*/
void
graphene_matrix_transform_vec3 (const graphene_matrix_t *m,
const graphene_vec3_t *v,
graphene_vec3_t *res)
{
graphene_simd4x4f_vec3_mul (&m->value, &v->value, &res->value);
}
/**
* graphene_matrix_transform_vec3:
* @m: ...
* @v: ...
* @res: (out caller-allocates): ...
*
* ...
*
* Since: 1.0
*/
void
graphene_matrix_transform_vec3 (const graphene_matrix_t *m,
const graphene_vec3_t *v,
graphene_vec3_t *res)
{
g_return_if_fail (m != NULL && v != NULL);
g_return_if_fail (res != NULL);
graphene_simd4x4f_vec3_mul (&m->value, &v->value, &res->value);
}
static void
matrix_project (gpointer data_)
{
MatrixBench *data = data_;
int i;
for (i = 0; i < N_ROUNDS; i++)
{
graphene_simd4f_t pback, qback, uback;
float t, x, y;
graphene_simd4x4f_vec3_mul (&(data->a[i]), &(data->pa[i]), &pback);
graphene_simd4x4f_vec3_mul (&(data->a[i]), &(data->qa[i]), &qback);
uback = graphene_simd4f_sub (data->pa[i], pback);
t = -1.0f * graphene_simd4f_get_z (pback) / graphene_simd4f_get_z (uback);
x = graphene_simd4f_get_x (pback) + t * graphene_simd4f_get_x (uback);
y = graphene_simd4f_get_y (pback) + t * graphene_simd4f_get_y (uback);
data->ra[i] = graphene_simd4f_init (x, y, 0.f, 0.f);
}
}
/**
* graphene_matrix_transform_point:
* @m: a #graphene_matrix_t
* @p: a #graphene_point_t
* @res: (out caller-allocates): return location for the
* transformed #graphene_point_t
*
* Transforms the given #graphene_point_t using the matrix @m.
*
* Since: 1.0
*/
void
graphene_matrix_transform_point (const graphene_matrix_t *m,
const graphene_point_t *p,
graphene_point_t *res)
{
graphene_simd4f_t vec3;
vec3 = graphene_simd4f_init (p->x, p->y, 0.0f, 0.0f);
graphene_simd4x4f_vec3_mul (&m->value, &vec3, &vec3);
res->x = graphene_simd4f_get_x (vec3);
res->y = graphene_simd4f_get_y (vec3);
}
