Skip to content

Commit

Permalink
Remove standard output messages (#129)
Browse files Browse the repository at this point in the history
  • Loading branch information
alxbilger authored Jan 15, 2024
1 parent fa3a036 commit 88ec684
Show file tree
Hide file tree
Showing 5 changed files with 0 additions and 43 deletions.
1 change: 0 additions & 1 deletion src/BeamAdapter/component/BeamInterpolation.inl
Original file line number Diff line number Diff line change
Expand Up @@ -788,7 +788,6 @@ void BeamInterpolation<DataTypes>::getSplinePoints(unsigned int edgeInList, cons
return;
}

//std::cout << " getSplinePoints : global_H_local0 ="<<global_H_local0<<" global_H_local1 ="<<global_H_local1<<std::endl;
const Real& _L = d_lengthList.getValue()[edgeInList];
this->getControlPointsFromFrame(global_H_local0, global_H_local1,_L,P0, P1,P2, P3);

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -177,7 +177,6 @@ void InterventionalRadiologyController<DataTypes>::loadMotionData(std::string fi
{
if (line.empty())
continue;
std::cout << line << std::endl;
std::istringstream values(line);
values >> result[0] >> result[1] >> result[2];
result[0] /= 1000;
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -671,8 +671,6 @@ void AdaptiveInflatableBeamForceField<DataTypes>::addForce (const MechanicalPara
}
}

//std::cout<<"F beam = "<<f<<std::endl;

if(d_computeMass.getValue())
{
/// add gravity:
Expand Down
5 changes: 0 additions & 5 deletions src/BeamAdapter/component/mapping/BeamLengthMapping.h
Original file line number Diff line number Diff line change
Expand Up @@ -166,11 +166,6 @@ class BeamLengthMapping : public Mapping<TIn, TOut>
void setInterpolation(BeamInterpolation< TIn >* interpolation)
{
l_adaptativebeamInterpolation.set(interpolation);

std::cout<<"l_adaptativebeamInterpolation set "<<l_adaptativebeamInterpolation.getName()<<std::endl;
std::cout<<"l_adaptativebeamInterpolation num Beam "<<l_adaptativebeamInterpolation->getNumBeams()<<std::endl;


}

virtual ~BeamLengthMapping(){}
Expand Down
34 changes: 0 additions & 34 deletions src/BeamAdapter/component/mapping/BeamLengthMapping.inl
Original file line number Diff line number Diff line change
Expand Up @@ -228,11 +228,6 @@ void BeamLengthMapping< TIn, TOut>::applyJ(const core::MechanicalParams* mparams
template <class TIn, class TOut>
void BeamLengthMapping< TIn, TOut>::applyJT(const core::MechanicalParams* mparams, Data<InVecDeriv>& dOut, const Data<VecDeriv>& dIn)
{

if(this->f_printLog.getValue() ){
std::cout<<" --------------------- applyJT --------------------"<<std::endl;
}

SOFA_UNUSED(mparams);

sofa::helper::ScopedAdvancedTimer timer("AdaptiveBeamMappingMechanicalApplyJT");
Expand Down Expand Up @@ -300,12 +295,6 @@ void BeamLengthMapping< TIn, TOut>::applyJT(const core::MechanicalParams* mparam
Transform DOF0Global_H_local0, DOF1Global_H_local1;
l_adaptativebeamInterpolation->getDOFtoLocalTransformInGlobalFrame(i, DOF0Global_H_local0, DOF1Global_H_local1, x_in);


if(this->f_printLog.getValue() ){
std::cout<<"DOF0Global_H_local0 = "<<DOF0Global_H_local0<<std::endl;
std::cout<<"DOF1Global_H_local1 = "<<DOF1Global_H_local1<<std::endl;
}

// rotate back the force to the local frame
SpatialVector f0, f1,f2,f3;
f0.setForce( DOF0Global_H_local0.getOrientation().inverseRotate(F0) );
Expand Down Expand Up @@ -444,9 +433,6 @@ void BeamLengthMapping< TIn, TOut>::applyJT(const core::ConstraintParams* cparam
template <class TIn, class TOut>
void BeamLengthMapping< TIn, TOut>::applyDJT(const MechanicalParams* mparams, core::MultiVecDerivId parentDfId, core::ConstMultiVecDerivId childDfId)
{
if(this->f_printLog.getValue() ){
std::cout<<" --------------------- applyDJT --------------------"<<std::endl;
}
const unsigned& geometricStiffness = d_geometricStiffness.getValue();
if( !geometricStiffness ) return;

Expand Down Expand Up @@ -607,10 +593,6 @@ void BeamLengthMapping< TIn, TOut>::applyDJT(const MechanicalParams* mparams, co
template <class TIn, class TOut>
void BeamLengthMapping<TIn, TOut>::updateK(const core::MechanicalParams* mparams, core::ConstMultiVecDerivId childForceId )
{

if(this->f_printLog.getValue() ){
std::cout<<" entering updateK ..........."<<std::endl;
}
const unsigned& geometricStiffness = d_geometricStiffness.getValue();
if( !geometricStiffness ) { K_geom.resize(0,0); return; }
//helper::ReadAccessor<Data<VecDeriv> > childForce( *childForceId[(const core::State<TOut>*)this->getToModels()[0]].read() );
Expand All @@ -633,11 +615,6 @@ void BeamLengthMapping<TIn, TOut>::updateK(const core::MechanicalParams* mparams
// if stabilized GS (geometricStiffness==2) -> keep only force in extension
if( childForce[i][0] < 0 || geometricStiffness==1 )
{

if(this->f_printLog.getValue() ){
std::cout<<" step1 ........... childForce[i][0] = "<<childForce[i][0]<<std::endl;
}

//1. get the indices of the Dofs of the beam a
unsigned int IdxNode[2];;
l_adaptativebeamInterpolation->getNodeIndices(i,IdxNode[0],IdxNode[1]);
Expand Down Expand Up @@ -932,17 +909,6 @@ void BeamLengthMapping<TIn, TOut>::computeDJtSpline(const Real &f_input, const V
// 1 compute the derivatives of these relations:
//F_IP[0] = IP[0]*f_input*(18.0 + B)/(72.0*IP[0].norm());

if(this->f_printLog.getValue() ){
std::cout<<"computeDJtSpline :"<<std::endl;
for(unsigned int k=0; k<4; k++){
for(unsigned int l=0; l<4;l++){
std::cout<<" "<<a[k][l];
}
std::cout<<" "<<std::endl;
}
}


Mat3 K_IP[4];

for(unsigned int k=0; k<4; k++){
Expand Down

0 comments on commit 88ec684

Please sign in to comment.