solverInconsistent< T > Class Template Reference

Solve 1D and 2D linear equations with one move variable than the number of equations. More...

#include <mathlib.h>

List of all members.

Static Public Member Functions
static boolc	d1linearequ (T &x, T &y, T const a0, T const a1, T const c)
	Find a solution for the equation a0*x+a1*y=c.
static boolc	d2linearequ (T &x, T &y, T &z, T const a00, T const a01, T const a02, T const c0, T const a10, T const a11, T const a12, T const c1)
	Find a solution for the equations { a00*x+a01*y=c0, a10*x+a11*y=c1 }.

---

Detailed Description

template<typename T>
class solverInconsistent< T >

Solve 1D and 2D linear equations with one move variable than the number of equations.

Find a point that satisfies the equations.

For example plane to plane intersection gives two equations in three unknowns.

Definition at line 161 of file mathlib.h.

---

Member Function Documentation

template<typename T >

boolc solverInconsistent< T >::d1linearequ	(	T &	x,
		T &	y,
		T const	a0,
		T const	a1,
		T const	c
	)			[inline, static]

Find a solution for the equation a0*x+a1*y=c.

Definition at line 1021 of file mathlib.h.

{
  if (zero<T>::test(a0))
  {
    x = 0;
    if (zero<T>::test(a1))
    {
      y = 0;
      return (c==0);
    }
    y = c/a1;
    return true;
  }

  if (zero<T>::test(a1))
  {
    y = 0;
    x = c/a0;
    return true;
  }

  x=1.0;
  y = (c-a0)/a1;
  return true;
}

template<typename T >

boolc solverInconsistent< T >::d2linearequ	(	T &	x,
		T &	y,
		T &	z,
		T const	a00,
		T const	a01,
		T const	a02,
		T const	c0,
		T const	a10,
		T const	a11,
		T const	a12,
		T const	c1
	)			[inline, static]

Find a solution for the equations { a00*x+a01*y=c0, a10*x+a11*y=c1 }.

Definition at line 1056 of file mathlib.h.

References solver< T >::d2linearequ().

Referenced by plane::intersects().

{
  if (zero<T>::test(a00))
  {
    if (zero<T>::test(a10))
    {
      x=0;
      return solver<T>::d2linearequ(y,z,a01,a02,a10,a11,c0,c1);
    }

    return d2linearequ(x,y,z,a10,a11,a12,c1,a00,a01,a02,c0);
  }

  T e[3];
  e[0] = a11-a10*a01/a00; 
  e[1] = a12-a10*a02/a00;
  e[2] = c1-a10*c0/a00;

  bool res = d1linearequ(y,z,e[0],e[1],e[2]);
  if (res==false)
    return false;

  x = (c0-a01*y-a02*z)/a00;

  return true;
}

---

The documentation for this class was generated from the following file:

• mathlib/mathlib.h