Files
create/src/Mod/Robot/App/kdl_cp/chainiksolverpos_nr_jl.cpp
2011-10-10 13:44:52 +00:00

87 lines
2.8 KiB
C++

// Copyright (C) 2007 Ruben Smits <ruben dot smits at mech dot kuleuven dot be>
// Copyright (C) 2008 Mikael Mayer
// Copyright (C) 2008 Julia Jesse
// Version: 1.0
// Author: Ruben Smits <ruben dot smits at mech dot kuleuven dot be>
// Maintainer: Ruben Smits <ruben dot smits at mech dot kuleuven dot be>
// URL: http://www.orocos.org/kdl
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
// Lesser General Public License for more details.
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
#include "chainiksolverpos_nr_jl.hpp"
#ifndef M_PI
#define M_PI 3.14159265358979323846 /* pi */
#endif
#ifndef M_PI_2
#define M_PI_2 1.57079632679489661923 /* pi/2 */
#endif
namespace KDL
{
ChainIkSolverPos_NR_JL::ChainIkSolverPos_NR_JL(const Chain& _chain, const JntArray& _q_min, const JntArray& _q_max, ChainFkSolverPos& _fksolver,ChainIkSolverVel& _iksolver,
unsigned int _maxiter, double _eps):
chain(_chain), q_min(chain.getNrOfJoints()), q_max(chain.getNrOfJoints()), fksolver(_fksolver),iksolver(_iksolver),delta_q(_chain.getNrOfJoints()),
maxiter(_maxiter),eps(_eps)
{
q_min = _q_min;
q_max = _q_max;
}
int ChainIkSolverPos_NR_JL::CartToJnt(const JntArray& q_init, const Frame& p_in, JntArray& q_out)
{
q_out = q_init;
unsigned int i;
for(i=0;i<maxiter;i++){
fksolver.JntToCart(q_out,f);
delta_twist = diff(f,p_in);
if(Equal(delta_twist,Twist::Zero(),eps))
break;
iksolver.CartToJnt(q_out,delta_twist,delta_q);
Add(q_out,delta_q,q_out);
for(unsigned int j=0; j<q_min.rows(); j++) {
if(q_out(j) < q_min(j))
//q_out(j) = q_min(j);
q_out(j) = q_out(j) + M_PI *2;
}
for(unsigned int j=0; j<q_max.rows(); j++) {
if(q_out(j) > q_max(j))
//q_out(j) = q_max(j);
q_out(j) = q_out(j) - M_PI *2;
}
}
if(i!=maxiter)
return 0;
else
return -3;
}
ChainIkSolverPos_NR_JL::~ChainIkSolverPos_NR_JL()
{
}
}