Program Listing for File ConstraintKinematicsSolver.h

Program Listing for File ConstraintKinematicsSolver.h#

Return to documentation for file (include/Karana/SOADyn/ConstraintKinematicsSolver.h)

/*
*Copyright(c)2024-2026KaranaDynamicsPtyLtd.Allrightsreserved.
*
*NOTICETOUSER:
*
*Thissourcecodeand/ordocumentation(the"LicensedMaterials")is
*theconfidentialandproprietaryinformationofKaranaDynamicsInc.
*UseoftheseLicensedMaterialsisgovernedbythetermsandconditions
*ofaseparatesoftwarelicenseagreementbetweenKaranaDynamicsandthe
*Licensee("LicenseAgreement").Unlessexpresslypermittedunderthat
*agreement,anyreproduction,modification,distribution,ordisclosure
*oftheLicensedMaterials,inwholeorinpart,toanythirdparty
*withoutthepriorwrittenconsentofKaranaDynamicsisstrictlyprohibited.
*
*THELICENSEDMATERIALSAREPROVIDED"ASIS"WITHOUTWARRANTYOFANYKIND.
*KARANADYNAMICSDISCLAIMSALLWARRANTIES,EXPRESSORIMPLIED,INCLUDING
*BUTNOTLIMITEDTOWARRANTIESOFMERCHANTABILITY,NON-INFRINGEMENT,AND
*FITNESSFORAPARTICULARPURPOSE.
*
*INNOEVENTSHALLKARANADYNAMICSBELIABLEFORANYDAMAGESWHATSOEVER,
*INCLUDINGBUTNOTLIMITEDTOLOSSOFPROFITS,DATA,ORUSE,EVENIF
*ADVISEDOFTHEPOSSIBILITYOFSUCHDAMAGES,WHETHERINCONTRACT,TORT,
*OROTHERWISEARISINGOUTOFORINCONNECTIONWITHTHELICENSEDMATERIALS.
*
*U.S.GovernmentEndUsers:TheLicensedMaterialsarea"commercialitem"
*asdefinedat48C.F.R.2.101,andareprovidedtotheU.S.Government
*onlyasacommercialenditemunderthetermsofthislicense.
*
*AnyuseoftheLicensedMaterialsinindividualorcommercialsoftwaremust
*include,intheuserdocumentationandinternalsourcecodecomments,
*thisNotice,Disclaimer,andU.S.GovernmentUseProvision.
*/


#pragmaonce

#include<string>
#include<vector>

#include"JacobianGenerator.h"
#include"Karana/KCore/LockingBase.h"
#include"Karana/Math/NonlinearSolver.h"

namespaceKarana::Dynamics{
namespacekc=Karana::Core;
namespacekm=Karana::Math;

template<int,int,int>classCoordBase_T;
classCoordData;
classLoopConstraint;
classSubGraph;

classConstraintKinematicsSolver:publickc::LockingBase{

/*foraccessto_velJacobian()*/
friendclassAlgorithms;

/*foraccesstocutjointjacobiangenerator*/
friendclassCECompoundBody;

public:
statickc::ks_ptr<ConstraintKinematicsSolver>create(std::string_viewname,
kc::ks_ptr<SubGraph>sg);

ConstraintKinematicsSolver(std::string_viewname,kc::ks_ptr<SubGraph>sg);

~ConstraintKinematicsSolver();

//----------------------------

doublesolveQ();

doublesolveU();

doublesolveUdot();

enumclassCKFrozenCoordType{
Q,
U,
Q_U,
};

voidfreezeCoord(kc::ks_ptr<CoordBase>coord,
size_tindex,
CKFrozenCoordTypetype=CKFrozenCoordType::Q_U);

voidunfreezeCoord(kc::ks_ptr<CoordBase>coord,
size_tindex,
CKFrozenCoordTypetype=CKFrozenCoordType::Q_U);

voidclearFrozenCoords();

usingCoordIndices=std::vector<unsignedint>;

CoordIndicesfrozenCoords(boolQ_type)const;//NOLINT(readability-identifier-naming)

km::VecgetNonfrozenQ()const;

voidsetNonfrozenQ(constEigen::Ref<constkm::Vec>&Q);

km::VecgetNonfrozenU()const;

voidsetNonfrozenU(constEigen::Ref<constkm::Vec>&U);

km::VecgetNonfrozenUdot()const;

voidsetNonfrozenUdot(constEigen::Ref<constkm::Vec>&u_dot);

voiddump(std::string_viewprefix="")const;

boolcheckVelJacobian(doubletolerance=1e-6)const;

protected:
CoordIndices
_updateFrozenCoords(CoordIndices&frozen_coords,size_tnew_frozen,size_ttotal_n);

CoordIndices
_updateNonFrozenCoords(CoordIndices&frozen_coords,size_tnew_unfrozen,size_ttotal_n);

kc::ks_ptr<km::NonlinearSolver>_nlSolver(){return_nl_solver;}

void_discard(kc::ks_ptr<Base>&base)override;

void_makeHealthy()override;

void_setup();

void_makeNlSolver();

void_updateDependentCoordsQ();

km::Mat_velJacobianNumDiff()const;

CoordIndices_recomputeNonFrozenCoord(constCoordIndices&frozen_coords,size_tn)const;

protected:
kc::ks_ptr<SubGraph>_sg;

kc::ks_ptr<MultiJacobianGenerator>_cutjoint_loops_multi_jacgen=nullptr;

kc::ks_ptr<MultiJacobianGenerator>_non_cutjoint_loops_multi_jacgen=nullptr;

kc::ks_ptr<km::NonlinearSolver>_nl_solver=nullptr;

CoordIndices_dependent_Q_coord_indices;

CoordIndices_frozen_indices_Q;

CoordIndices_frozen_indices_U;

CoordIndices_nonfrozen_indices_Q;

CoordIndices_nonfrozen_indices_U;
};

}//namespaceKarana::Dynamics