ASubtraction.cpp

Go to the documentation of this file.

/**
 * @file ASubtraction.cpp
 * Surface Modeling Library
 * CS497
 *
 * @author William Nagel
 */

#include "ASubtraction.h"

double ASubtraction::h(double f, double g)
{
  if ((f <= m_r1) && (g <= m_r2))
    return (((f - m_r1) * (f - m_r1)) / (m_r1 * m_r1)) +
           (((g - m_r2) * (g - m_r2)) / (m_r2 * m_r2)) - 1;
  else 
    return 0.0;
}

Intervald ASubtraction::h(Intervald f, Intervald g)
{
  if ((f <= m_r1) && (g <= m_r2))
    {
      Intervald r1(m_r1);
      Intervald r2(m_r2);
      return ((f - r1).squared() / r1.squared()) +
             ((g - r2).squared() / r2.squared()) - 
             Intervald(1.0);
    }
  else
    return Intervald(0.0);
}

double ASubtraction::hf(double f, double g)
{
  if ((f <= m_r1) && (g <= m_r2))
    return ((2.0 * f - 2.0 * m_r1) / (m_r1 * m_r1));
  else
    return 0.0;
}

Intervald ASubtraction::hf(Intervald f, Intervald g)
{
  if((f <= m_r1) && (g <= m_r2))
    {
      Intervald r1(m_r1);
      Intervald two(2.0);
      return ((two * f) - (two * r1)) / r1.squared();
    }
  else
    return Intervald(0.0);
}

double ASubtraction::hg(double f, double g)
{
  if ((f <= m_r1) && (g <= m_r2))
    return ((2.0 * g - 2.0 * m_r2) / (m_r2 * m_r2));
  else
    return 0.0;
}

Intervald ASubtraction::hg(Intervald f, Intervald g)
{
  if ((f <= m_r1) && (g <= m_r2))
    {
      Intervald r2(m_r2);
      Intervald two(2.0);
      return ((two * g) - (two * r2)) / r2.squared();
    }
  else
    return Intervald(0.0);
}

double ASubtraction::hff(double f, double g)
{
  if ((f <= m_r1) && (g <= m_r2))
    return (2.0 / (m_r1 * m_r1));
  else
    return 0.0;
}

Intervald ASubtraction::hff(Intervald f, Intervald g)
{
  if ((f <= m_r1) && (g <= m_r2))
    {
      Intervald r1(m_r1);
      return Intervald(2.0) / r1.squared();
    }
  else
    return Intervald(0.0);
}

double ASubtraction::hgg(double f, double g)
{
  if ((f <= m_r1) && (g <= m_r2))
    return (2.0 / (m_r2 * m_r2));
  else 
    return 0.0;
}

Intervald ASubtraction::hgg(Intervald f, Intervald g)
{
  if ((f <= m_r1) && (g <= m_r2))
    {
      Intervald r2(m_r2);
      return Intervald(2.0) / r2.squared();
    }
  else
    return 0.0;
}

double ASubtraction::hfg(double f, double g)
{
  return 0.0;
}

Intervald ASubtraction::hfg(Intervald f, Intervald g)
{
  return Intervald(0.0);
}

---