force.h

Go to the documentation of this file.

#pragma once
 
#include "Common/Float.h"
 
namespace AR {
 
 
    struct Force {
    public:
        Float acc_in[3];   
        Float acc_pert[3]; 
        Float pot_in;      
        Float pot_pert;    
#ifdef AR_TTL
        Float gtgrad[3];   
 
 
        Force(): acc_in{Float(0.0),Float(0.0),Float(0.0)}, 
                 acc_pert{Float(0.0),Float(0.0),Float(0.0)},
                 pot_in(0.0), pot_pert(0.0),
                 gtgrad{Float(0.0),Float(0.0),Float(0.0)} {}
#else
 
        Force(): acc_in{Float(0.0),Float(0.0),Float(0.0)}, 
                 acc_pert{Float(0.0),Float(0.0),Float(0.0)},
                 pot_in(0.0), pot_pert(0.0) {}
#endif
 
 
        void writeBinary(FILE *_fout) {
            fwrite(this, sizeof(*this),1,_fout);
        }
 
 
        void readBinary(FILE *_fin) {
            size_t rcount = fread(this, sizeof(*this), 1, _fin);
            if (rcount<1) {
                std::cerr<<"Error: Data reading fails! requiring data number is 1, only obtain "<<rcount<<".\n";
                abort();
            }
        }
 
        void clear() {
            acc_in[0] = acc_in[1] = acc_in[2] = 0.0;
            acc_pert[0] = acc_pert[1] = acc_pert[2] = 0.0;
            pot_in = pot_pert = 0.0;
#ifdef AR_TTL
            gtgrad[0] = gtgrad[1] = gtgrad[2] = 0.0;
#endif            
        }
    };
 
 
}