void time_loop(const int active_source) { for (int i=0;i<max_num_iter;i++) { velocity_update(i,active_source); displacement_update(); record(i,active_source); } }
void velocity_init(struct element *e) { struct element *p = e->parent; /* Basis velocity, in parent coordinates */ cross_product_mat4(&e->a, &e->tv); print("angular velocity", &e->a); print("velocity", &e->tv); velocity_update(e); }