From 88ec68443c72871bc8529a51b3177c0389ac0d2c Mon Sep 17 00:00:00 2001
From: Alex Bilger <alxbilger@users.noreply.github.com>
Date: Mon, 15 Jan 2024 13:41:46 +0100
Subject: [PATCH] Remove standard output messages (#129)

---
 .../component/BeamInterpolation.inl           |  1 -
 .../InterventionalRadiologyController.inl     |  1 -
 .../AdaptiveInflatableBeamForceField.inl      |  2 --
 .../component/mapping/BeamLengthMapping.h     |  5 ---
 .../component/mapping/BeamLengthMapping.inl   | 34 -------------------
 5 files changed, 43 deletions(-)

diff --git a/src/BeamAdapter/component/BeamInterpolation.inl b/src/BeamAdapter/component/BeamInterpolation.inl
index 44cf183cf..d5bd84fc7 100644
--- a/src/BeamAdapter/component/BeamInterpolation.inl
+++ b/src/BeamAdapter/component/BeamInterpolation.inl
@@ -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);
 
diff --git a/src/BeamAdapter/component/controller/InterventionalRadiologyController.inl b/src/BeamAdapter/component/controller/InterventionalRadiologyController.inl
index 1cc8a3410..5b96e4748 100644
--- a/src/BeamAdapter/component/controller/InterventionalRadiologyController.inl
+++ b/src/BeamAdapter/component/controller/InterventionalRadiologyController.inl
@@ -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;
diff --git a/src/BeamAdapter/component/forcefield/AdaptiveInflatableBeamForceField.inl b/src/BeamAdapter/component/forcefield/AdaptiveInflatableBeamForceField.inl
index 01c47b353..615e845f4 100644
--- a/src/BeamAdapter/component/forcefield/AdaptiveInflatableBeamForceField.inl
+++ b/src/BeamAdapter/component/forcefield/AdaptiveInflatableBeamForceField.inl
@@ -671,8 +671,6 @@ void AdaptiveInflatableBeamForceField<DataTypes>::addForce (const MechanicalPara
         }
     }
 
-    //std::cout<<"F beam = "<<f<<std::endl;
-
     if(d_computeMass.getValue())
     {
         /// add gravity:
diff --git a/src/BeamAdapter/component/mapping/BeamLengthMapping.h b/src/BeamAdapter/component/mapping/BeamLengthMapping.h
index e58e2bab5..46cb3eae5 100644
--- a/src/BeamAdapter/component/mapping/BeamLengthMapping.h
+++ b/src/BeamAdapter/component/mapping/BeamLengthMapping.h
@@ -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(){}
diff --git a/src/BeamAdapter/component/mapping/BeamLengthMapping.inl b/src/BeamAdapter/component/mapping/BeamLengthMapping.inl
index 35385090e..582bb0194 100644
--- a/src/BeamAdapter/component/mapping/BeamLengthMapping.inl
+++ b/src/BeamAdapter/component/mapping/BeamLengthMapping.inl
@@ -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");
@@ -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) );
@@ -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;
 
@@ -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() );
@@ -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]);
@@ -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++){