matblock.h

Go to the documentation of this file.

/*
 * This work is licensed under a Creative Commons 
 * Attribution-Noncommercial-Share Alike 3.0 United States License.
 * 
 *    http://creativecommons.org/licenses/by-nc-sa/3.0/us/
 * 
 * You are free:
 * 
 *    to Share - to copy, distribute, display, and perform the work
 *    to Remix - to make derivative works
 * 
 * Under the following conditions:
 * 
 *    Attribution. You must attribute the work in the manner specified by the
 *    author or licensor (but not in any way that suggests that they endorse you
 *    or your use of the work).
 * 
 *    Noncommercial. You may not use this work for commercial purposes.
 * 
 *    Share Alike. If you alter, transform, or build upon this work, you may
 *    distribute the resulting work only under the same or similar license to
 *    this one.
 * 
 * For any reuse or distribution, you must make clear to others the license
 * terms of this work. The best way to do this is by including this header.
 * 
 * Any of the above conditions can be waived if you get permission from the
 * copyright holder.
 * 
 * Apart from the remix rights granted under this license, nothing in this
 * license impairs or restricts the author's moral rights.
 */


#ifndef MATBLOCK_H 
#define MATBLOCK_H 


#include <jwsc/config.h>

#include <stdlib.h>
#include <inttypes.h>

#include "jwsc/base/error.h"
#include "jwsc/math/complex.h"
#include "jwsc/matrix/matrix.h"


