CARLsim  4.1.0
CARLsim: a GPU-accelerated SNN simulator
linear_algebra.h
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2014 Regents of the University of California. All rights reserved.
3  *
4  * Redistribution and use in source and binary forms, with or without
5  * modification, are permitted provided that the following conditions
6  * are met:
7  *
8  * 1. Redistributions of source code must retain the above copyright
9  * notice, this list of conditions and the following disclaimer.
10  *
11  * 2. Redistributions in binary form must reproduce the above copyright
12  * notice, this list of conditions and the following disclaimer in the
13  * documentation and/or other materials provided with the distribution.
14  *
15  * 3. The names of its contributors may not be used to endorse or promote
16  * products derived from this software without specific prior written
17  * permission.
18  *
19  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
20  * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
21  * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
22  * A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
23  * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
24  * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
25  * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
26  * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
27  * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
28  * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
29  * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
30  *
31  * *********************************************************************************************** *
32  * CARLsim
33  * created by: (MDR) Micah Richert, (JN) Jayram M. Nageswaran
34  * maintained by: (MA) Mike Avery <averym@uci.edu>, (MB) Michael Beyeler <mbeyeler@uci.edu>,
35  * (KDC) Kristofor Carlson <kdcarlso@uci.edu>
36  * (TSC) Ting-Shuo Chou <tingshuc@uci.edu>
37  *
38  * CARLsim available from http://socsci.uci.edu/~jkrichma/CARLsim/
39  * Ver 10/11/2014
40  */
41 
42 #ifndef _LINEAR_ALGEBRA_H_
43 #define _LINEAR_ALGEBRA_H_
44 
45 #include <ostream> // print struct info
46 
47 unsigned long int get_time_ms64();
48 
57 struct Point3D {
58 public:
59  Point3D(int _x, int _y, int _z) : x(1.0*_x), y(1.0*_y), z(1.0*_z) {}
60  Point3D(double _x, double _y, double _z) : x(_x), y(_y), z(_z) {}
61 
62  // print struct info
63  friend std::ostream& operator<<(std::ostream &strm, const Point3D &p) {
64  strm.precision(2);
65  return strm << "Point3D=(" << p.x << "," << p.y << "," << p.z << ")";
66  }
67 
68  // overload operators
69  Point3D operator+(const double a) const { return Point3D(x+a,y+a,z+a); }
70  Point3D operator+(const Point3D& p) const { return Point3D(x+p.x,y+p.y,z+p.z); }
71  Point3D operator-(const double a) const { return Point3D(x-a,y-a,z-a); }
72  Point3D operator-(const Point3D& p) const { return Point3D(x-p.x,y-p.y,z-p.z); }
73  Point3D operator*(const double a) const { return Point3D(x*a,y*a,z*a); }
74  Point3D operator*(const Point3D& p) const { return Point3D(x*p.x,y*p.y,z*p.z); }
75  Point3D operator/(const double a) const { return Point3D(x/a,y/a,z/a); }
76  Point3D operator/(const Point3D& p) const { return Point3D(x/p.x,y/p.y,z/p.z); }
77  bool operator==(const Point3D& p) const { return Equals(p); }
78  bool operator!=(const Point3D& p) const { return !Equals(p); }
79  bool operator<(const Point3D& p) const { return (CompareTo(p)<0); }
80  bool operator>(const Point3D& p) const { return (CompareTo(p)>0); }
81  bool operator<=(const Point3D& p) const { return (CompareTo(p)<=0); }
82  bool operator>=(const Point3D& p) const { return (CompareTo(p)>=0); }
83 
84  // coordinates
85  double x, y, z;
86 
87 private:
88  bool Equals(const Point3D& p) const { return (x==p.x && y==p.y && z==p.z); }
89  int CompareTo(const Point3D& p) const { return (x>p.x&&y>p.y&&z>p.z) ? 1 : ( (x<p.x&&y<p.y&&z<p.z) ? -1 : 0); }
90 };
91 
92 double dist(Point3D& p1, Point3D& p2);
93 
95 double norm2(const Point3D& p);
96 
98 double norm(const Point3D& p);
99 
101 //bool isPointOnGrid(Point3D& p, Grid3D& g);
102 
103 #endif
Point3D::operator<<
friend std::ostream & operator<<(std::ostream &strm, const Point3D &p)
Definition: linear_algebra.h:63
dist
double dist(Point3D &p1, Point3D &p2)
Definition: linear_algebra.cpp:5
Point3D::operator==
bool operator==(const Point3D &p) const
Definition: linear_algebra.h:77
Point3D::operator+
Point3D operator+(const Point3D &p) const
Definition: linear_algebra.h:70
Point3D::operator<
bool operator<(const Point3D &p) const
Definition: linear_algebra.h:79
Point3D::operator>=
bool operator>=(const Point3D &p) const
Definition: linear_algebra.h:82
Point3D::operator*
Point3D operator*(const double a) const
Definition: linear_algebra.h:73
Point3D::operator+
Point3D operator+(const double a) const
Definition: linear_algebra.h:69
Point3D::operator-
Point3D operator-(const Point3D &p) const
Definition: linear_algebra.h:72
Point3D::operator-
Point3D operator-(const double a) const
Definition: linear_algebra.h:71
Point3D::operator/
Point3D operator/(const Point3D &p) const
Definition: linear_algebra.h:76
Point3D::operator!=
bool operator!=(const Point3D &p) const
Definition: linear_algebra.h:78
Point3D::operator*
Point3D operator*(const Point3D &p) const
Definition: linear_algebra.h:74
get_time_ms64
unsigned long int get_time_ms64()
Point3D::operator<=
bool operator<=(const Point3D &p) const
Definition: linear_algebra.h:81
Point3D::operator/
Point3D operator/(const double a) const
Definition: linear_algebra.h:75
Point3D::operator>
bool operator>(const Point3D &p) const
Definition: linear_algebra.h:80
Point3D
a point in 3D space
Definition: linear_algebra.h:57
Point3D::y
double y
Definition: linear_algebra.h:85
Point3D::Point3D
Point3D(double _x, double _y, double _z)
Definition: linear_algebra.h:60
norm
double norm(const Point3D &p)
Definition: linear_algebra.cpp:14
norm2
double norm2(const Point3D &p)
calculate norm^2
Definition: linear_algebra.cpp:9
Point3D::Point3D
Point3D(int _x, int _y, int _z)
Definition: linear_algebra.h:59
Point3D::x
double x
Definition: linear_algebra.h:85
Point3D::z
double z
Definition: linear_algebra.h:85