helper.h

Go to the documentation of this file.

#ifndef MOJA_MODULES_CBM_HELPER_H_
#define MOJA_MODULES_CBM_HELPER_H_
 
#include <boost/algorithm/string.hpp> 
 
#include "moja/modules/cbm/_modules.cbm_exports.h"
 
namespace moja {
    namespace modules {
        namespace cbm {
 
            class CBM_API Helper {
            public:
                Helper() = delete;
                ~Helper() = delete;
 
                static double calculateMerchFactor(double volume, double a, double b);
 
                static double calculateNonMerchFactor(
                    double merchStemwood, double a_nonmerch,
                    double b_nonmerch, double k_nonmerch);
 
                static double calculateSaplingFactor(double stemwood, double k_sap,
                    double a_sap, double b_sap);
 
                static double modelTerm(double vol, double a1, double a2, double a3);
 
                static bool runMoss(int gcID, std::string mossLeadingSpecies, std::string speciesName);
            };
 
            inline double Helper::calculateMerchFactor(double volume, double a, double b) {
                return a * pow(volume, b);
            }
 
            inline double Helper::calculateNonMerchFactor(
                double merchStemwood, double a_nonmerch, double b_nonmerch, double k_nonmerch) {
 
                if (merchStemwood > 0.0) {
                    return k_nonmerch + a_nonmerch * pow(merchStemwood, b_nonmerch);
                }
 
                return exp(a_nonmerch);
            }
 
            inline double Helper::calculateSaplingFactor(
                double stemwood, double k_sap, double a_sap, double b_sap) {
 
                if (stemwood > 0.0) {
                    return k_sap + a_sap * pow(stemwood, b_sap);
                }
 
                return exp(a_sap);
            }
 
            inline double Helper::modelTerm(double vol, double a1, double a2, double a3) {
                return exp(a1 + a2 * vol + a3 * log(vol + 5.0));
            }
 
            inline bool Helper::runMoss(int gcID, std::string mossLeadingSpecies, std::string speciesName) {
                boost::algorithm::to_lower(mossLeadingSpecies);
                boost::algorithm::to_lower(speciesName);
 
                return (gcID > 0 && boost::contains(speciesName, mossLeadingSpecies));
            }
        }
    }
}
#endif // MOJA_MODULES_CBM_HELPER_H_