#ifdef __cplusplus
namespace jwsc {
extern "C" {
#endif


typedef struct
{
    uint32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    int8_t*** elts;
}
Matblock_i8;

typedef struct
{
    uint32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    int16_t*** elts;
}
Matblock_i16;

typedef struct
{
    uint32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    int32_t*** elts;
}
Matblock_i32;

typedef struct
{
    uint32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    int64_t*** elts;
}
Matblock_i64;


typedef struct
{
    int32_t num_mats;

    int32_t num_rows;

    int32_t num_cols;

    uint8_t*** elts;
}
Matblock_u8;

typedef struct
{
    int32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    uint16_t*** elts;
}
Matblock_u16;

typedef struct
{
    int32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    uint32_t*** elts;
}
Matblock_u32;

typedef struct
{
    int32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    uint64_t*** elts;
}
Matblock_u64;

typedef struct
{
    uint32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    float*** elts;
}
Matblock_f;

typedef struct
{
    uint32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    double*** elts;
}
Matblock_d;

typedef struct
{
    uint32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    Complex_f*** elts;
}
Matblock_cf;

typedef struct
{
    uint32_t num_mats;

    uint32_t num_rows;

    uint32_t num_cols;

    Complex_d*** elts;
}
Matblock_cd;


typedef enum
{
    MATBLOCK_MAT_MATRIX,
    MATBLOCK_ROW_MATRIX,
    MATBLOCK_COL_MATRIX
}
Matblock_matrix;


void create_matblock_u8
(
    Matblock_u8** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols
);

void create_matblock_u32
(
    Matblock_u32** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols
);

void create_matblock_i32
(
    Matblock_i32** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols
);

void create_matblock_i64
(
    Matblock_i64** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols
);

void create_matblock_f
(
    Matblock_f** m_out, 
    uint32_t     num_mats, 
    uint32_t     num_rows, 
    uint32_t     num_cols
);

void create_matblock_d
(
    Matblock_d** m_out, 
    uint32_t     num_mats, 
    uint32_t     num_rows, 
    uint32_t     num_cols
);

void create_matblock_cf
(
    Matblock_cf** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols
);

void create_matblock_cd
(
    Matblock_cd** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols
);

void create_init_matblock_u8
(
    Matblock_u8** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols,
    uint8_t       val
);

void create_init_matblock_u32
(
    Matblock_u32** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols,
    uint32_t       val
);

void create_init_matblock_i32
(
    Matblock_i32** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols,
    int32_t        val
);

void create_init_matblock_i64
(
    Matblock_i64** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols,
    int64_t        val
);

void create_init_matblock_f
(
    Matblock_f** m_out, 
    uint32_t     num_mats, 
    uint32_t     num_rows, 
    uint32_t     num_cols,
    float        val
);

void create_init_matblock_d
(
    Matblock_d** m_out, 
    uint32_t     num_mats, 
    uint32_t     num_rows, 
    uint32_t     num_cols,
    double       val
);

void create_init_matblock_cf
(
    Matblock_cf** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols,
    Complex_f     val
);

void create_init_matblock_cd
(
    Matblock_cd** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols,
    Complex_d     val
);

void create_zero_matblock_u8
(
    Matblock_u8** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols
);

void create_zero_matblock_u32
(
    Matblock_u32** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols
);

void create_zero_matblock_i32
(
    Matblock_i32** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols
);

void create_zero_matblock_i64
(
    Matblock_i64** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols
);

void create_zero_matblock_f
(
    Matblock_f** m_out, 
    uint32_t     num_mats, 
    uint32_t     num_rows, 
    uint32_t     num_cols
);

void create_zero_matblock_d
(
    Matblock_d** m_out, 
    uint32_t     num_mats, 
    uint32_t     num_rows, 
    uint32_t     num_cols
);

void create_zero_matblock_cf
(
    Matblock_cf** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols
);

void create_zero_matblock_cd
(
    Matblock_cd** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols
);

void create_random_matblock_u32
(
    Matblock_u32** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols,
    uint32_t       min,
    uint32_t       max
);

void create_random_matblock_i32
(
    Matblock_i32** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols,
    int32_t        min,
    int32_t        max
);

void create_random_matblock_i64
(
    Matblock_i64** m_out, 
    uint32_t       num_mats, 
    uint32_t       num_rows, 
    uint32_t       num_cols,
    int64_t        min,
    int64_t        max
);

void create_random_matblock_f
(
    Matblock_f** m_out, 
    uint32_t     num_mats, 
    uint32_t     num_rows, 
    uint32_t     num_cols,
    float        min,
    float        max
);

void create_random_matblock_d
(
    Matblock_d** m_out, 
    uint32_t     num_mats, 
    uint32_t     num_rows, 
    uint32_t     num_cols,
    double       min,
    double       max
);

void create_random_matblock_cf
(
    Matblock_cf** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols,
    Complex_f     min,
    Complex_f     max
);

void create_random_matblock_cd
(
    Matblock_cd** m_out, 
    uint32_t      num_mats, 
    uint32_t      num_rows, 
    uint32_t      num_cols,
    Complex_d     min,
    Complex_d     max
);

Error* copy_matblock_u32(Matblock_u32** m_out, const Matblock_u32* m_in);

Error* copy_matblock_i32(Matblock_i32** m_out, const Matblock_i32* m_in);

Error* copy_matblock_i64(Matblock_i64** m_out, const Matblock_i64* m_in);

Error* copy_matblock_f(Matblock_f** m_out, const Matblock_f* m_in);

Error* copy_matblock_d(Matblock_d** m_out, const Matblock_d* m_in);

Error* copy_matblock_cf(Matblock_cf** m_out, const Matblock_cf* m_in);

Error* copy_matblock_cd(Matblock_cd** m_out, const Matblock_cd* m_in);

Error* copy_matblock_block_u32
(
    Matblock_u32**      m_out, 
    const Matblock_u32* m_in,
    uint32_t            mat_offset,
    uint32_t            row_offset,
    uint32_t            col_offset,
    uint32_t            num_mats,
    uint32_t            num_rows,
    uint32_t            num_cols
);

Error* copy_matblock_block_i32
(
    Matblock_i32**      m_out, 
    const Matblock_i32* m_in,
    uint32_t            mat_offset,
    uint32_t            row_offset,
    uint32_t            col_offset,
    uint32_t            num_mats,
    uint32_t            num_rows,
    uint32_t            num_cols
);

Error* copy_matblock_block_i64
(
    Matblock_i64**      m_out, 
    const Matblock_i64* m_in,
    uint32_t            mat_offset,
    uint32_t            row_offset,
    uint32_t            col_offset,
    uint32_t            num_mats,
    uint32_t            num_rows,
    uint32_t            num_cols
);

Error* copy_matblock_block_f
(
    Matblock_f**      m_out, 
    const Matblock_f* m_in,
    uint32_t          mat_offset,
    uint32_t          row_offset,
    uint32_t          col_offset,
    uint32_t          num_mats,
    uint32_t          num_rows,
    uint32_t          num_cols
);

Error* copy_matblock_block_d
(
    Matblock_d**      m_out, 
    const Matblock_d* m_in,
    uint32_t          mat_offset,
    uint32_t          row_offset,
    uint32_t          col_offset,
    uint32_t          num_mats,
    uint32_t          num_rows,
    uint32_t          num_cols
);

Error* copy_matblock_block_cf
(
    Matblock_cf**      m_out, 
    const Matblock_cf* m_in,
    uint32_t           mat_offset,
    uint32_t           row_offset,
    uint32_t           col_offset,
    uint32_t           num_mats,
    uint32_t           num_rows,
    uint32_t           num_cols
);

Error* copy_matblock_block_cd
(
    Matblock_cd**      m_out, 
    const Matblock_cd* m_in,
    uint32_t           mat_offset,
    uint32_t           row_offset,
    uint32_t           col_offset,
    uint32_t           num_mats,
    uint32_t           num_rows,
    uint32_t           num_cols
);

Error* copy_matblock_block_into_matblock_u32
(
    Matblock_u32*       m_1, 
    uint32_t            mat_offset_1,
    uint32_t            row_offset_1,
    uint32_t            col_offset_1,
    const Matblock_u32* m_2,
    uint32_t            mat_offset_2,
    uint32_t            row_offset_2,
    uint32_t            col_offset_2,
    uint32_t            num_mats,
    uint32_t            num_rows,
    uint32_t            num_cols
);

Error* copy_matblock_block_into_matblock_i32
(
    Matblock_i32*       m_1, 
    uint32_t            mat_offset_1,
    uint32_t            row_offset_1,
    uint32_t            col_offset_1,
    const Matblock_i32* m_2,
    uint32_t            mat_offset_2,
    uint32_t            row_offset_2,
    uint32_t            col_offset_2,
    uint32_t            num_mats,
    uint32_t            num_rows,
    uint32_t            num_cols
);

Error* copy_matblock_block_into_matblock_i64
(
    Matblock_i64*       m_1, 
    uint32_t            mat_offset_1,
    uint32_t            row_offset_1,
    uint32_t            col_offset_1,
    const Matblock_i64* m_2,
    uint32_t            mat_offset_2,
    uint32_t            row_offset_2,
    uint32_t            col_offset_2,
    uint32_t            num_mats,
    uint32_t            num_rows,
    uint32_t            num_cols
);

Error* copy_matblock_block_into_matblock_f
(
    Matblock_f*       m_1, 
    uint32_t          mat_offset_1,
    uint32_t          row_offset_1,
    uint32_t          col_offset_1,
    const Matblock_f* m_2,
    uint32_t          mat_offset_2,
    uint32_t          row_offset_2,
    uint32_t          col_offset_2,
    uint32_t          num_mats,
    uint32_t          num_rows,
    uint32_t          num_cols
);

Error* copy_matblock_block_into_matblock_d
(
    Matblock_d*       m_1, 
    uint32_t          mat_offset_1,
    uint32_t          row_offset_1,
    uint32_t          col_offset_1,
    const Matblock_d* m_2,
    uint32_t          mat_offset_2,
    uint32_t          row_offset_2,
    uint32_t          col_offset_2,
    uint32_t          num_mats,
    uint32_t          num_rows,
    uint32_t          num_cols
);

Error* copy_matblock_block_into_matblock_cf
(
    Matblock_cf*       m_1, 
    uint32_t           mat_offset_1,
    uint32_t           row_offset_1,
    uint32_t           col_offset_1,
    const Matblock_cf* m_2,
    uint32_t           mat_offset_2,
    uint32_t           row_offset_2,
    uint32_t           col_offset_2,
    uint32_t           num_mats,
    uint32_t           num_rows,
    uint32_t           num_cols
);

Error* copy_matblock_block_into_matblock_cd
(
    Matblock_cd*       m_1, 
    uint32_t           mat_offset_1,
    uint32_t           row_offset_1,
    uint32_t           col_offset_1,
    const Matblock_cd* m_2,
    uint32_t           mat_offset_2,
    uint32_t           row_offset_2,
    uint32_t           col_offset_2,
    uint32_t           num_mats,
    uint32_t           num_rows,
    uint32_t           num_cols
);

Error* copy_matrix_into_matblock_u32
(
    Matblock_u32**      mb_out,
    const Matblock_u32* mb_in,
    const Matrix_u32*   m,
    Matblock_matrix     orient,
    uint32_t            index
);

Error* copy_matrix_into_matblock_i32
(
    Matblock_i32**      mb_out,
    const Matblock_i32* mb_in,
    const Matrix_i32*   m,
    Matblock_matrix     orient,
    uint32_t            index
);

Error* copy_matrix_into_matblock_i64
(
    Matblock_i64**      mb_out,
    const Matblock_i64* mb_in,
    const Matrix_i64*   m,
    Matblock_matrix     orient,
    uint32_t            index
);

Error* copy_matrix_into_matblock_f
(
    Matblock_f**      mb_out,
    const Matblock_f* mb_in,
    const Matrix_f*   m,
    Matblock_matrix   orient,
    uint32_t          index
);

Error* copy_matrix_into_matblock_d
(
    Matblock_d**      mb_out,
    const Matblock_d* mb_in,
    const Matrix_d*   m,
    Matblock_matrix   orient,
    uint32_t          index
);

Error* copy_matrix_into_matblock_cf
(
    Matblock_cf**      mb_out,
    const Matblock_cf* mb_in,
    const Matrix_cf*   m,
    Matblock_matrix    orient,
    uint32_t           index
);

Error* copy_matrix_into_matblock_cd
(
    Matblock_cd**      mb_out,
    const Matblock_cd* mb_in,
    const Matrix_cd*   m,
    Matblock_matrix    orient,
    uint32_t           index
);

Error* copy_matrix_from_matblock_u32
(
    Matrix_u32**        m_out,
    const Matblock_u32* mb,
    Matblock_matrix     orient,
    uint32_t            index
);

Error* copy_matrix_from_matblock_i32
(
    Matrix_i32**        m_out,
    const Matblock_i32* mb,
    Matblock_matrix     orient,
    uint32_t            index
);

Error* copy_matrix_from_matblock_i64
(
    Matrix_i64**        m_out,
    const Matblock_i64* mb,
    Matblock_matrix     orient,
    uint32_t            index
);

Error* copy_matrix_from_matblock_f
(
    Matrix_f**        m_out,
    const Matblock_f* mb,
    Matblock_matrix   orient,
    uint32_t          index
);

Error* copy_matrix_from_matblock_d
(
    Matrix_d**        m_out,
    const Matblock_d* mb,
    Matblock_matrix   orient,
    uint32_t          index
);

Error* copy_matrix_from_matblock_cf
(
    Matrix_cf**        m_out,
    const Matblock_cf* mb,
    Matblock_matrix    orient,
    uint32_t           index
);

Error* copy_matrix_from_matblock_cd
(
    Matrix_cd**        m_out,
    const Matblock_cd* mb,
    Matblock_matrix    orient,
    uint32_t           index
);

void free_matblock_u8(Matblock_u8* m);

void free_matblock_u32(Matblock_u32* m);

void free_matblock_i32(Matblock_i32* m);

void free_matblock_i64(Matblock_i64* m);

void free_matblock_f(Matblock_f* m);

void free_matblock_d(Matblock_d* m);

void free_matblock_cf(Matblock_cf* m);

void free_matblock_cd(Matblock_cd* m);

#ifdef __cplusplus
}
}
#endif


#endif