CMSIS DSP Library from CMSIS 2.0. See http://www.onarm.com/cmsis/ for full details

Dependents:   K22F_DSP_Matrix_least_square BNO055-ELEC3810 1BNO055 ECE4180Project--Slave2 ... more

Committer:
simon
Date:
Thu Mar 10 15:07:50 2011 +0000
Revision:
0:1014af42efd9

        

Who changed what in which revision?

UserRevisionLine numberNew contents of line
simon 0:1014af42efd9 1 /* ----------------------------------------------------------------------
simon 0:1014af42efd9 2 * Copyright (C) 2010 ARM Limited. All rights reserved.
simon 0:1014af42efd9 3 *
simon 0:1014af42efd9 4 * $Date: 29. November 2010
simon 0:1014af42efd9 5 * $Revision: V1.0.3
simon 0:1014af42efd9 6 *
simon 0:1014af42efd9 7 * Project: CMSIS DSP Library
simon 0:1014af42efd9 8 * Title: arm_mat_trans_q31.c
simon 0:1014af42efd9 9 *
simon 0:1014af42efd9 10 * Description: Q31 matrix transpose.
simon 0:1014af42efd9 11 *
simon 0:1014af42efd9 12 * Target Processor: Cortex-M4/Cortex-M3
simon 0:1014af42efd9 13 *
simon 0:1014af42efd9 14 * Version 1.0.3 2010/11/29
simon 0:1014af42efd9 15 * Re-organized the CMSIS folders and updated documentation.
simon 0:1014af42efd9 16 *
simon 0:1014af42efd9 17 * Version 1.0.2 2010/11/11
simon 0:1014af42efd9 18 * Documentation updated.
simon 0:1014af42efd9 19 *
simon 0:1014af42efd9 20 * Version 1.0.1 2010/10/05
simon 0:1014af42efd9 21 * Production release and review comments incorporated.
simon 0:1014af42efd9 22 *
simon 0:1014af42efd9 23 * Version 1.0.0 2010/09/20
simon 0:1014af42efd9 24 * Production release and review comments incorporated.
simon 0:1014af42efd9 25 *
simon 0:1014af42efd9 26 * Version 0.0.5 2010/04/26
simon 0:1014af42efd9 27 * incorporated review comments and updated with latest CMSIS layer
simon 0:1014af42efd9 28 *
simon 0:1014af42efd9 29 * Version 0.0.3 2010/03/10
simon 0:1014af42efd9 30 * Initial version
simon 0:1014af42efd9 31 * -------------------------------------------------------------------- */
simon 0:1014af42efd9 32
simon 0:1014af42efd9 33 #include "arm_math.h"
simon 0:1014af42efd9 34
simon 0:1014af42efd9 35 /**
simon 0:1014af42efd9 36 * @ingroup groupMatrix
simon 0:1014af42efd9 37 */
simon 0:1014af42efd9 38
simon 0:1014af42efd9 39 /**
simon 0:1014af42efd9 40 * @addtogroup MatrixTrans
simon 0:1014af42efd9 41 * @{
simon 0:1014af42efd9 42 */
simon 0:1014af42efd9 43
simon 0:1014af42efd9 44 /*
simon 0:1014af42efd9 45 * @brief Q31 matrix transpose.
simon 0:1014af42efd9 46 * @param[in] *pSrc points to the input matrix
simon 0:1014af42efd9 47 * @param[out] *pDst points to the output matrix
simon 0:1014af42efd9 48 * @return The function returns either <code>ARM_MATH_SIZE_MISMATCH</code>
simon 0:1014af42efd9 49 * or <code>ARM_MATH_SUCCESS</code> based on the outcome of size checking.
simon 0:1014af42efd9 50 */
simon 0:1014af42efd9 51
simon 0:1014af42efd9 52 arm_status arm_mat_trans_q31(
simon 0:1014af42efd9 53 const arm_matrix_instance_q31 * pSrc,
simon 0:1014af42efd9 54 arm_matrix_instance_q31 * pDst)
simon 0:1014af42efd9 55 {
simon 0:1014af42efd9 56 q31_t *pIn = pSrc->pData; /* input data matrix pointer */
simon 0:1014af42efd9 57 q31_t *pOut = pDst->pData; /* output data matrix pointer */
simon 0:1014af42efd9 58 q31_t *px; /* Temporary output data matrix pointer */
simon 0:1014af42efd9 59 uint16_t nRows = pSrc->numRows; /* number of nRows */
simon 0:1014af42efd9 60 uint16_t nColumns = pSrc->numCols; /* number of nColumns */
simon 0:1014af42efd9 61 uint16_t blkCnt, i = 0u, row = nRows; /* loop counters */
simon 0:1014af42efd9 62 arm_status status; /* status of matrix transpose */
simon 0:1014af42efd9 63
simon 0:1014af42efd9 64
simon 0:1014af42efd9 65 #ifdef ARM_MATH_MATRIX_CHECK
simon 0:1014af42efd9 66 /* Check for matrix mismatch condition */
simon 0:1014af42efd9 67 if((pSrc->numRows != pDst->numCols) || (pSrc->numCols != pDst->numRows))
simon 0:1014af42efd9 68 {
simon 0:1014af42efd9 69 /* Set status as ARM_MATH_SIZE_MISMATCH */
simon 0:1014af42efd9 70 status = ARM_MATH_SIZE_MISMATCH;
simon 0:1014af42efd9 71 }
simon 0:1014af42efd9 72 else
simon 0:1014af42efd9 73 #endif
simon 0:1014af42efd9 74 {
simon 0:1014af42efd9 75 /* Matrix transpose by exchanging the rows with columns */
simon 0:1014af42efd9 76 /* row loop */
simon 0:1014af42efd9 77 do
simon 0:1014af42efd9 78 {
simon 0:1014af42efd9 79 /* Apply loop unrolling and exchange the columns with row elements */
simon 0:1014af42efd9 80 blkCnt = nColumns >> 2u;
simon 0:1014af42efd9 81
simon 0:1014af42efd9 82 /* The pointer px is set to starting address of the column being processed */
simon 0:1014af42efd9 83 px = pOut + i;
simon 0:1014af42efd9 84
simon 0:1014af42efd9 85 /* First part of the processing with loop unrolling. Compute 4 outputs at a time.
simon 0:1014af42efd9 86 ** a second loop below computes the remaining 1 to 3 samples. */
simon 0:1014af42efd9 87 while(blkCnt > 0u)
simon 0:1014af42efd9 88 {
simon 0:1014af42efd9 89 /* Read and store the input element in the destination */
simon 0:1014af42efd9 90 *px = *pIn++;
simon 0:1014af42efd9 91
simon 0:1014af42efd9 92 /* Update the pointer px to point to the next row of the transposed matrix */
simon 0:1014af42efd9 93 px += nRows;
simon 0:1014af42efd9 94
simon 0:1014af42efd9 95 /* Read and store the input element in the destination */
simon 0:1014af42efd9 96 *px = *pIn++;
simon 0:1014af42efd9 97
simon 0:1014af42efd9 98 /* Update the pointer px to point to the next row of the transposed matrix */
simon 0:1014af42efd9 99 px += nRows;
simon 0:1014af42efd9 100
simon 0:1014af42efd9 101 /* Read and store the input element in the destination */
simon 0:1014af42efd9 102 *px = *pIn++;
simon 0:1014af42efd9 103
simon 0:1014af42efd9 104 /* Update the pointer px to point to the next row of the transposed matrix */
simon 0:1014af42efd9 105 px += nRows;
simon 0:1014af42efd9 106
simon 0:1014af42efd9 107 /* Read and store the input element in the destination */
simon 0:1014af42efd9 108 *px = *pIn++;
simon 0:1014af42efd9 109
simon 0:1014af42efd9 110 /* Update the pointer px to point to the next row of the transposed matrix */
simon 0:1014af42efd9 111 px += nRows;
simon 0:1014af42efd9 112
simon 0:1014af42efd9 113 /* Decrement the column loop counter */
simon 0:1014af42efd9 114 blkCnt--;
simon 0:1014af42efd9 115 }
simon 0:1014af42efd9 116
simon 0:1014af42efd9 117 /* Perform matrix transpose for last 3 samples here. */
simon 0:1014af42efd9 118 blkCnt = nColumns % 0x4u;
simon 0:1014af42efd9 119
simon 0:1014af42efd9 120 while(blkCnt > 0u)
simon 0:1014af42efd9 121 {
simon 0:1014af42efd9 122 /* Read and store the input element in the destination */
simon 0:1014af42efd9 123 *px = *pIn++;
simon 0:1014af42efd9 124
simon 0:1014af42efd9 125 /* Update the pointer px to point to the next row of the transposed matrix */
simon 0:1014af42efd9 126 px += nRows;
simon 0:1014af42efd9 127
simon 0:1014af42efd9 128 /* Decrement the column loop counter */
simon 0:1014af42efd9 129 blkCnt--;
simon 0:1014af42efd9 130 }
simon 0:1014af42efd9 131
simon 0:1014af42efd9 132 i++;
simon 0:1014af42efd9 133
simon 0:1014af42efd9 134 /* Decrement the row loop counter */
simon 0:1014af42efd9 135 row--;
simon 0:1014af42efd9 136
simon 0:1014af42efd9 137 }
simon 0:1014af42efd9 138 while(row > 0u); /* row loop end */
simon 0:1014af42efd9 139
simon 0:1014af42efd9 140 /* set status as ARM_MATH_SUCCESS */
simon 0:1014af42efd9 141 status = ARM_MATH_SUCCESS;
simon 0:1014af42efd9 142 }
simon 0:1014af42efd9 143
simon 0:1014af42efd9 144 /* Return to application */
simon 0:1014af42efd9 145 return (status);
simon 0:1014af42efd9 146 }
simon 0:1014af42efd9 147
simon 0:1014af42efd9 148 /**
simon 0:1014af42efd9 149 * @} end of MatrixTrans group
simon 0:1014af42efd9 150 */