Function SECTION

tests/unit/channels.cpp:48–94

TEST_CASE( "createKrausMap", TEST_CATEGORY ) {
 
    SECTION( LABEL_CORRECTNESS ) {

        int numQubits = GENERATE( range(1,4) ); // 4qb superop = 8qb matrix = 16 statevec
        int numMatrices = GENERATE( 1, 2, 10 );
        CAPTURE( numQubits, numMatrices );

        KrausMap map = createKrausMap(numQubits, numMatrices);

        // verify dimensions
        REQUIRE( map.numQubits == numQubits );
        REQUIRE( map.numMatrices == numMatrices );
        REQUIRE( map.numRows == getPow2(numQubits) );
        
        // verify superoperator dimensions
        REQUIRE( map.superop.numQubits == numQubits );
        REQUIRE( map.superop.numRows == getPow2(2 * numQubits) );

        // verify default fields
        REQUIRE( *(map.isApproxCPTP) == -1 );
        REQUIRE( *(map.superop.wasGpuSynced) == 0 );

        // verify pointers
        REQUIRE( map.matrices != nullptr );
        REQUIRE( map.superop.cpuElems != nullptr );
        REQUIRE( map.superop.cpuElemsFlat != nullptr );
        if (getQuESTEnv().isGpuAccelerated)
            REQUIRE( map.superop.gpuElemsFlat != nullptr );
        else
        REQUIRE( map.superop.gpuElemsFlat == nullptr );

        // verify that all matrices default to zero
        bool isZero = true;
        for (qindex i=0; i<map.numMatrices && isZero; i++)
            for (qindex r=0; r<map.numRows && isZero; r++)
                for (qindex c=0; c<map.numRows && isZero; c++)
                    isZero = (map.matrices[i][r][c] == qcomp(0,0));
        REQUIRE( isZero );

        // verify superoperator defaults to zero
        isZero = true;
        for (qindex r=0; r<map.superop.numRows && isZero; r++)
            for (qindex c=0; c<map.superop.numRows && isZero; c++)
                isZero = (map.superop.cpuElems[r][c] == qcomp(0,0));
        REQUIRE( isZero );

        destroyKrausMap(map);
    }

    SECTION( LABEL_VALIDATION ) {