C++ (Cpp) DenseVector::l2Norm Examples

Example #1

File: P_NormTest.cpp Project: fast-project/lama

BOOST_AUTO_TEST_CASE_TEMPLATE( Norm, T, test_types ) 
{
    typedef T ValueType;

    CONTEXTLOOP()
    {
        GETCONTEXT( context );

        for ( IndexType size = 0; size < 4; ++size )
        {
            DistributionPtr dist( new BlockDistribution( size, comm ) );

            const ValueType VAL = 1.0;

            DenseVector<ValueType> repVector( size, VAL );

            DenseVector<ValueType> distVector;

            distVector.setContext( context );

            distVector = repVector;

            distVector.redistribute( dist );
    
            Scalar l1norm = distVector.l1Norm();

            ValueType expectedL1Norm = static_cast<ValueType>( size );

            BOOST_CHECK_CLOSE( l1norm.getValue<ValueType>(), expectedL1Norm, 1 );
    
            Scalar l2norm = distVector.l2Norm();

            ValueType expectedL2Norm = static_cast<ValueType>( size );
            expectedL2Norm = std::sqrt( expectedL2Norm );
    
            BOOST_CHECK_CLOSE( l2norm.getValue<ValueType>(), expectedL2Norm, 1 );
    
            Scalar maxNorm = distVector.maxNorm();

            ValueType expectedMaxNorm = static_cast<ValueType>( VAL );

            if ( size == 0 )
            {
                expectedMaxNorm = static_cast<ValueType>( 0 );
            }

            BOOST_CHECK_CLOSE( maxNorm.getValue<ValueType>(), expectedMaxNorm, 1 );
        }
    }
}