OpenXM/src/hgm/so3/src/so3_nc.c - diff

Return to so3_nc.c CVS log

Up to [local] / OpenXM / src / hgm / so3 / src

version 1.4, 2013/02/12 23:51:51

version 1.5, 2013/03/05 05:26:07

Line 1

/* $OpenXM$ */

/* $OpenXM: OpenXM/src/hgm/so3/src/so3_nc.c,v 1.4 2013/02/12 23:51:51 takayama Exp $ */

#include <stdio.h>

#include <math.h>

#ifdef USE_GSL_LIB

Line 183 Ref: Note: See the Corollary 1 of the paper rotation.c

int so3_func(double t, const double y[], double f[], void *params)

{

extern double SO3_A[3];

#ifndef STANDALONE

mh_check_intr(100);

#endif

f[0] = SO3_A[0]*y[1]+SO3_A[1]*y[2]+SO3_A[2]*y[3] -SO3_R*y[0];

f[1] = SO3_A[0]*y[0]+SO3_A[2]*y[2]+SO3_A[1]*y[3] - (2/t)*y[1] -SO3_R*y[1];

f[2] = SO3_A[1]*y[0]+SO3_A[2]*y[1]+SO3_A[0]*y[3] - (2/t)*y[2] -SO3_R*y[2];