mbs_euler_implicit.c File Reference

This file implements the functions of the Euler implicit integration method in C. More...

#include <math.h>
#include "integrator.h"
#include "mbs_matrix.h"
#include "mbs_message.h"
#include "mbs_1D_array.h"
#include "mbs_dirdyn_struct.h"
#include "mbs_euler_implicit.h"

Functions
void	initialize_euler_implicit (MbsData *mbs_data, MbsDirdyn *mbs_dd, int(*derivs)(double, double[], double[], MbsData *, MbsDirdyn *), void(*loop_fct)(MbsData *, MbsDirdyn *), int(*save_realtime_update)(MbsDirdyn *dd, MbsData *mbs_data))
	Initialize all needed fields in MbsDirdyn for Euler Implicit integrator. More...
int	loop_euler_implicit (double t0, double tf, MbsData *mbs_data, MbsDirdyn *mbs_dd)
	Integrate the motion from t0 to tf. More...
void	finish_euler_implicit (MbsData *mbs_data, MbsDirdyn *dd)
	Finalize and clear memory allocated for Euler Implicit integrator. More...
void	error_euler_implicit (MbsData *mbs_data, MbsDirdyn *dd, int err)
	If an error occurs during time integration, this function prints the error message and call finish_euler_implicit. More...

Detailed Description

This file implements the functions of the Euler implicit integration method in C.

Creation date: April 2018

Author

Sebastien Timmermans

@source Arnold M. et al, Linearly implicit time integration methods in real-time applications: DAEs and stiff ODEs Multibody System Dynamics, 2007, 17:99–117

(c) Universite catholique de Louvain

Function Documentation

error_euler_implicit()

void error_euler_implicit	(	MbsData *	mbs_data,
		MbsDirdyn *	dd,
		int	err
	)

If an error occurs during time integration, this function prints the error message and call finish_euler_implicit.

Parameters

mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure
err	The error number

finish_euler_implicit()

void finish_euler_implicit	(	MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd
	)

Finalize and clear memory allocated for Euler Implicit integrator.

Parameters

mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure

initialize_euler_implicit()

void initialize_euler_implicit	(	MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd,
		int(*)(double, double[], double[], MbsData *, MbsDirdyn *)	derivs,
		void(*)(MbsData *, MbsDirdyn *)	loop_fct,
		int(*)(MbsDirdyn *dd, MbsData *mbs_data)	save_realtime_update
	)

Initialize all needed fields in MbsDirdyn for Euler Implicit integrator.

Parameters

mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure
derivs	The function computing f'
loop_fct	The function to be called at the end of each time step, or waypoints.

loop_euler_implicit()

int loop_euler_implicit	(	double	t0,
		double	tf,
		MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd
	)

Integrate the motion from t0 to tf.

Parameters

t0	The initial time
tf	The final time
mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure

Returns

Error status, <0 in case of failure.

user loop

user loop