linalg::vector_view Class Reference

A vector that doesn't own its data; rather, points to data owned by another vector. More...

#include <linalg.hpp>

Inheritance diagram for linalg::vector_view:

[legend]

Collaboration diagram for linalg::vector_view:

[legend]

List of all members.

Public Member Functions
	~vector_view ()
vector_view &	operator= (const vector &w)
vector_view &	operator= (const vector_view &w)
	Assignment, useful for modifying chunks of vectors all at once.
Private Member Functions
	vector_view ()
	vector_view (gsl_vector *y, const slice &s)
	See above.
	vector_view (gsl_matrix *A, const slice &a, const size_t j)
	See above.
	vector_view (gsl_matrix *A, const size_t i, const slice &b)
	See above.
Friends
class	vector
class	matrix

---

Detailed Description

A vector that doesn't own its data; rather, points to data owned by another vector.

---

Constructor & Destructor Documentation

linalg::vector_view::~vector_view

(

)

                           {
    delete x;
    x = 0;
  }

linalg::vector_view::vector_view

(

)

[private]

vector_views are invisible to the user and can only be created using slices and matrix or vector indexing.

                          {
    x = 0;
  }

linalg::vector_view::vector_view	(	gsl_vector *	y,
		const slice &	s
	)			[private]

See above.

                                                       :vector(new gsl_vector){
    if(s.end() > y->size - 1){
      indexOutOfRange exc;
      exc.file = __FILE__;
      exc.line = __LINE__;
      exc.reason = "invalid vector slice upper range.";
      exc.i = s.end();
      exc.n = y->size;
      throw exc;
    }
    x -> size = (s.end() - s.begin())/s.stride()+1;
    x -> data = gsl_vector_ptr(y,s.begin());
    x -> stride = s.stride();
    x -> block = 0;
    x -> owner = 0;
  }

Here is the call graph for this function:

linalg::vector_view::vector_view	(	gsl_matrix *	A,
		const slice &	a,
		const size_t	j
	)			[private]

See above.

                                                                        :vector(new gsl_vector){
    if(a.end() > A -> size1 - 1 or j-1 > A -> size2-1){
      indexOutOfRange exc;
      exc.file = __FILE__;
      exc.line = __LINE__;
      exc.reason = "invalid matrix view range.";
      exc.i = a.end();
      exc.j = j;
      exc.n = A -> size1;
      exc.m = A -> size2;
      throw exc;
    }
    x -> size = (a.end() - a.begin())/a.stride()+1;
    x -> data = gsl_matrix_ptr(A, a.begin(), j-1);
    x -> stride = a.stride()*(A -> tda); //Note that a longer step is
                                         //necessary here.
    x -> block = 0;
    x -> owner = 0;
  }

Here is the call graph for this function:

linalg::vector_view::vector_view	(	gsl_matrix *	A,
		const size_t	i,
		const slice &	b
	)			[private]

See above.

                                                                        :vector(new gsl_vector){
    if(b.end() > A -> size2 - 1 or i-1 > A-> size1-1){
      indexOutOfRange exc;
      exc.file = __FILE__;
      exc.line = __LINE__;
      exc.reason = "invalid matrix view range.";
      exc.i = i;
      exc.j = b.end();
      exc.n = A -> size1;
      exc.m = A -> size2;
      throw exc;
    }
    x -> size = (b.end() - b.begin())/b.stride()+1;
    x -> data = gsl_matrix_ptr(A, i-1, b.begin());
    x -> stride = b.stride();
    x -> block = 0;
    x -> owner = 0;
  }

Here is the call graph for this function:

---

Member Function Documentation

vector_view & linalg::vector_view::operator=

(

const vector &

w

)

Assignment to a vector; will share same data in memory as that other vector.

Reimplemented from linalg::vector.

                                                    {
    if(x->size != w.size()){
      badArgument exc;
      exc.reason =  "Cannot assign to vector view: incomformant sizes.";
      exc.file = __FILE__;
      exc.line = __LINE__;
      throw exc;
    }
    gsl_vector_memcpy(x,w.x);
    return *this;
  }

Here is the call graph for this function:

vector_view & linalg::vector_view::operator=

(

const vector_view &

w

)

Assignment, useful for modifying chunks of vectors all at once.

                                                         {
    if(x->size != w.size()){
      badArgument exc;
      exc.reason =  "Cannot assign to vector view: incomformant sizes.";
      exc.file = __FILE__;
      exc.line = __LINE__;
      throw exc;
    }
    gsl_vector_memcpy(x,w.x);
    return *this;
    
  }

Here is the call graph for this function:

---

Friends And Related Function Documentation

friend class vector [friend]

friend class matrix [friend]

Reimplemented from linalg::vector.

---

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

• include/linalg.hpp
• linalg.cpp