19#ifndef CCGSL_MATRIX_FLOAT_HPP
20#define CCGSL_MATRIX_FLOAT_HPP
22#include<gsl/gsl_matrix.h>
24#include<gsl/gsl_permute_matrix_float.h>
75 if( n1 > 0 and n2 > 0 )
ccgsl_pointer = gsl_matrix_float_alloc( n1, n2 );
79 try {
count =
new size_t; }
catch( std::bad_alloc& e ){
81 if( n1 > 0 and n2 > 0 ) gsl_matrix_float_free(
ccgsl_pointer );
98#ifdef __GXX_EXPERIMENTAL_CXX0X__
104 size_t const n1 = initializer_list.size();
105 size_t const n2 = initializer_list.begin()->size();
106 for(
auto l : initializer_list ){
107 if( l.size() != n2 ){
108 gsl::exception e(
"matrix_float rows have unequal sizes", __FILE__, __LINE__,
115 try {
count =
new size_t; }
catch( std::bad_alloc& e ){
117 if( n1 > 0 and n2 > 0 ) gsl_matrix_float_free(
ccgsl_pointer );
124 for(
auto row : initializer_list ){
126 for(
auto x :
row ){
set( r, c, x ); ++c; }
168 gsl_matrix_float_memcpy( copy.
get(),
get() );
213#ifdef __GXX_EXPERIMENTAL_CXX0X__
220 std::swap(
count, v.count );
221 v.ccgsl_pointer =
nullptr;
258 template<
typename container,
typename content,
bool reverse_t>
class iterator_base {
285 static content something;
290 }
else if(
v->ccgsl_pointer == 0 ){
309 static content something_base;
315 }
else if(
v->ccgsl_pointer == 0 ){
326#ifdef __GXX_EXPERIMENTAL_CXX0X__
327 return std::move(
ptr );
358 }
else if(
v->ccgsl_pointer == 0 ){
375 }
else if(
v->ccgsl_pointer == 0 ){
496 :
public iterator_base<matrix_float const,vector_float const,reverse_t>{
709 if(
ccgsl_pointer == 0 ) gsl_error(
"null vector_float", __FILE__, __LINE__, GSL_EFAULT );
710#ifndef GSL_RANGE_CHECK_OFF
712 gsl_error(
"matrix_float size2 and tda do not match", __FILE__, __LINE__, GSL_EBADLEN );
723 if(
ccgsl_pointer == 0 ) gsl_error(
"null vector_float", __FILE__, __LINE__, GSL_EFAULT );
724#ifndef GSL_RANGE_CHECK_OFF
726 gsl_error(
"matrix_float size2 and tda do not match", __FILE__, __LINE__, GSL_EBADLEN );
745 gsl_matrix_float_tricpy( Uplo, Diag,
get(), src.
get() );
754 gsl_matrix_float_transpose_tricpy( Uplo, Diag,
get(), src.
get() );
766 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
767 *m = gsl_matrix_float_submatrix(
get(), i, j, n1, n2 ).matrix;
776 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
777 *w = gsl_matrix_float_row(
get(), i ).vector;
786 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
787 *w = gsl_matrix_float_column(
get(), j ).vector;
795 diagonal(){ gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
796 *w = gsl_matrix_float_diagonal(
get() ).vector;
805 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
806 *w = gsl_matrix_float_subdiagonal(
get(), k ).vector;
815 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
816 *w = gsl_matrix_float_superdiagonal(
get(), k ).vector;
827 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
828 *w = gsl_matrix_float_subrow(
get(), i, offset,
n ).vector;
839 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
840 *w = gsl_matrix_float_subcolumn(
get(), j, offset,
n ).vector;
851 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
852 *m = gsl_matrix_float_view_array( base, n1, n2 ).matrix;
864 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
865 *m = gsl_matrix_float_view_array_with_tda( base, n1, n2, tda ).matrix;
876 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
877 *m = gsl_matrix_float_view_vector( v.
get(), n1, n2 ).matrix;
889 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
890 *m = gsl_matrix_float_view_vector_with_tda( v.
get(), n1, n2, tda ).matrix;
903 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
904 *m = gsl_matrix_float_const_submatrix(
get(), i, j, n1, n2 ).matrix;
913 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
914 *w = gsl_matrix_float_const_row(
get(), i ).vector;
923 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
924 *w = gsl_matrix_float_const_column(
get(), j ).vector;
932 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
933 *w = gsl_matrix_float_const_diagonal(
get() ).vector;
942 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
943 *w = gsl_matrix_float_const_subdiagonal(
get(), k ).vector;
952 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
953 *w = gsl_matrix_float_const_superdiagonal(
get(), k ).vector;
964 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
965 *w = gsl_matrix_float_const_subrow(
get(), i, offset,
n ).vector;
976 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
977 *w = gsl_matrix_float_const_subcolumn(
get(), j, offset,
n ).vector;
989 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
990 *m = gsl_matrix_float_const_submatrix(
get(), i, j, n1, n2 ).matrix;
999 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
1000 *w = gsl_matrix_float_const_row(
get(), i ).vector;
1009 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
1010 *w = gsl_matrix_float_const_column(
get(), j ).vector;
1018 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
1019 *w = gsl_matrix_float_const_diagonal(
get() ).vector;
1028 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
1029 *w = gsl_matrix_float_const_subdiagonal(
get(), k ).vector;
1038 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
1039 *w = gsl_matrix_float_const_superdiagonal(
get(), k ).vector;
1050 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
1051 *w = gsl_matrix_float_const_subrow(
get(), i, offset,
n ).vector;
1062 gsl_vector_float* w =
static_cast<gsl_vector_float*
>( malloc(
sizeof( gsl_vector_float ) ) );
1063 *w = gsl_matrix_float_const_subcolumn(
get(), j, offset,
n ).vector;
1074 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
1075 *m = gsl_matrix_float_const_view_array( base, n1, n2 ).matrix;
1088 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
1089 *m = gsl_matrix_float_const_view_array_with_tda( base, n1, n2, tda ).matrix;
1100 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
1101 *m = gsl_matrix_float_const_view_vector( v.
get(), n1, n2 ).matrix;
1114 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
1115 *m = gsl_matrix_float_const_view_vector_with_tda( v.
get(), n1, n2, tda ).matrix;
1126 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
1127 *m = gsl_matrix_float_const_view_array( base, n1, n2 ).matrix;
1140 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
1141 *m = gsl_matrix_float_const_view_array_with_tda( base, n1, n2, tda ).matrix;
1152 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
1153 *m = gsl_matrix_float_const_view_vector( v.
get(), n1, n2 ).matrix;
1166 gsl_matrix_float* m =
static_cast<gsl_matrix_float*
>( malloc(
sizeof( gsl_matrix_float ) ) );
1167 *m = gsl_matrix_float_const_view_vector_with_tda( v.
get(), n1, n2, tda ).matrix;
1211#ifdef __GXX_EXPERIMENTAL_CXX0X__
1243 float max()
const {
return gsl_matrix_float_max(
get() ); }
1248 float min()
const {
return gsl_matrix_float_min(
get() ); }
1254 void minmax(
float* min_out,
float* max_out )
const {
1255 gsl_matrix_float_minmax(
get(), min_out, max_out ); }
1261 void minmax(
float& min_out,
float& max_out )
const {
1262 gsl_matrix_float_minmax(
get(), &min_out, &max_out ); }
1280 int scale(
float const x ){
return gsl_matrix_float_scale(
get(), x ); }
1291 int isnull()
const {
return gsl_matrix_float_isnull(
get() ); }
1296 int ispos()
const {
return gsl_matrix_float_ispos(
get() ); }
1301 int isneg()
const {
return gsl_matrix_float_isneg(
get() ); }
1313 float get(
size_t const i,
size_t const j )
const {
return gsl_matrix_float_get(
get(), i, j ); }
1320 void set(
size_t const i,
size_t const j,
float x ){ gsl_matrix_float_set(
get(), i, j, x ); }
1327 float*
ptr(
size_t const i,
size_t const j ){
return gsl_matrix_float_ptr(
get(), i, j ); }
1334 float const*
const_ptr(
size_t const i,
size_t const j )
const {
1335 return gsl_matrix_float_const_ptr(
get(), i, j ); }
1341 int fread( FILE* stream ){
return gsl_matrix_float_fread( stream,
get() ); }
1347 int fwrite( FILE* stream )
const {
return gsl_matrix_float_fwrite( stream,
get() ); }
1353 int fscanf( FILE* stream ){
return gsl_matrix_float_fscanf( stream,
get() ); }
1360 int fprintf( FILE* stream,
char const* format )
const {
1361 return gsl_matrix_float_fprintf( stream,
get(), format ); }
1371 ccgsl_pointer = gsl_matrix_float_alloc_from_block(
b.get(), offset, n1, n2, d2 );
1373 try {
count =
new size_t; }
catch( std::bad_alloc& e ){
1389 ccgsl_pointer = gsl_matrix_float_alloc_from_matrix( m.
get(), k1, k2, n1, n2 );
1391 try {
count =
new size_t; }
catch( std::bad_alloc& e ){
1409 int swap_rows(
size_t const i,
size_t const j ){
return gsl_matrix_float_swap_rows(
get(), i, j ); }
1417 return gsl_matrix_float_swap_columns(
get(), i, j ); }
1424 int swap_rowcol(
size_t const i,
size_t const j ){
return gsl_matrix_float_swap_rowcol(
get(), i, j ); }
1436 return gsl_matrix_float_transpose_memcpy(
get(), src.
get() ); }
1443 void max_index(
size_t* imax,
size_t* jmax )
const {
1444 gsl_matrix_float_max_index(
get(), imax, jmax ); }
1450 void min_index(
size_t* imin,
size_t* jmin )
const {
1451 gsl_matrix_float_min_index(
get(), imin, jmin ); }
1459 void minmax_index(
size_t* imin,
size_t* jmin,
size_t* imax,
size_t* jmax )
const {
1460 gsl_matrix_float_minmax_index(
get(), imin, jmin, imax, jmax ); }
1468 gsl_matrix_float_max_index(
get(), &imax, &jmax ); }
1475 gsl_matrix_float_min_index(
get(), &imin, &jmin ); }
1483 void minmax_index(
size_t& imin,
size_t& jmin,
size_t& imax,
size_t& jmax )
const {
1484 gsl_matrix_float_minmax_index(
get(), &imin, &jmin, &imax, &jmax ); }
1492 return gsl_matrix_float_mul_elements(
get(),
b.get() ); }
1500 return gsl_matrix_float_div_elements(
get(),
b.get() ); }
1506 return gsl_matrix_float_norm1(
get() ); }
1513 return gsl_matrix_float_scale_rows(
get(), x.
get() ); }
1520 return gsl_matrix_float_scale_columns(
get(), x.
get() ); }
1527 return gsl_matrix_float_add_diagonal(
get(), x ); }
1535 return gsl_matrix_float_get_row( v.
get(),
get(), i ); }
1543 return gsl_matrix_float_get_col( v.
get(),
get(), j ); }
1551 return gsl_matrix_float_set_row(
get(), i, v.
get() ); }
1559 return gsl_matrix_float_set_col(
get(), j, v.
get() ); }
1574#ifndef GSL_RANGE_CHECK_OFF
1575#ifndef __GXX_EXPERIMENTAL_CXX0X__
1580 gsl_error(
"trying to read beyond last row of matrix_float",
1582#ifdef __GXX_EXPERIMENTAL_CXX0X__
1591 gsl_vector_float_view w = gsl_matrix_float_row(
ccgsl_pointer, i );
1609 gsl_vector_float_view w = gsl_matrix_float_row(
ccgsl_pointer, i );
1620 return gsl_permute_matrix_float( p.
get(),
get() ); }
1626 return vector_float ( gsl_vector_float_alloc_row_from_matrix( m.get(), i ) ); }
1628 return vector_float ( gsl_vector_float_alloc_col_from_matrix( m.get(), i ) ); }
This class handles vector_floats as shared handles.
This class is used to handle gsl exceptions so that gsl can use these rather than the GSL error handl...
@ GSL_EFAILED
generic failure
@ GSL_EINVAL
invalid argument supplied by user
@ GSL_EBADLEN
matrix, vector lengths are not conformant
@ GSL_EFAULT
invalid pointer
A class template for the const iterators.
const_iterator_t(iterator_t< reverse_t > const &i)
A copy constructor.
const_iterator_t< reverse_t > & operator=(const_iterator_t< reverse_t > const &i)
We can assign one output iterator from another.
const_iterator_t< reverse_t > & operator++()
The prefix ++ operator.
bool operator!=(const_iterator_t< reverse_t > const &i) const
Comparison with const iterator.
bool operator==(iterator_t< reverse_t > const &i) const
Comparison with non-const iterator.
bool operator!=(iterator_t< reverse_t > const &i) const
Comparison with non-const iterator.
const_iterator_t< reverse_t > operator++(int)
The postfix ++ operator.
const_iterator_t()
The default constructor.
const_iterator_t< reverse_t > operator--(int)
The postfix – operator.
bool operator==(const_iterator_t< reverse_t > const &i) const
Comparison with const iterator.
const_iterator_t< reverse_t > & operator--()
The prefix – operator.
const_iterator_t(matrix_float const *v, size_t position)
This constructor allows vector_float to create non-default iterators.
We create a suitable class for iterator types here.
iterator_base(container *v, size_t position)
This constructor allows vector_float to create non-default iterators.
container * v
Store a pointer to a matrix_float we can iterate over: 0 if no matrix_float.
std::bidirectional_iterator_tag iterator_category
An iterator must have an iterator category.
bool operator==(iterator_base< container, content, reverse_t > const &i) const
The == operator.
void increment()
Increment the iterator.
pointer operator->() const
Dereference the pointer.
bool operator!=(iterator_base< container, content, reverse_t > const &i) const
The != operator.
reference operator*() const
Dereference the pointer.
value_type reference
An iterator must have a reference type.
iterator_base()
The iterator is default constructible.
vector_float_ptr pointer
An iterator must have a pointer typea.
size_t position
Mark position of iterator within matrix_float.
vector_float value_type
An iterator must have a value type.
void decrement()
Decrement the iterator.
A class template for the two non-const iterators.
iterator_t()
The default constructor.
bool operator!=(const_iterator_t< reverse_t > const &i) const
Comparison with const iterator.
iterator_t< reverse_t > & operator=(iterator_t< reverse_t > const &i)
We can assign one output iterator from another.
iterator_t(matrix_float *v, size_t position)
This constructor allows vector_float to create non-default iterators.
iterator_t< reverse_t > & operator++()
The prefix ++ operator.
iterator_t< reverse_t > & operator--()
The prefix – operator.
iterator_t< reverse_t > operator--(int)
The postfix – operator.
bool operator==(const_iterator_t< reverse_t > const &i) const
Comparison with const iterator.
iterator_t< reverse_t > operator++(int)
The postfix ++ operator.
This class handles matrix_float objects as shared handles.
int scale_rows(vector_float const &x)
C++ version of gsl_matrix_float_scale_rows().
static matrix_float const view_array_with_tda(float const *base, size_t const n1, size_t const n2, size_t const tda)
Another C++ version of gsl_matrix_float_const_view_array_with_tda().
static matrix_float const view_vector(vector_float const &v, size_t const n1, size_t const n2)
Another C++ version of gsl_matrix_float_const_view_vector().
void reset()
Stop sharing ownership of the shared pointer.
void set(size_t const i, size_t const j, float x)
C++ version of gsl_matrix_float_set().
size_t use_count() const
Find how many matrix_float objects share this pointer.
~matrix_float()
The destructor only deletes the pointers if count reaches zero.
vector_float subrow(size_t const i, size_t const offset, size_t const n)
C++ version of gsl_matrix_float_subrow().
int isnull() const
C++ version of gsl_matrix_float_isnull().
int swap_rows(size_t const i, size_t const j)
C++ version of gsl_matrix_float_swap_rows().
int fscanf(FILE *stream)
C++ version of gsl_matrix_float_fscanf().
vector_float const const_row(size_t const i) const
C++ version of gsl_matrix_float_const_row().
const_iterator_t< false > const_iterator
The const_iterator type.
int permute(permutation &p)
Permute the columns of this by permutation p.
vector_float operator[](size_t const i)
This function allows us to use a matrix_float like an array.
static matrix_float const const_view_vector_with_tda(vector_float const &v, size_t const n1, size_t const n2, size_t const tda)
C++ version of gsl_matrix_float_const_view_vector_with_tda().
int fread(FILE *stream)
C++ version of gsl_matrix_float_fread().
int ispos() const
C++ version of gsl_matrix_float_ispos().
size_t size2() const
The number of columns of the matrix_float.
static matrix_float const view_array(float const *base, size_t const n1, size_t const n2)
Another C++ version of gsl_matrix_float_const_view_array().
vector_float const subcolumn(size_t const j, size_t const offset, size_t const n) const
Another C++ version of gsl_matrix_float_const_subcolumn().
float get(size_t const i, size_t const j) const
C++ version of gsl_matrix_float_get().
vector_float subdiagonal(size_t const k)
C++ version of gsl_matrix_float_subdiagonal().
size_t * count
The shared reference count.
matrix_float(gsl_matrix_float *v)
Could construct from a gsl_matrix_float.
reverse_iterator rend()
Get iterator pointing beyond last vector_float element.
void minmax(float &min_out, float &max_out) const
C++ version of gsl_matrix_float_minmax().
bool unique() const
Find if this is the only object sharing the gsl_matrix_float.
int isnonneg() const
C++ version of gsl_matrix_float_isnonneg().
vector_float column(size_t const j)
C++ version of gsl_matrix_float_column().
matrix_float()
The default constructor is only really useful for assigning to.
static matrix_float view_array_with_tda(float *base, size_t const n1, size_t const n2, size_t const tda)
C++ version of gsl_matrix_float_view_array_with_tda().
int set_col(size_t const j, vector_float const &v)
C++ version of gsl_matrix_float_set_col().
static matrix_float view_vector_with_tda(vector_float &v, size_t const n1, size_t const n2, size_t const tda)
C++ version of gsl_matrix_float_view_vector_with_tda().
int isneg() const
C++ version of gsl_matrix_float_isneg().
static matrix_float const const_view_array(float const *base, size_t const n1, size_t const n2)
C++ version of gsl_matrix_float_const_view_array().
matrix_float(matrix_float const &v)
The copy constructor.
int get_row(vector_float &v, size_t const i) const
C++ version of gsl_matrix_float_get_row().
static matrix_float const const_view_vector(vector_float const &v, size_t const n1, size_t const n2)
C++ version of gsl_matrix_float_const_view_vector().
vector_float const diagonal() const
Another C++ version of gsl_matrix_float_const_diagonal().
matrix_float clone() const
The clone function.
int memcpy(matrix_float const &src)
C++ version of gsl_matrix_float_memcpy().
matrix_float(matrix_float &m, size_t const k1, size_t const k2, size_t const n1, size_t const n2)
C++ version of gsl_matrix_float_alloc_from_matrix().
vector_float const const_subcolumn(size_t const j, size_t const offset, size_t const n) const
C++ version of gsl_matrix_float_const_subcolumn().
void tricpy(CBLAS_UPLO_t Uplo, CBLAS_DIAG_t Diag, matrix_float const &src)
Copy the upper or lower triangular part of matrix src to this.
const_reverse_iterator rend() const
Get iterator pointing beyond last vector_float element.
vector_float const subrow(size_t const i, size_t const offset, size_t const n) const
Another C++ version of gsl_matrix_float_const_subrow().
int get_col(vector_float &v, size_t const j) const
C++ version of gsl_matrix_float_get_col().
void max_index(size_t &imax, size_t &jmax) const
C++ version of gsl_matrix_float_max_index().
iterator_t< true > reverse_iterator
The reverse_iterator type.
int transpose()
C++ version of gsl_matrix_float_transpose().
int add_constant(float const x)
C++ version of gsl_matrix_float_add_constant().
matrix_float & operator=(matrix_float &&v)
Move operator.
int fwrite(FILE *stream) const
C++ version of gsl_matrix_float_fwrite().
int add_diagonal(float const x)
C++ version of gsl_matrix_float_add_diagonal().
int fprintf(FILE *stream, char const *format) const
C++ version of gsl_matrix_float_fprintf().
const_iterator begin() const
Get iterator pointing to first vector_float element.
gsl_matrix_float const * get() const
Get the gsl_matrix_float.
vector_float const row(size_t const i) const
Another C++ version of gsl_matrix_float_const_row().
const_iterator_t< true > const_reverse_iterator
The const_reverse_t type.
vector_float const const_superdiagonal(size_t const k) const
C++ version of gsl_matrix_float_const_superdiagonal().
void transpose_tricpy(CBLAS_UPLO_t Uplo, CBLAS_DIAG_t Diag, matrix_float const &src)
Copy the upper or lower triangular part of matrix src to this.
vector_float const const_diagonal() const
C++ version of gsl_matrix_float_const_diagonal().
int scale_columns(vector_float const &x)
C++ version of gsl_matrix_float_scale_columns().
vector_float const superdiagonal(size_t const k) const
Another C++ version of gsl_matrix_float_const_superdiagonal().
vector_float diagonal()
C++ version of gsl_matrix_float_diagonal().
vector_float const const_subrow(size_t const i, size_t const offset, size_t const n) const
C++ version of gsl_matrix_float_const_subrow().
matrix_float submatrix(size_t const i, size_t const j, size_t const n1, size_t const n2)
C++ version of gsl_matrix_float_submatrix().
matrix_float(std::initializer_list< std::initializer_list< float > > initializer_list)
Could construct from a std::initializer_list in C++11.
bool owns_data
Used to allow a vector that does not own its data.
reverse_iterator rbegin()
Get iterator pointing to first vector_float element.
static matrix_float view_vector(vector_float &v, size_t const n1, size_t const n2)
C++ version of gsl_matrix_float_view_vector().
vector_float const const_column(size_t const j) const
C++ version of gsl_matrix_float_const_column().
gsl_matrix_float * ccgsl_pointer
The shared pointer.
size_t size1() const
The number of rows of the matrix_float.
void minmax_index(size_t &imin, size_t &jmin, size_t &imax, size_t &jmax) const
C++ version of gsl_matrix_float_minmax_index().
gsl_matrix_float * get()
Get the gsl_matrix_float.
int mul_elements(matrix_float const &b)
C++ version of gsl_matrix_float_mul_elements().
static matrix_float view_array(float *base, size_t const n1, size_t const n2)
C++ version of gsl_matrix_float_view_array().
float * data()
Give access to the data block_float.
matrix_float(block_float &b, size_t const offset, size_t const n1, size_t const n2, size_t const d2)
C++ version of gsl_matrix_float_alloc_from_block().
vector_float row(size_t const i)
C++ version of gsl_matrix_float_row().
vector_float subcolumn(size_t const j, size_t const offset, size_t const n)
C++ version of gsl_matrix_float_subcolumn().
void minmax(float *min_out, float *max_out) const
C++ version of gsl_matrix_float_minmax().
iterator begin()
Get iterator pointing to first vector_float element.
void set_identity()
C++ version of gsl_matrix_float_set_identity().
vector_float const column(size_t const j) const
Another C++ version of gsl_matrix_float_const_column().
matrix_float const submatrix(size_t const i, size_t const j, size_t const n1, size_t const n2) const
Another C++ version of gsl_matrix_float_const_submatrix().
matrix_float const const_submatrix(size_t const i, size_t const j, size_t const n1, size_t const n2) const
C++ version of gsl_matrix_float_const_submatrix().
float * ptr(size_t const i, size_t const j)
C++ version of gsl_matrix_float_ptr().
int transpose_memcpy(matrix_float const &src)
C++ version of gsl_matrix_float_transpose_memcpy().
void wrap_gsl_matrix_float_without_ownership(gsl_matrix_float *v)
This function is intended mainly for internal use.
matrix_float & operator=(matrix_float const &v)
The assignment operator.
void set_zero()
C++ version of gsl_matrix_float_set_zero().
size_t size_type
A container must have a size_type.
vector_float const operator[](size_t const i) const
This function allows us to use a matrix_float like an array.
int swap_rowcol(size_t const i, size_t const j)
C++ version of gsl_matrix_float_swap_rowcol().
int swap_columns(size_t const i, size_t const j)
C++ version of gsl_matrix_float_swap_columns().
int div_elements(matrix_float const &b)
C++ version of gsl_matrix_float_div_elements().
static matrix_float const view_vector_with_tda(vector_float const &v, size_t const n1, size_t const n2, size_t const tda)
Another C++ version of gsl_matrix_float_const_view_vector_with_tda().
iterator_t< false > iterator
The iterator type.
int add(matrix_float const &b)
C++ version of gsl_matrix_float_add().
float const * const_ptr(size_t const i, size_t const j) const
C++ version of gsl_matrix_float_const_ptr().
matrix_float(matrix_float &&v)
Move constructor.
iterator end()
Get iterator pointing beyond last vector_float element.
int sub(matrix_float const &b)
C++ version of gsl_matrix_float_sub().
vector_float superdiagonal(size_t const k)
C++ version of gsl_matrix_float_superdiagonal().
static matrix_float const const_view_array_with_tda(float const *base, size_t const n1, size_t const n2, size_t const tda)
C++ version of gsl_matrix_float_const_view_array_with_tda().
matrix_float(size_t const n1, size_t const n2)
This constructor creates a new matrix_float with n1 rows and n2 columns.
float min() const
C++ version of gsl_matrix_float_min().
float norm1() const
C++ version of gsl_matrix_float_norm1().
float const * data() const
Give access to the data block_float.
vector_float const subdiagonal(size_t const k) const
Another C++ version of gsl_matrix_float_const_subdiagonal().
float max() const
C++ version of gsl_matrix_float_max().
int scale(float const x)
C++ version of gsl_matrix_float_scale().
static matrix_float calloc(size_t const n1, size_t const n2)
C++ version of gsl_matrix_float_calloc().
int set_row(size_t const i, vector_float const &v)
C++ version of gsl_matrix_float_set_row().
const_reverse_iterator rbegin() const
Get iterator pointing to first vector_float element.
void min_index(size_t &imin, size_t &jmin) const
C++ version of gsl_matrix_float_min_index().
void swap(matrix_float &m)
Swap two matrix_float objects.
const_iterator end() const
Get iterator pointing beyond last vector_float element.
vector_float const const_subdiagonal(size_t const k) const
C++ version of gsl_matrix_float_const_subdiagonal().
void set_all(float x)
C++ version of gsl_matrix_float_set_all().
This class handles GSL permutation objects.
gsl_permutation * get()
Get the gsl_permutation.
This class handles vector_float objects as shared handles.
static vector_float alloc_row_from_matrix(matrix_float &m, size_t const i)
C++ version of gsl_vector_float_alloc_row_from_matrix().
gsl_vector_float * get()
Get the gsl_vector_float.
vector_float()
The default constructor is only really useful for assigning to.
static vector_float alloc_col_from_matrix(matrix_float &m, size_t const j)
C++ version of gsl_vector_float_alloc_col_from_matrix().
void wrap_gsl_vector_float_without_ownership(gsl_vector_float *v)
This function is intended mainly for internal use.
size_t n(workspace const &w)
C++ version of gsl_rstat_n().
double b(int order, double qq)
C++ version of gsl_sf_mathieu_b().
gsl_sf_result result
Typedef for gsl_sf_result.
The gsl package creates an interface to the GNU Scientific Library for C++.
This is a pointer-like type for iterator return values.
vector_float_ptr(vector_float const &v)
Typically we have to construct from a vector_float.
vector_float & operator*()
Dereference operator.
vector_float * operator->()
Dereference operator.