bwt.c

这是在PCA下的基于IPP库示例代码例子,在网上下了IPP的库之后,设置相关参数就可以编译该代码.

/*//////////////////////////////// "bwt.c" /////////////////////////////
//
//                  INTEL CORPORATION PROPRIETARY INFORMATION
//     This software is supplied under the terms of a license agreement or
//     nondisclosure agreement with Intel Corporation and may not be copied
//     or disclosed except in accordance with the terms of that agreement.
//          Copyright(c) 2005 Intel Corporation. All Rights Reserved.
//
//     higher level functions of IPP DC ( data compression ) domain for
//     BWT forward/backward transform
*/

#include "bwt.h"

/*************************************************************************************************
 * Function:
 * int ForwardBWT( Ipp8u * pSrc, int SrcLen, Ipp8u *pDst, int *pDstLen, Ipp8u *addBuff )
 *
 * Function for relization of forward BWT
 * Gets the pSrc as a source vector of length SrcLen, does BWT and writes result to pDst
 * ( and length of pDst in pDstLen )
 *
 * Main Parameters:
 *          pSrc         - pointer to Source buffer,
 *          pDst         - pointer to Destination buffer,
 *          SrcLen       - Length of source buffer,
 *          pDstLen      - Length of destination buffer
 *          addBuff      - Additional buffer for BWT transform
 *
 * Returns:
 *    ippStsNullPtrErr  If One or several pointer(s) is 0
 *    ippStsSizeErr     Length of the source vector is less or equal zero
 *    ippStsNoErr       No errors
 *
 ************************************************************************************************/
int ForwardBWT( Ipp8u *pSrc, int SrcLen, Ipp8u *pDst, int *pDstLen, Ipp8u *addBuff ) {

    IppStatus st;           /* to store and analize return statuses of ippdc function */
    int firstMeetAt = 0;    /* to store the first occurence index */

    /* checks the pointer is not null */
    if( pSrc == 0 || pDst == 0 )
        return ippStsNullPtrErr;

    /* checks the pointer is not null */
    if( addBuff == 0 )
        return ippStsNullPtrErr;

    /* size of  vector after bwt is 4 bytes longer, due to the
       first 4 bytes contain the 'first occurence index' value */
    *pDstLen = SrcLen + sizeof(int);

    /* does the forward BWT using ippdc function */
    if( ippStsNoErr != ( st = ippsBWTFwd_8u( pSrc, pDst + sizeof(int), SrcLen, &firstMeetAt, addBuff ) ) ) {
        fprintf( stderr, "Error <%d> while trying forward BWT. Exiting.\n", st );
        return st;
    }

    /* store the first occurence index in first 4 bytes of pSrc vector */
    ((int*)pDst)[0] = firstMeetAt;

    return 0; /* No errors */
}

/*************************************************************************************************
 * Function:
 * int InverseBWT( Ipp8u * pSrc, int SrcLen, Ipp8u *pDst, int *pDstLen, Ipp8u *addBuff )
 *
 * Function for relization of backward BWT
 * Gets the pSrc as a source vector of length SrcLen, does BWT and writes result to pDst
 * ( and length of pDst in pDstLen )
 *
 * Main Parameters:
 *          pSrc         - pointer to Source buffer,
 *          pDst         - pointer to Destination buffer,
 *          SrcLen       - Length of source buffer,
 *          pDstLen      - Length of destination buffer
 *          addBuff      - Additional buffer for BWT transform
 *
 * Returns:
 *    ippStsNullPtrErr  If One or several pointer(s) is 0
 *    ippStsSizeErr     Length of the source vector is less or equal zero
 *    ippStsNoErr       No errors
 *
 ************************************************************************************************/
int InverseBWT( Ipp8u * pSrc, int SrcLen, Ipp8u *pDst, int *pDstLen, Ipp8u *addBuff ) {

    IppStatus st;           /* to store and analize return statuses of ippdc function */
    int bwtFirstShot = 0;   /* to store the first occurence index */

    /* checks the pointer is not null */
    if( addBuff == 0)
        return ippStsNullPtrErr;

    /* checks the pointer is not null */
    if( pSrc == 0 || pDst == 0 )
        return ippStsNullPtrErr;

    /* reads first four bytes to get the first occurence index for BWT */
    bwtFirstShot = ((int*)pSrc)[0];

    /* resulting length is 4 bytes shorter */
    *pDstLen = SrcLen - sizeof(int);

    /* does the backward BWT using ippdc function */
    if( ippStsNoErr != ( st = ippsBWTInv_8u( pSrc + ( sizeof(int) ), pDst, SrcLen - sizeof(int), bwtFirstShot, addBuff ) ) ) {
        fprintf(stderr, "Error <%d> while trying do BWT. Exiting.\n", st);
        return st;
    }

    return 0; /* No errors */
}