AIntersection.cpp

Go to the documentation of this file.

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

#include "AIntersection.h"

AIntersection::AIntersection(Implicit *f, Implicit *g, int cont)
{
  m_f = f; 
  m_g = g;
  m_cont = cont;
  m_r1 = .3200;
  m_r2 = .3200;
}

double AIntersection::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 RFunction::h(f,g);
}

Intervald AIntersection::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 RFunction::h(f,g);
}

double AIntersection::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 RFunction::hf(f,g);
}

Intervald AIntersection::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 RFunction::hf(f,g);
}

double AIntersection::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 RFunction::hg(f,g);
}

Intervald AIntersection::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 RFunction::hg(f,g);
}

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

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

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

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

double AIntersection::hfg(double f, double g)
{
  if ((f <= m_r1) && (g <= m_r2))
    return 0.0;
  else
    return RFunction::hfg(f,g);
}

Intervald AIntersection::hfg(Intervald f, Intervald g)
{
  if ((f <= m_r1) && (g <= m_r2))
    return Intervald(0.0);
  else
    return RFunction::hfg(f,g);
}

---