2023-05-10 19:55:40 +08:00

90 lines
2.6 KiB
C

/**
* \file mp.h
*
* \brief multicore matrix operations
*
* \author Martin Albrecht <martinralbrecht@googlemail.com>
*/
#ifndef M4RI_MP_H
#define M4RI_MP_H
/*******************************************************************
*
* M4RI: Linear Algebra over GF(2)
*
* Copyright (C) 2014 Martin Albrecht <martinralbrecht@googlemail.com>
*
* Distributed under the terms of the GNU General Public License (GPL)
* version 2 or higher.
*
* This code 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.
*
* The full text of the GPL is available at:
*
* http://www.gnu.org/licenses/
*
********************************************************************/
#include <m4ri/mzd.h>
/**
* \brief Matrix multiplication via the cubic multiplication algorithm on
* multiple cores, i.e. compute C = AB.
*
* This is the wrapper function including bounds checks. See _mzd_mul_mp4 for
* implementation details.
*
* \param C Preallocated product matrix, may be NULL for automatic creation.
* \param A Input matrix A
* \param B Input matrix B
* \param cutoff Minimal dimension for recursion.
*/
mzd_t *mzd_mul_mp(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff);
/**
* \brief Matrix multiplication and in-place addition via the cubic matrix
* multiplication algorithm on multiple cores, i.e. compute C = C+ AB.
*
* This is the wrapper function including bounds checks. See _mzd_addmul_mp4 for
* implementation details.
*
* \param C product matrix
* \param A Input matrix A
* \param B Input matrix B
* \param cutoff Minimal dimension for recursion.
*/
mzd_t *mzd_addmul_mp(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff);
/**
* \brief Matrix multiplication and in-place addition via cubic matrix
* multiplication algorithm on up to four cores, i.e. compute C = C+ AB.
*
* \param C Product matrix
* \param A Input matrix A
* \param B Input matrix B
* \param cutoff Minimal dimension for recursion.
*/
mzd_t *_mzd_addmul_mp4(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff);
/**
* \brief Matrix multiplication via cubic matrix multiplication algorithm on up
* to four cores, i.e. compute C = C+ AB.
*
* \param C Product matrix
* \param A Input matrix A
* \param B Input matrix B
* \param cutoff Minimal dimension for recursion.
*/
mzd_t *_mzd_mul_mp4(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff);
#endif //__M4RI_HAVE_OPENMP