1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83
|
//------------------------------------------------------------------------------
// GB_transpose_template: C=op(cast(A')), transpose, typecast, and apply op
//------------------------------------------------------------------------------
// SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2025, All Rights Reserved.
// SPDX-License-Identifier: Apache-2.0
//------------------------------------------------------------------------------
{
//--------------------------------------------------------------------------
// get A and C
//--------------------------------------------------------------------------
#undef GB_AB_TYPE
#ifdef GB_BIND_1ST
// A is the name of the matrix passed in to this kernel, but it takes
// the place of the B matrix for C=op(x,A'). As a result, the B macros
// must be used to access its contents.
#define GB_AB_TYPE GB_B_TYPE
#else
// for bind2nd, unary ops, and mere typecasting, use the A macros to
// access the A matrix.
#define GB_AB_TYPE GB_A_TYPE
#endif
#ifndef GB_ISO_TRANSPOSE
const GB_AB_TYPE *restrict Ax = (GB_AB_TYPE *) A->x ;
GB_C_TYPE *restrict Cx = (GB_C_TYPE *) C->x ;
#endif
//--------------------------------------------------------------------------
// C = op (cast (A'))
//--------------------------------------------------------------------------
#ifdef GB_JIT_KERNEL
{
#if GB_C_IS_FULL
{
// A and C are both full
#include "template/GB_transpose_full.c"
}
#elif GB_C_IS_BITMAP
{
// A and C are both bitmap
#include "template/GB_transpose_bitmap.c"
}
#else
{
// A is sparse or hypersparse; C is sparse
#include "template/GB_transpose_sparse.c"
}
#endif
}
#else
{
if (Workspaces == NULL)
{
// A and C are both full or both bitmap
if (A->b == NULL)
{
// A and C are both full
#include "template/GB_transpose_full.c"
}
else
{
// A and C are both bitmap
#include "template/GB_transpose_bitmap.c"
}
}
else
{
// A is sparse or hypersparse; C is sparse
#include "template/GB_transpose_sparse.c"
}
}
#endif
}
#undef GB_ISO_TRANSPOSE
#undef GB_AB_TYPE
|