mirror of
https://git.code.sf.net/p/quake/quakeforge
synced 2024-11-15 01:11:27 +00:00
b9cff7aae0
This gives the resultant point the correct sign. Though the projective divide would take care of the sign, this makes reading the point a little less confusing (still need to sort out automatic blade reversals for the likes of e31).
98 lines
2.4 KiB
R
98 lines
2.4 KiB
R
#include <AutoreleasePool.h>
|
|
#include "algebra.h"
|
|
#include "basisblade.h"
|
|
#include "basisgroup.h"
|
|
#include "metric.h"
|
|
#include "multivector.h"
|
|
#include "util.h"
|
|
|
|
@static AutoreleasePool *autorelease_pool;
|
|
@static void
|
|
arp_start (void)
|
|
{
|
|
autorelease_pool = [[AutoreleasePool alloc] init];
|
|
}
|
|
|
|
@static void
|
|
arp_end (void)
|
|
{
|
|
[autorelease_pool release];
|
|
autorelease_pool = nil;
|
|
}
|
|
|
|
int
|
|
main ()
|
|
{
|
|
arp_start ();
|
|
|
|
BasisBlade *a = [[BasisBlade basis:1] retain];
|
|
BasisBlade *b = [[BasisBlade basis:2] retain];
|
|
BasisBlade *c = [[BasisBlade basis:4] retain];
|
|
BasisBlade *d = [[BasisBlade basis:8] retain];
|
|
BasisBlade *blades[] = {a, b, c, d};
|
|
static string names[] = {"a", "b", "c", "d"};
|
|
|
|
// printf ("a: %@\n", a);
|
|
// printf ("b: %@\n", b);
|
|
// printf ("c: %@\n", c);
|
|
// printf ("d: %@\n", d);
|
|
|
|
arp_end ();
|
|
#if 0
|
|
arp_start ();
|
|
for (int i = 0; i < 4; i++) {
|
|
arp_end ();
|
|
arp_start ();
|
|
BasisBlade *vec = blades[i];
|
|
printf ("%s: %@\n", names[i], vec);
|
|
for (int j = 0; j < 4; j++) {
|
|
BasisBlade *bvec = [vec outerProduct:blades[j]];
|
|
if (![bvec scale]) {
|
|
continue;
|
|
}
|
|
printf ("%s^%s: %@\n", names[i], names[j], bvec);
|
|
for (int k = 0; k < 4; k++) {
|
|
BasisBlade *tvec = [bvec outerProduct:blades[k]];
|
|
if (![tvec scale]) {
|
|
continue;
|
|
}
|
|
printf ("%s^%s^%s: %@\n", names[i], names[j], names[k],
|
|
tvec);
|
|
for (int l = 0; l < 4; l++) {
|
|
BasisBlade *qvec = [tvec outerProduct:blades[l]];
|
|
if (![qvec scale]) {
|
|
continue;
|
|
}
|
|
printf ("%s^%s^%s^%s: %@\n",
|
|
names[i], names[j], names[k], names[l],
|
|
qvec);
|
|
}
|
|
}
|
|
}
|
|
}
|
|
arp_end ();
|
|
#endif
|
|
arp_start ();
|
|
|
|
Metric *m = [Metric R:3,0,1];
|
|
BasisBlade *ad = [a geometricProduct:d metric:m];
|
|
BasisBlade *prod = [ad geometricProduct:ad metric:m];
|
|
printf ("%s%s %s%s: %@\n",
|
|
names[0], names[3], names[0], names[3], prod);
|
|
|
|
Algebra *alg = [Algebra R:3, 0, 1];
|
|
double plane1_vals[4] = {1, 0, 0, 8};
|
|
double plane2_vals[4] = {0, 1, 0, 8};
|
|
double origin_vals[4] = {0, 0, 0, 1};
|
|
MultiVector *plane1 = [alg group:0 values:plane1_vals];
|
|
MultiVector *plane2 = [alg group:0 values:plane2_vals];
|
|
MultiVector *origin = [alg group:3 values:origin_vals];
|
|
|
|
MultiVector *line = [plane1 wedge:plane2];
|
|
MultiVector *point = [[line dot:origin] product:[line reverse]];
|
|
printf ("plane1:%@\nplane2:%@\nline:%@\norigin:%@\n", plane1, plane2, line, origin);
|
|
printf ("point:%@\n", point);
|
|
|
|
arp_end ();
|
|
return 0;
|
|
}
|