vectors_3d.c File Reference

Functions related to 3D vector operations. More...

#include <stdio.h>
#include <math.h>
#include <assert.h>
#include "geometry_datatypes.h"

Include dependency graph for vectors_3d.c:

Functions
vec_3d	vector_sub (const vec_3d *a, const vec_3d *b)
	Subtract one vector from another.
vec_3d	vector_add (const vec_3d *a, const vec_3d *b)
	Add one vector to another.
float	dot_prod (const vec_3d *a, const vec_3d *b)
	Obtain the dot product of two 3D vectors.
vec_3d	vector_prod (const vec_3d *a, const vec_3d *b)
	Compute the vector product of two 3d vectors.
const char *	print_vector (const vec_3d *a, const char *name)
	Print formatted vector on stdout.
float	vector_norm (const vec_3d *a)
	Compute the norm a vector.
vec_3d	unit_vec (const vec_3d *a)
	Obtain unit vector in the same direction as given vector.
mat_3x3	get_cross_matrix (const vec_3d *a)
	The cross product of vectors can be represented as a matrix multiplication operation.
double	get_angle (const vec_3d *a, const vec_3d *b)
	Obtain the angle between two given vectors.
static void	test ()
	Testing function.
int	main (void)
	Main function.

Detailed Description

Functions related to 3D vector operations.

Author

Krishna Vedala

Function Documentation

main()

int main

(

void

)

Main function.

Returns

0 on exit

{
    test();
 
    return 0;
}

Here is the call graph for this function:

test()

static void test ( void  )

static

Testing function.

Returns

void

{
    vec_3d a = {1., 2., 3.};
    vec_3d b = {1., 1., 1.};
    float d;
 
    // printf("%s", print_vector(&a, "a"));
    // printf("%s", print_vector(&b, "b"));
 
    d = vector_norm(&a);
    // printf("|a| = %.4g\n", d);
    assert(fabsf(d - 3.742f) < 0.01);
    d = vector_norm(&b);
    // printf("|b| = %.4g\n", d);
    assert(fabsf(d - 1.732f) < 0.01);
 
    d = dot_prod(&a, &b);
    // printf("Dot product: %f\n", d);
    assert(fabsf(d - 6.f) < 0.01);
 
    vec_3d c = vector_prod(&a, &b);
    // printf("Vector product ");
    // printf("%s", print_vector(&c, "c"));
    assert(fabsf(c.x - (-1.f)) < 0.01);
    assert(fabsf(c.y - (2.f)) < 0.01);
    assert(fabsf(c.z - (-1.f)) < 0.01);
 
    double alpha = get_angle(&a, &b);
    // printf("The angle is %f\n", alpha);
    assert(fabsf(alpha - 0.387597) < 0.01);
}

Here is the call graph for this function: