From 938fef4a37a7b7c61b4b6ff74cb4cfd2f100c427 Mon Sep 17 00:00:00 2001 From: Georgey Date: Wed, 5 Jul 2017 11:40:43 +0530 Subject: Added linux shared libraries and header files for int and ecos functions --- thirdparty/linux/include/coin/spla.h | 98 ++++++++++++++++++++++++++++++++++++ 1 file changed, 98 insertions(+) create mode 100644 thirdparty/linux/include/coin/spla.h (limited to 'thirdparty/linux/include/coin/spla.h') diff --git a/thirdparty/linux/include/coin/spla.h b/thirdparty/linux/include/coin/spla.h new file mode 100644 index 0000000..5e8e173 --- /dev/null +++ b/thirdparty/linux/include/coin/spla.h @@ -0,0 +1,98 @@ +/* + * ECOS - Embedded Conic Solver. + * Copyright (C) 2012-2015 A. Domahidi [domahidi@embotech.com], + * Automatic Control Lab, ETH Zurich & embotech GmbH, Zurich, Switzerland. + * + * 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, see . + */ + +/* + * Sparse linear algebra library for solver, i.e. no memory manager + * such as malloc is accessed by this module. + */ + +#ifndef __SPLA_H__ +#define __SPLA_H__ + +#include "glblopts.h" + +/* Data structure for sparse matrices */ +typedef struct spmat{ + idxint* jc; + idxint* ir; + pfloat* pr; + idxint n; + idxint m; + idxint nnz; +} spmat; + + +/* SPARSE MATRIX OPERATIONS PROVIDED BY THIS MODULE -------------------- */ + +/* + * Sparse matrix-vector multiply for operations + * + * y = A*x (if a > 0 && newVector == 1) + * y += A*x (if a > 0 && newVector == 0) + * y = -A*x (if a < 0 && newVector == 1) + * y -= A*x (if a < 0 && newVector == 0) + * + * where A is a sparse matrix and both x and y are assumed to be dense. + */ +void sparseMV(spmat* A, pfloat* x, pfloat* y, idxint a, idxint newVector); + + +/* + * Sparse matrix-transpose-vector multiply with subtraction. + * + * If newVector > 0, then this computes y = -A'*x, + * otherwise y -= A'*x, + * + * where A is a sparse matrix and both x and y are assumed to be dense. + * If skipDiagonal == 1, then the contributions of diagonal elements are + * not counted. + * + * NOTE: The product is calculating without explicitly forming the + * transpose. + */ +void sparseMtVm(spmat* A, pfloat* x, pfloat* y, idxint newVector, idxint skipDiagonal); + + +/* + * Vector addition y += x of size n. + */ +void vadd(idxint n, pfloat* x, pfloat* y); + +/* + * Vector subtraction with scaling: y -= a*x of size n. + */ +void vsubscale(idxint n, pfloat a, pfloat* x, pfloat* y); + +/* + * 2-norm of a vector. + */ +pfloat norm2(pfloat* v, idxint n); + +/* + * inf-norm of a vector. + */ +pfloat norminf(pfloat* v, idxint n); + +/* + * ECOS dot product z = x'*y of size n. + */ +pfloat eddot(idxint n, pfloat* x, pfloat* y); + + +#endif -- cgit