Search  

---

C:/temp/src/j2k/Test/VectorTest.cpp

Go to the documentation of this file.

///////////////////////////////////////////////////////////////////////////
//
// vectst.cpp - Test vector functions
//   Tests the following functions:
//    Add vectors c = a + b
//    Subtract vectors c = a - b
//    Multiply vectors c = a * b
//    Dot product c = sum( a * b )
//    Sum of vector s = sum( a )
//    Scale and copy a vector
//    Comparison ( a == b ) and ( a != b )
//
//   Inputs:
//    None
//
//   Outputs:
//    Print of test vectors and results
//
///////////////////////////////////////////////////////////////////////////
#include <j2k/Fred/Math/Vector.hpp>

#error "Not Working Yet! Test to be modified"

int main()
{
  // Initialize vector with data
  static float x[] = { 9.0f, 8.0f, 7.0f, 6.0f, 5.0f };
  static float y[] = { 1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f };
  Vector<float> v1( ELEMENTS(x), x );
  Vector<float> v2( ELEMENTS(y), y );
  Vector<float> v3;

  // Sum of two float vectors
  v3 = add( v1, v2 );
  cout << "Vector sum: v3 = v1 + v2\n";
  cout << "v1   v2   v3\n";
  for( int i = 0 ; i < v3.length(); i++ )
    cout << v1[i] << " " << v2[i] << " " << v3[i] << endl;

  // Difference of two float vectors
  v3 = sub( v1, v2 );
  cout << "\nVector subtract: v3 = v1 - v2\n";
  cout << "v1   v2   v3\n";
  for( i = 0 ; i < v3.length(); i++ )
  cout << v1[i] << " " << v2[i] << " " << v3[i] << endl;

  // Product of two float vectors
  v3 = pwisemult( v1, v2 );
  cout << "\nVector product: v3 = v1 * v2\n";
  cout << "v1   v2   v3\n";
  for( i = 0 ; i < v3.length(); i++ )
  cout << v1[i] << " " << v2[i] << " " << v3[i] << endl;

  // New vector v4 from subrange of vector v2
  Vector<float> v4 = v2( 2, 5 );

  // Offset product of two float vectors
  v3 = pwisemult( v1, v4 );
  cout << "\nVector product: v3 = v1 * v4\n";
  cout << "v1   v4   v3\n";
  for( i = 0 ; i < v3.length(); i++ )
  cout << v1[i] << " " << v4[i] << " " << v3[i] << endl;

  // Dot product of two float vectors
  float fDot = mult( v1, v2 );
  cout << "\nDot product of v1 and v2 = " << fDot << endl;

  // Sum of all elements of a float vector
  float fSum = sum( v1 );
  cout << "\nSum of v1 = " << fSum << endl;

  // Format floating point numbers
  cout.precision( 1 );
  cout.setf( ios::fixed, ios::floatfield );

  // Convert a vector to a different type
  Vector<int> vInt( v1.length() );
  conv( vInt, v1 );
  cout << "\nConvert float to int: vInt = v1\n";
  cout << "v1  vInt\n";
  for( i = 0 ; i < vInt.length(); i++ )
  cout << v1[i] << "  " << vInt[i] << endl;

  // Scale and convert a vector to a different type
  conv( vInt, scale( 0.5, v1 ) );
  cout << "\nScale and convert float to int: vInt = 0.5 * v1\n";
  cout << "v1  vInt\n";
  for( i = 0 ; i < vInt.length(); i++ )
  cout << v1[i] << "  " << vInt[i] << endl;

  // Set elements of a vector to a constant
  cout << "\nAll elements of a vector set to 5.0\n";
  v2 = 5.0f;
  cout << v2;

  // Compare two vectors
  v3 = v2;
  cout << v3;
  cout << ( ( v2 == v3 ) ? "Vectors equal\n" : "Vectors not equal\n" );

  // Change vector slightly
  v3[2] = -v3[2];
  cout << v3;
  cout << ( ( v2 == v3 ) ? "Vectors equal\n" : "Vectors not equal\n" );

  return 0;
}

---