summaryrefslogtreecommitdiff
path: root/2.3-1/thirdparty/includes/GSL/gsl/gsl_spmatrix.h
diff options
context:
space:
mode:
Diffstat (limited to '2.3-1/thirdparty/includes/GSL/gsl/gsl_spmatrix.h')
-rw-r--r--2.3-1/thirdparty/includes/GSL/gsl/gsl_spmatrix.h150
1 files changed, 0 insertions, 150 deletions
diff --git a/2.3-1/thirdparty/includes/GSL/gsl/gsl_spmatrix.h b/2.3-1/thirdparty/includes/GSL/gsl/gsl_spmatrix.h
deleted file mode 100644
index 741d9e69..00000000
--- a/2.3-1/thirdparty/includes/GSL/gsl/gsl_spmatrix.h
+++ /dev/null
@@ -1,150 +0,0 @@
-/* gsl_spmatrix.h
- *
- * Copyright (C) 2012-2014 Patrick Alken
- *
- * This program is free software; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 3 of the License, or (at
- * your option) any later version.
- *
- * This program is distributed in the hope that it will be useful, but
- * WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
- * General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
- */
-
-#ifndef __GSL_SPMATRIX_H__
-#define __GSL_SPMATRIX_H__
-
-#include <stdlib.h>
-
-#include <gsl/gsl_math.h>
-#include <gsl/gsl_vector.h>
-#include <gsl/gsl_matrix.h>
-
-#undef __BEGIN_DECLS
-#undef __END_DECLS
-#ifdef __cplusplus
-# define __BEGIN_DECLS extern "C" {
-# define __END_DECLS }
-#else
-# define __BEGIN_DECLS /* empty */
-# define __END_DECLS /* empty */
-#endif
-
-__BEGIN_DECLS
-
-/*
- * Binary tree data structure for storing sparse matrix elements
- * in triplet format. This is used for efficiently detecting
- * duplicates and element retrieval via gsl_spmatrix_get
- */
-typedef struct
-{
- void *tree; /* tree structure */
- void *node_array; /* preallocated array of tree nodes */
- size_t n; /* number of tree nodes in use (<= nzmax) */
-} gsl_spmatrix_tree;
-
-/*
- * Triplet format:
- *
- * If data[n] = A_{ij}, then:
- * i = A->i[n]
- * j = A->p[n]
- *
- * Compressed column format:
- *
- * If data[n] = A_{ij}, then:
- * i = A->i[n]
- * A->p[j] <= n < A->p[j+1]
- * so that column j is stored in
- * [ data[p[j]], data[p[j] + 1], ..., data[p[j+1] - 1] ]
- */
-
-typedef struct
-{
- size_t size1; /* number of rows */
- size_t size2; /* number of columns */
-
- size_t *i; /* row indices of size nzmax */
- double *data; /* matrix elements of size nzmax */
-
- /*
- * p contains the column indices (triplet) or column pointers (compcol)
- *
- * triplet: p[n] = column number of element data[n]
- * comp. col: p[j] = index in data of first non-zero element in column j
- * comp. row: p[i] = index in data of first non-zero element in row i
- */
- size_t *p;
-
- size_t nzmax; /* maximum number of matrix elements */
- size_t nz; /* number of non-zero values in matrix */
-
- gsl_spmatrix_tree *tree_data; /* binary tree for sorting triplet data */
-
- /*
- * workspace of size MAX(size1,size2)*MAX(sizeof(double),sizeof(size_t))
- * used in various routines
- */
- void *work;
-
- size_t sptype; /* sparse storage type */
-} gsl_spmatrix;
-
-#define GSL_SPMATRIX_TRIPLET (0)
-#define GSL_SPMATRIX_CCS (1)
-
-#define GSL_SPMATRIX_ISTRIPLET(m) ((m)->sptype == GSL_SPMATRIX_TRIPLET)
-#define GSL_SPMATRIX_ISCCS(m) ((m)->sptype == GSL_SPMATRIX_CCS)
-
-/*
- * Prototypes
- */
-
-gsl_spmatrix *gsl_spmatrix_alloc(const size_t n1, const size_t n2);
-gsl_spmatrix *gsl_spmatrix_alloc_nzmax(const size_t n1, const size_t n2,
- const size_t nzmax, const size_t flags);
-void gsl_spmatrix_free(gsl_spmatrix *m);
-int gsl_spmatrix_realloc(const size_t nzmax, gsl_spmatrix *m);
-int gsl_spmatrix_set_zero(gsl_spmatrix *m);
-size_t gsl_spmatrix_nnz(const gsl_spmatrix *m);
-int gsl_spmatrix_compare_idx(const size_t ia, const size_t ja,
- const size_t ib, const size_t jb);
-
-/* spcopy.c */
-int gsl_spmatrix_memcpy(gsl_spmatrix *dest, const gsl_spmatrix *src);
-
-/* spgetset.c */
-double gsl_spmatrix_get(const gsl_spmatrix *m, const size_t i,
- const size_t j);
-int gsl_spmatrix_set(gsl_spmatrix *m, const size_t i, const size_t j,
- const double x);
-
-/* spcompress.c */
-gsl_spmatrix *gsl_spmatrix_compcol(const gsl_spmatrix *T);
-void gsl_spmatrix_cumsum(const size_t n, size_t *c);
-
-/* spoper.c */
-int gsl_spmatrix_scale(gsl_spmatrix *m, const double x);
-int gsl_spmatrix_minmax(const gsl_spmatrix *m, double *min_out,
- double *max_out);
-int gsl_spmatrix_add(gsl_spmatrix *c, const gsl_spmatrix *a,
- const gsl_spmatrix *b);
-int gsl_spmatrix_d2sp(gsl_spmatrix *S, const gsl_matrix *A);
-int gsl_spmatrix_sp2d(gsl_matrix *A, const gsl_spmatrix *S);
-
-/* spprop.c */
-int gsl_spmatrix_equal(const gsl_spmatrix *a, const gsl_spmatrix *b);
-
-/* spswap.c */
-int gsl_spmatrix_transpose_memcpy(gsl_spmatrix *dest, const gsl_spmatrix *src);
-
-__END_DECLS
-
-#endif /* __GSL_SPMATRIX_H__ */