summaryrefslogtreecommitdiff
path: root/gsl-1.9/matrix/view_source.c
diff options
context:
space:
mode:
Diffstat (limited to 'gsl-1.9/matrix/view_source.c')
-rw-r--r--gsl-1.9/matrix/view_source.c273
1 files changed, 273 insertions, 0 deletions
diff --git a/gsl-1.9/matrix/view_source.c b/gsl-1.9/matrix/view_source.c
new file mode 100644
index 0000000..f5bcc68
--- /dev/null
+++ b/gsl-1.9/matrix/view_source.c
@@ -0,0 +1,273 @@
+/* matrix/view_source.c
+ *
+ * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2001 Gerard Jungman, Brian Gough
+ *
+ * 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 2 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.
+ */
+
+QUALIFIED_VIEW (_gsl_matrix,view)
+FUNCTION (gsl_matrix, view_array) (QUALIFIER ATOMIC * array,
+ const size_t n1, const size_t n2)
+{
+ QUALIFIED_VIEW (_gsl_matrix,view) view = NULL_MATRIX_VIEW;
+
+ if (n1 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n1 must be positive integer",
+ GSL_EINVAL, view);
+ }
+ else if (n2 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n2 must be positive integer",
+ GSL_EINVAL, view);
+ }
+
+ {
+ TYPE(gsl_matrix) m = NULL_MATRIX;
+
+ m.data = (ATOMIC *)array;
+ m.size1 = n1;
+ m.size2 = n2;
+ m.tda = n2;
+ m.block = 0;
+ m.owner = 0;
+
+ view.matrix = m;
+ return view;
+ }
+}
+
+QUALIFIED_VIEW (_gsl_matrix,view)
+FUNCTION(gsl_matrix, view_array_with_tda) (QUALIFIER ATOMIC * base,
+ const size_t n1,
+ const size_t n2,
+ const size_t tda)
+{
+ QUALIFIED_VIEW (_gsl_matrix,view) view = NULL_MATRIX_VIEW;
+
+ if (n1 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n1 must be positive integer",
+ GSL_EINVAL, view);
+ }
+ else if (n2 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n2 must be positive integer",
+ GSL_EINVAL, view);
+ }
+ else if (n2 > tda)
+ {
+ GSL_ERROR_VAL ("matrix dimension n2 must not exceed tda",
+ GSL_EINVAL, view);
+ }
+
+
+ {
+ TYPE(gsl_matrix) m = NULL_MATRIX;
+
+ m.data = (ATOMIC *)base;
+ m.size1 = n1;
+ m.size2 = n2;
+ m.tda = tda;
+ m.block = 0;
+ m.owner = 0;
+
+ view.matrix = m;
+ return view;
+ }
+}
+
+QUALIFIED_VIEW (_gsl_matrix,view)
+FUNCTION(gsl_matrix, view_vector) (QUALIFIED_TYPE(gsl_vector) * v,
+ const size_t n1,
+ const size_t n2)
+{
+ QUALIFIED_VIEW (_gsl_matrix,view) view = NULL_MATRIX_VIEW;
+
+ if (n1 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n1 must be positive integer",
+ GSL_EINVAL, view);
+ }
+ else if (n2 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n2 must be positive integer",
+ GSL_EINVAL, view);
+ }
+ else if (v->stride != 1)
+ {
+ GSL_ERROR_VAL ("vector must have unit stride",
+ GSL_EINVAL, view);
+ }
+ else if (n1 * n2 > v->size)
+ {
+ GSL_ERROR_VAL ("matrix size exceeds size of original",
+ GSL_EINVAL, view);
+ }
+
+ {
+ TYPE(gsl_matrix) m = NULL_MATRIX;
+
+ m.data = v->data;
+ m.size1 = n1;
+ m.size2 = n2;
+ m.tda = n2;
+ m.block = v->block;
+ m.owner = 0;
+
+ view.matrix = m;
+ return view;
+ }
+}
+
+
+QUALIFIED_VIEW (_gsl_matrix,view)
+FUNCTION(gsl_matrix, view_vector_with_tda) (QUALIFIED_TYPE(gsl_vector) * v,
+ const size_t n1,
+ const size_t n2,
+ const size_t tda)
+{
+ QUALIFIED_VIEW (_gsl_matrix,view) view = NULL_MATRIX_VIEW;
+
+ if (n1 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n1 must be positive integer",
+ GSL_EINVAL, view);
+ }
+ else if (n2 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n2 must be positive integer",
+ GSL_EINVAL, view);
+ }
+ else if (v->stride != 1)
+ {
+ GSL_ERROR_VAL ("vector must have unit stride",
+ GSL_EINVAL, view);
+ }
+ else if (n2 > tda)
+ {
+ GSL_ERROR_VAL ("matrix dimension n2 must not exceed tda",
+ GSL_EINVAL, view);
+ }
+ else if (n1 * tda > v->size)
+ {
+ GSL_ERROR_VAL ("matrix size exceeds size of original",
+ GSL_EINVAL, view);
+ }
+
+ {
+ TYPE(gsl_matrix) m = NULL_MATRIX;
+
+ m.data = v->data;
+ m.size1 = n1;
+ m.size2 = n2;
+ m.tda = tda;
+ m.block = v->block;
+ m.owner = 0;
+
+ view.matrix = m;
+ return view;
+ }
+}
+
+#ifdef JUNK
+int
+FUNCTION (gsl_matrix, view_from_matrix) (TYPE(gsl_matrix) * m,
+ TYPE(gsl_matrix) * mm,
+ const size_t k1,
+ const size_t k2,
+ const size_t n1,
+ const size_t n2)
+{
+ if (n1 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n1 must be positive integer",
+ GSL_EINVAL, 0);
+ }
+ else if (n2 == 0)
+ {
+ GSL_ERROR_VAL ("matrix dimension n2 must be positive integer",
+ GSL_EINVAL, 0);
+ }
+ else if (k1 + n1 > mm->size1)
+ {
+ GSL_ERROR_VAL ("submatrix dimension 1 exceeds size of original",
+ GSL_EINVAL, 0);
+ }
+ else if (k2 + n2 > mm->size2)
+ {
+ GSL_ERROR_VAL ("submatrix dimension 2 exceeds size of original",
+ GSL_EINVAL, 0);
+ }
+
+ m->data = mm->data + k1 * mm->tda + k2 ;
+ m->size1 = n1;
+ m->size2 = n2;
+ m->tda = mm->tda;
+ m->block = mm->block;
+ m->owner = 0;
+
+ return GSL_SUCCESS;
+}
+
+int
+FUNCTION (gsl_vector, view_row_from_matrix) (TYPE(gsl_vector) * v,
+ TYPE(gsl_matrix) * m,
+ const size_t i)
+{
+ const size_t column_length = m->size1;
+
+ if (i >= column_length)
+ {
+ GSL_ERROR ("row index is out of range", GSL_EINVAL);
+ }
+
+ if (v->block != 0)
+ {
+ GSL_ERROR ("vector already has memory allocated to it", GSL_ENOMEM);
+ }
+
+ v->data = m->data + MULTIPLICITY * i * m->tda ;
+ v->size = m->size2;
+ v->stride = 1;
+
+ return GSL_SUCCESS;
+}
+
+int
+FUNCTION (gsl_vector, view_col_from_matrix) (TYPE(gsl_vector) * v,
+ TYPE(gsl_matrix) * m,
+ const size_t j)
+{
+ const size_t row_length = m->size2;
+
+ if (j >= row_length)
+ {
+ GSL_ERROR_VAL ("column index is out of range", GSL_EINVAL, 0);
+ }
+
+ if (v->block != 0)
+ {
+ GSL_ERROR ("vector already has memory allocated to it", GSL_ENOMEM);
+ }
+
+ v->data = m->data + MULTIPLICITY * j ;
+ v->size = m->size1;
+ v->stride = m->tda;
+
+ return GSL_SUCCESS;
+}
+#endif /* JUNK */
+
generated by cgit v1.2.3 (git 2.25.1) at 2024-05-01 00:37:21 +0000