qf/quakeforge
0
0
Fork 0
mirror of https://git.code.sf.net/p/quake/quakeforge synced 2025-02-25 21:21:14 +00:00
Code Issues Projects Releases Packages Wiki Activity

[qfcc] Optimize (a+a)/2 in geometric algebra expressions

Browse source 
Really, I need to look into getting the optimizations more general, but
`a × ⋆a` is now just a single cross product instead off cross, add and
divide.
This commit is contained in:
Bill Currie 2024-09-30 19:11:40 +09:00
parent b00d4d3b5b
commit 9797dcb5de
1 changed files with 24 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

18
tools/qfcc/source/expr_algebra.c
View file

@ -2634,6 +2634,18 @@ commutator_product (const expr_t *e1, const expr_t *e2)
new_int_expr (2, false)); new_int_expr (2, false));
} }
static bool
is_two (const expr_t *e)
{
if (is_integral_val (e)) {
return expr_integral (e) == 2;
}
if (is_floating_val (e)) {
return expr_floating (e) == 2;
}
return false;
}
static const expr_t * static const expr_t *
multivector_divide (const expr_t *e1, const expr_t *e2) multivector_divide (const expr_t *e1, const expr_t *e2)
{ {
@ -2650,6 +2662,7 @@ multivector_divide (const expr_t *e1, const expr_t *e2)
auto layout = &algebra->layout; auto layout = &algebra->layout;
auto stype = algebra->type; auto stype = algebra->type;
const expr_t *a[layout->count] = {}; const expr_t *a[layout->count] = {};
bool is_half = is_two (e2);
e1 = mvec_expr (e1, algebra); e1 = mvec_expr (e1, algebra);
e2 = promote_scalar (algebra->type, e2); e2 = promote_scalar (algebra->type, e2);
mvec_scatter (a, e1, algebra); mvec_scatter (a, e1, algebra);
@ -2659,6 +2672,10 @@ multivector_divide (const expr_t *e1, const expr_t *e2)
continue; continue;
} }
auto den = e2; auto den = e2;
if (is_half && a[i]->type == ex_expr && a[i]->expr.op == '+'
&& a[i]->expr.e1 == a[i]->expr.e2) {
a[i] = a[i]->expr.e1;
} else {
auto ct = get_type (a[i]); auto ct = get_type (a[i]);
int width = type_width (ct); int width = type_width (ct);
if (width > 1) { if (width > 1) {
@ -2667,6 +2684,7 @@ multivector_divide (const expr_t *e1, const expr_t *e2)
a[i] = typed_binary_expr (ct, '/', a[i], den); a[i] = typed_binary_expr (ct, '/', a[i], den);
a[i] = edag_add_expr (a[i]); a[i] = edag_add_expr (a[i]);
} }
}
return mvec_gather (a, algebra); return mvec_gather (a, algebra);
} }