physics.h

Go to the documentation of this file.

#ifndef BJ_PHYSICS_H
#define BJ_PHYSICS_H
 
#include <banjo/math.h>
#include <banjo/vec.h>


#define BJ_GRAVITATIONAL_CONSTANT_SI BJ_F(6.67430e-11)
 

static BJ_INLINE bj_real bj_galileo_position(
    bj_real position,
    bj_real velocity,
    bj_real acceleration,
    bj_real time
) {
    return BJ_F(0.5) * acceleration * time * time + velocity * time + position;
}

static BJ_INLINE bj_real bj_galileo_velocity(
    bj_real velocity,
    bj_real acceleration,
    bj_real time
) {
    return acceleration * time + velocity;
}

static BJ_INLINE bj_real bj_newton_gravitation(
    bj_real m1,
    bj_real m2,
    bj_real r,
    bj_real g
) {
    return g * (m1 * m2) / (r * r);
}

static BJ_INLINE bj_real bj_newton_plummer_gravitation(
    bj_real m1,
    bj_real m2,
    bj_real r,
    bj_real g,
    bj_real eps
) {
    const bj_real r2    = r * r;
    const bj_real e2    = eps * eps;
    const bj_real denom = bj_pow(r2 + e2, BJ_F(1.5));
    return (denom > BJ_FZERO) ? (g * m1 * m2 * r) / denom : BJ_FZERO;
}
 

 
#endif