 |
Robotran C Documentation
|
Go to the documentation of this file.
29 int mbs_eig_0(
double **A,
int n,
double *eval_a,
double *eval_b,
double **evec_a,
double **evec_b);
43 int mbs_schur_0(
double **A,
int n,
double **T,
double **Z,
double *WR,
double *WI,
int reorder);
107 int mbs_svdDcmp_0(
double **A,
int x,
int y,
double **U,
double **S,
double **VT);
int w_methods(double y[], double dydx[], int n, double x, double h, double yout[], int(*derivs)(double, double[], double[], MbsData *, MbsDirdyn *), MbsData *s, MbsDirdyn *mbs_dd)
W Methods integrator implementation.
Definition: w_methods.c:27
#define GAM
Definition: rosenbrock.c:45
#define SAFETY
Definition: rosenbrock.c:36
int mbs_estim_jac_acc(double x, double htry, double y[], double dydx[], int compute_dfdx, double dfdx[], double **dfdy, int n, int(*derivs)(double, double[], double[], MbsData *, MbsDirdyn *), MbsData *s, MbsDirdyn *dd)
Evaluate the Jacobian of a function accelerations using finite difference.
Definition: integrator.c:30
double * qd
Array with the current values of the generalized velocities.
Definition: mbs_data.h:188
#define E4
Definition: rosenbrock.c:62
int mbs_cholsl_0(double **A, int n, double **B, int nb)
Solves a symmetric positive definite system of linear equations AX=B (B being the right hand side wit...
Definition: cholesky.c:40
#define C3X
Definition: rosenbrock.c:65
This header defines global functions for W methods integrator in C.
void lubksb_0(double **a, int n, int *indx, double b[])
Definition: lubksb.c:31
#define B1
Definition: rosenbrock.c:55
double tsim
The time value.
Definition: mbs_data.h:208
double * y
Pointer to the state vector, this is a dvec_0 of size nState.
Definition: mbs_dirdyn_struct.h:153
#define GROW
Definition: rosenbrock.c:37
#define W_S
Definition: mbs_w_methods.h:27
int nState
size of the state vector
Definition: mbs_dirdyn_struct.h:152
#define C32
Definition: rosenbrock.c:51
#define ERRCON
Definition: rosenbrock.c:41
int nqu
Number of independent articulations.
Definition: mbs_data.h:165
#define C43
Definition: rosenbrock.c:54
#define A32
Definition: rosenbrock.c:48
#define E1
Definition: rosenbrock.c:59
#define E2
Definition: rosenbrock.c:60
#define A21
Definition: rosenbrock.c:46
#define PSHRNK
Definition: rosenbrock.c:40
This header defines global functions for rosenbrock integrator in C.
int mbs_invLU_0(double **A, double **Am1, int n)
Inverse the square matrix A using a LU decomposition.
Definition: inverse_lapack.c:28
int mbs_freeze_jac(int freeze_index, int *next_freeze_index, double x, double h, double y[], double dydx[], int compute_dfdx, double dfdx[], double **dfdy, double dydx_freeze[], double dfdx_freeze[], double **dfdy_freeze, int n, int(*derivs)(double, double[], double[], MbsData *, MbsDirdyn *), MbsData *s, MbsDirdyn *mbs_dd)
Freeze the Jacobian of an integrator structure until mbs_dd->options->n_freeze.
Definition: integrator.c:107
void mbs_close_kin(MbsData *s, MbsAux *mbs_aux)
compute the dependent velocities that solves the constraints.
Definition: mbs_close_loops.c:140
void * integrator_struct
pointer to store integrator structure
Definition: mbs_dirdyn_struct.h:183
#define A3X
Definition: rosenbrock.c:68
#define C31
Definition: rosenbrock.c:50
#define C1X
Definition: rosenbrock.c:63
int mbs_svdDcmp_0(double **A, int x, int y, double **U, double **S, double **VT)
Realize the singular value decomposition.
Definition: svdcmp_lapack.c:78
int ludcmp_0(double **a, int n, int *indx, double *d)
LU decomposition of a matrix with index starting at 0.
Definition: ludcmp.c:82
double tsim
current simulation time
Definition: mbs_dirdyn_struct.h:149
void mbs_msg(char *msg,...)
Send a message.
Definition: mbs_message.c:76
Runge-Kutta 4 structure for dirdyn.
Definition: mbs_rk4.h:31
Rosenbrock integration structure for dirdyn.
Definition: mbs_rosenbrock.h:29
int * qu
Array with the indices of the independent articulations (free)
Definition: mbs_data.h:174
#define B2
Definition: rosenbrock.c:56
#define C4X
Definition: rosenbrock.c:66
int n_freeze
number of time step when the jacobian is freezed (computed once at the start of the n_freeze time ste...
Definition: mbs_dirdyn_struct.h:134
int mbs_eig_0(double **A, int n, double *eval_a, double *eval_b, double **evec_a, double **evec_b)
Compute the eigen values and vector of a matrix with index starting a 0.
Definition: eig_lapack.c:32
#define B3
Definition: rosenbrock.c:57
double * q
Array with the current values of the generalized coordinates.
Definition: mbs_data.h:187
int rosenbrock(int n, int(*derivs)(double, double[], double[], MbsData *, MbsDirdyn *), double *x, double y[], double eps, double hmax, double htry, long nmax, double dydx[], double yscal[], double *hnext, MbsData *s, MbsDirdyn *dd, double *hdid)
Fourth-order Rosenbrock step for integrating stiff problems, with monitoring of local truncation erro...
Definition: rosenbrock.c:103
#define E3
Definition: rosenbrock.c:61
int rk4(double y[], double dydx[], int n, double x, double h, double yout[], int(*derivs)(double, double[], double[], MbsData *, MbsDirdyn *), MbsData *s, MbsDirdyn *dd)
Runge Kutta 4 integrator implementation Given values for the variables y[1..n] and their derivatives ...
Definition: rk4.c:21
W Methods structure for dirdyn.
Definition: mbs_w_methods.h:52
#define _MBS_ERR_MOD_SPEC_14
Generic error number Module specific errors range from -11 to -19 please read the error message a...
Definition: mbs_errors_names.h:180
int mbs_close_geo(MbsData *s, MbsAux *mbs_aux)
compute a position of the multibody system which solves the constraints
Definition: mbs_close_loops.c:43
static double dmaxarg2
Definition: rosenbrock.c:100
#define _MBS_ERR_MOD_SPEC_18
Genericerror number Module specific errors range from -11 to -19 please read the error message an...
Definition: mbs_errors_names.h:204
Definition: mbs_data.h:149
MbsAux * mbs_aux
Definition: mbs_dirdyn_struct.h:147
int mbs_choldc_0(double **A, int n)
Compute the Cholesky decomposition of A, i.e.
Definition: cholesky.c:25
#define A31
Definition: rosenbrock.c:47
int dirdynared(MbsAux *mbs_aux, MbsData *s)
Definition: mbs_dirdynared.c:25
int ThetaSC(double y[], double dydx[], int n, double x, double h, double yout[], int(*derivs)(double, double[], double[], MbsData *, MbsDirdyn *), MbsData *s, MbsDirdyn *mbs_dd)
Unknown integrator.
Definition: thetaSC.c:26
static double dmaxarg1
Definition: rosenbrock.c:100
#define C2X
Definition: rosenbrock.c:64
This header defines specific integrators functions in C.
static double dminarg1
Definition: rosenbrock.c:101
int mbs_over_under_determined(double **A, double *b, int x, int y)
Solve a linear system of equations, A*x=b, using dgelss lapack function.
Definition: svdcmp_lapack.c:18
#define A2X
Definition: rosenbrock.c:67
#define SHRNK
Definition: rosenbrock.c:39
int mbs_schur_0(double **A, int n, double **T, double **Z, double *WR, double *WI, int reorder)
A : matrix for the Eigenvalue problem in row major n : size of the matrix WR : Result of the Eigenval...
Definition: schur_lapack.c:45
#define C41
Definition: rosenbrock.c:52
void copy_dvec_0(double *vec_src, double *vec_dest, int l_v)
Copy the content of a vector (of sizel_v) of doubles to a second vector, starting with index 0.
Definition: mbs_1D_array.c:906
#define C42
Definition: rosenbrock.c:53
General structure of the direct dynamic module.
Definition: mbs_dirdyn_struct.h:144
This header defines global functions for Runge Kutta 4 integrator in C.
double * qdd
Array with the current values of the generalized acceleration.
Definition: mbs_data.h:189
#define C21
Definition: rosenbrock.c:49
int mbs_rank_0(double **A, int m, int n)
Compute the rank of matrix_in.
Definition: rank_lapack.c:38
#define B4
Definition: rosenbrock.c:58
MbsDirdynOptions * options
structure defining the option of a direct dynamic
Definition: mbs_dirdyn_struct.h:146
static double dminarg2
Definition: rosenbrock.c:101
int nhu
Number of independent constraints.
Definition: mbs_data.h:171
#define PGROW
Definition: rosenbrock.c:38
1.8.17