//
//  pd3d.h
//  Playdate Simulator
//
//  Created by Dave Hayden on 8/25/15.
//  Copyright (c) 2015 Panic, Inc. All rights reserved.
//

#ifndef __Playdate_Simulator__pd3d__
#define __Playdate_Simulator__pd3d__

#include <stdio.h>
#include <stdint.h>
#include <math.h>

#if TARGET_PLAYDATE
static inline float _lib3d_sqrtf(float x)
{
	float result;
	asm ("vsqrt.f32 %0, %1" : "=t"(result) : "t"(x));
	return result;
}
#define sqrtf(f) _lib3d_sqrtf(f)
#endif

typedef struct
{
	float x;
	float y;
	float z;
} Point3D;

static inline Point3D Point3DMake(float x, float y, float z)
{
	return (Point3D){ .x = x, .y = y, .z = z };
}

static inline int Point3D_equals(Point3D a, Point3D b)
{
	return (a.x == b.x) && (a.y == b.y) && (a.z == b.z);
}

typedef struct
{
	float dx;
	float dy;
	float dz;
} Vector3D;

static inline Vector3D Vector3DMake(float dx, float dy, float dz)
{
	return (Vector3D){ .dx = dx, .dy = dy, .dz = dz };
}

static inline Vector3D Vector3DCross(Vector3D a, Vector3D b)
{
	return (Vector3D){ .dx = a.dy * b.dz - a.dz * b.dy, .dy = a.dz * b.dx - a.dx * b.dz, .dz = a.dx * b.dy - a.dy * b.dx };
}

static inline float Vector3DDot(Vector3D a, Vector3D b)
{
	return a.dx * b.dx + a.dy * b.dy + a.dz * b.dz;
}

Vector3D pnormal(Point3D* p1, Point3D* p2, Point3D* p3);

static inline float Vector3D_lengthSquared(Vector3D* v)
{
	return v->dx * v->dx + v->dy * v->dy + v->dz * v->dz;
}

static inline float Vector3D_length(Vector3D* v)
{
	return sqrtf(Vector3D_lengthSquared(v));
}

static inline Point3D Point3D_addVector(Point3D a, Vector3D v)
{
	return Point3DMake(a.x + v.dx, a.y + v.dy, a.z + v.dz);
}

// XXX - which side are we multiplying on?

typedef struct
{
	int isIdentity : 1;
	int inverting : 1;
	float m[3][3];
	float dx;
	float dy;
	float dz;
} Matrix3D;

extern Matrix3D identityMatrix;

Matrix3D Matrix3DMake(float m11, float m12, float m13, float m21, float m22, float m23, float m31, float m32, float m33, int inverting);
Matrix3D Matrix3DMakeTranslate(float dx, float dy, float dz);
Matrix3D Matrix3D_multiply(Matrix3D l, Matrix3D r);
Point3D Matrix3D_apply(Matrix3D m, Point3D p);

float Matrix3D_getDeterminant(Matrix3D* m);

#endif