example.cc


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

// -*-C++-*-

#include "vecmathlib.h"

#include <iostream>

using namespace std;
using namespace vecmathlib;

int main(int argc, char **argv) {
  // Declare a double precision vector with an architecture-dependent
  // number of elements
  float64_vec x;
  // Set each element separately. This is inefficient and should be
  // avoided if possible, but we want to demonstrate it here anyway.
  for (int i = 0; i < float64_vec::size; ++i)
    x.set_elt(i, double(i));
  float64_vec y = x + float64_vec(1.0);
  y = sqrt(y);
  float64_vec z = log(y);

  // Boolean vectors are closely related to either double or float
  // vectors, thus we need to make a distinction
  bool64_vec b = x < y;
  // Integer vectors are closely related to either double or float,
  // thus we need to make a distinction -- there is "long_vec"
  // corresponding to "float64_vec", and there is "int_vec"
  // correpsonding to "float_vec".
  int64_vec i = convert_int(y);

  cout << "x=" << x << "\n";
  cout << "y=" << y << "\n";
  cout << "z=" << z << "\n";
  cout << "b=" << b << "\n";
  cout << "i=" << i << "\n";

  return 0;
}