[pre-commit.ci lite] apply automatic fixes

Author: pre-commit-ci-lite[bot]

arch/arch_device_cuda.h

@@ -539,7 +539,7 @@ namespace arch{
                printf("The blocksize (%u) does not match with any of the predetermined block sizes! The error occurred in %s at line %d\n", blocksize, __FILE__, __LINE__);
                exit(1);
             }
-         }      
+         }
          /* Check for kernel launch errors */
          CHK_ERR(cudaPeekAtLastError());
          /* Synchronize and free the thread data allocation */

arch/gpu_base.hpp

@@ -281,7 +281,7 @@ struct GPUMemoryManager {
          std::cerr << "Error: Pointer name '" << name << "' not found.\n";
          return false;
       }
-      
+
       if (gpuMemoryPointers[name] != nullptr) {
          CHK_ERR( gpuFree(gpuMemoryPointers[name]) );
       }

backgroundfield/backgroundfield.cpp

@@ -124,7 +124,7 @@ void setBackgroundField(
                for(uint fComponent=0;fComponent<3;fComponent++){
                   T3DFunction valueFunction = std::bind(bgFunction, std::placeholders::_1, std::placeholders::_2, std::placeholders::_3, (coordinate)fComponent, 0, (coordinate)0);
                   BgBGrid.get(x,y,z)->at(fsgrids::bgbfield::BGBXVOL+fComponent) += volumeAverage(valueFunction,accuracy,start.data(),end);
-                  
+
                   //Compute derivatives. Note that we scale by dx[] as the arrays are assumed to contain differences, not true derivatives!
                   for(uint dComponent=0;dComponent<3;dComponent++){
                      T3DFunction derivFunction = std::bind(bgFunction, std::placeholders::_1, std::placeholders::_2, std::placeholders::_3, (coordinate)fComponent, 1, (coordinate)dComponent);
@@ -158,4 +158,3 @@ void setBackgroundFieldToZero(
       }
    }
 }
-

backgroundfield/backgroundfield.h

@@ -69,7 +69,7 @@ template<long unsigned int numFields> void setPerturbedFieldToZero(
     This function is used along with a field function to append the given function values to the
     perturbed B grid.
     A special other use case is, if using a Magnetosphere with dipole type 4 (vector dipole), when
-    the corrective terms to vanish the dipole towards the inflow boundary are stored in 
+    the corrective terms to vanish the dipole towards the inflow boundary are stored in
     the backgroundfield FSgrid object at offset fsgrids::bgbfield::BGBXVDCORR
 */
 template<long unsigned int numFields> void setPerturbedField(
@@ -125,12 +125,11 @@ template<long unsigned int numFields> void setPerturbedField(
                                  dx[faceCoord2[fComponent]]
                                 );
 
-	    }
-	    // Derivatives or volume averages are not calculated for the perBField
-	 }
+            }
+            // Derivatives or volume averages are not calculated for the perBField
+         }
       }
    }
 }
 
 #endif
-

backgroundfield/constantfield.cpp

@@ -47,9 +47,3 @@ double ConstantField::operator()( double x, double y, double z, coordinate compo
       return 0.0;
    }
 }
-
-
-
-
-
-

backgroundfield/constantfield.hpp

@@ -34,7 +34,7 @@ class ConstantField {
    bool _initialized;
    double _B[3]; // constant backgroundfield
 public:
-  
+
    ConstantField(){
       this->_initialized = false;
    }
@@ -44,4 +44,3 @@ class ConstantField {
 };
 
 #endif
-

backgroundfield/dipole.cpp

@@ -42,40 +42,40 @@ void Dipole::initialize(const double moment,const double center_x, const double
 
 
 double Dipole::operator()( double x, double y, double z, coordinate component, unsigned int derivative, coordinate dcomponent) const {
-   const double minimumR=1e-3*physicalconstants::R_E; //The dipole field is defined to be outside of Earth, and units are in meters     
+   const double minimumR=1e-3*physicalconstants::R_E; //The dipole field is defined to be outside of Earth, and units are in meters
    if(this->initialized==false) {
       return 0.0;
    }
    double r[3];
-   
+
    r[0]= x-center[0];
    r[1]= y-center[1];
    r[2]= z-center[2];
-   
+
    double r2 = r[0]*r[0]+r[1]*r[1]+r[2]*r[2];
-   
+
    if(r2<minimumR*minimumR) {
       //  r2=minimumR*minimumR;
       return 0.0; //set zero field inside dipole
    }
-   
+
    const double r5 = (r2*r2*sqrt(r2));
    const double rdotq=q[0]*r[0] + q[1]*r[1] +q[2]*r[2];
-   
+
    const double B=( 3*r[component]*rdotq-q[component]*r2)/r5;
-   
+
    if(derivative == 0) {
       //Value of B
       return B;
    } else if(derivative == 1) {
-      //first derivatives       
+      //first derivatives
       unsigned int sameComponent;
       if(dcomponent==component) {
          sameComponent=1;
       } else {
          sameComponent=0;
       }
-      
+
       return -5*B*r[dcomponent]/r2+
          (3*q[dcomponent]*r[component] -
           2*q[component]*r[dcomponent] +
@@ -84,9 +84,3 @@ double Dipole::operator()( double x, double y, double z, coordinate component, u
 
    return 0; // dummy, but prevents gcc from yelling
 }
-
-
-
-
-
-

backgroundfield/dipole.hpp

@@ -35,12 +35,11 @@ class Dipole {
    double q[3];      // Dipole moment; set to (0,0,moment)
    double center[3]; // Coordinates where the dipole sits; set to (0,0,0)
 public:
-   
+
    Dipole() {}
 
    void initialize(const double moment,const double center_x, const double center_y, const double center_z, const double tilt_angle);
    double operator()(double x, double y, double z, coordinate component, unsigned int derivative=0, coordinate dcomponent=X) const;
 };
 
 #endif
-

backgroundfield/fieldfunction.hpp

@@ -32,4 +32,3 @@ Background magnetic field class of Vlasiator.
 
 typedef std::function<double(double x, double y, double z, coordinate component, unsigned int derivative, coordinate dcomponent)> FieldFunction;
 #endif
-

backgroundfield/integratefunction.cpp

@@ -47,7 +47,7 @@ double lineAverage(
    const double acc = accuracy*L;
    const double a = r1[line];
    const double b = r1[line] + L;
-   
+
    switch (line) {
       case X:
       {
@@ -77,7 +77,7 @@ double lineAverage(
 
 
 double surfaceAverage(
-   const T3DFunction& f1, 
+   const T3DFunction& f1,
    coordinate face, double accuracy,
    const double r1[3],
    double L1,
@@ -97,7 +97,7 @@ double surfaceAverage(
       case Y:
       {
          T2DFunction f = std::bind(f1,std::placeholders::_1,r1[1],std::placeholders::_2);
-         value = Romberg(f, r1[0],r1[0]+L1, r1[2],r1[2]+L2, acc)*norm; 
+         value = Romberg(f, r1[0],r1[0]+L1, r1[2],r1[2]+L2, acc)*norm;
       }
       break;
       case Z:
@@ -127,4 +127,3 @@ double volumeAverage(
    value= Romberg(f1, r1[0],r2[0], r1[1],r2[1], r1[2],r2[2], acc)*norm;
    return value;
 }
-