IASqrt.cpp

#include "GPSymbRegEvalOp.h"
 
#include "IASqrt.h"
#include "IATree.h"
#include <cmath>
 
void IASqrt::execute(void* result, Tree::Tree& tree)
{
    if (!calculateBoundsIfRootNode(result, tree)) {
        return;
    }
    
    Tree::Primitives::Sqrt::execute(result, tree);
}
 
IABounds IASqrt::calculateBounds(Tree::IATree& tree) {
    // get next argument
    tree.iNode_++;
    IAPrimitiveP prim = boost::dynamic_pointer_cast<IAPrimitive>(tree[tree.iNode_]->primitive_);
    // calculate bounds
    IABounds bounds = prim->calculateBounds(tree);
 
    if (!bounds.isWithinBounds) {
        // subtree not within bounds, delete expression
        return bounds;
    } else if (bounds.lowerBound < 0) {
        IABounds retVal;
        retVal.isWithinBounds = false;
        return retVal;
    } else {
        // subtree is ok
        IABounds retVal;
        retVal.isWithinBounds = true;
        retVal.isLowerBoundInclusive = bounds.isLowerBoundInclusive;
        retVal.lowerBound = sqrt(bounds.lowerBound);
        retVal.isUpperBoundInclusive = bounds.isUpperBoundInclusive;
        retVal.upperBound = sqrt(bounds.upperBound);
        return retVal;
    }
}