1 //=================================================================================================
2 /*!
3 // \file src/mathtest/operations/smatdmatmult/UCbHDb.cpp
4 // \brief Source file for the UCbHDb sparse matrix/dense matrix multiplication math test
5 //
6 // Copyright (C) 2012-2020 Klaus Iglberger - All Rights Reserved
7 //
8 // This file is part of the Blaze library. You can redistribute it and/or modify it under
9 // the terms of the New (Revised) BSD License. Redistribution and use in source and binary
10 // forms, with or without modification, are permitted provided that the following conditions
11 // are met:
12 //
13 // 1. Redistributions of source code must retain the above copyright notice, this list of
14 // conditions and the following disclaimer.
15 // 2. Redistributions in binary form must reproduce the above copyright notice, this list
16 // of conditions and the following disclaimer in the documentation and/or other materials
17 // provided with the distribution.
18 // 3. Neither the names of the Blaze development group nor the names of its contributors
19 // may be used to endorse or promote products derived from this software without specific
20 // prior written permission.
21 //
22 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY
23 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
24 // OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT
25 // SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
26 // INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED
27 // TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR
28 // BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
29 // CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
30 // ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH
31 // DAMAGE.
32 */
33 //=================================================================================================
34
35
36 //*************************************************************************************************
37 // Includes
38 //*************************************************************************************************
39
40 #include <cstdlib>
41 #include <iostream>
42 #include <blaze/math/CompressedMatrix.h>
43 #include <blaze/math/DynamicMatrix.h>
44 #include <blaze/math/HermitianMatrix.h>
45 #include <blaze/math/UpperMatrix.h>
46 #include <blazetest/mathtest/Creator.h>
47 #include <blazetest/mathtest/operations/smatdmatmult/OperationTest.h>
48 #include <blazetest/system/MathTest.h>
49
50 #ifdef BLAZE_USE_HPX_THREADS
51 # include <hpx/hpx_main.hpp>
52 #endif
53
54
55 //=================================================================================================
56 //
57 // MAIN FUNCTION
58 //
59 //=================================================================================================
60
61 //*************************************************************************************************
main()62 int main()
63 {
64 std::cout << " Running 'UCbHDb'..." << std::endl;
65
66 using blazetest::mathtest::ScalarB;
67
68 try
69 {
70 // Matrix type definitions
71 using UCb = blaze::UpperMatrix< blaze::CompressedMatrix<ScalarB> >;
72 using HDb = blaze::HermitianMatrix< blaze::DynamicMatrix<ScalarB> >;
73
74 // Creator type definitions
75 using CUCb = blazetest::Creator<UCb>;
76 using CHDb = blazetest::Creator<HDb>;
77
78 // Running tests with small matrices
79 for( size_t i=0UL; i<=6UL; ++i ) {
80 for( size_t j=0UL; j<=UCb::maxNonZeros( i ); ++j ) {
81 RUN_SMATDMATMULT_OPERATION_TEST( CUCb( i, j ), CHDb( i ) );
82 }
83 }
84
85 // Running tests with large matrices
86 RUN_SMATDMATMULT_OPERATION_TEST( CUCb( 31UL, 7UL ), CHDb( 31UL ) );
87 RUN_SMATDMATMULT_OPERATION_TEST( CUCb( 67UL, 7UL ), CHDb( 67UL ) );
88 RUN_SMATDMATMULT_OPERATION_TEST( CUCb( 127UL, 13UL ), CHDb( 127UL ) );
89 RUN_SMATDMATMULT_OPERATION_TEST( CUCb( 32UL, 8UL ), CHDb( 32UL ) );
90 RUN_SMATDMATMULT_OPERATION_TEST( CUCb( 64UL, 8UL ), CHDb( 64UL ) );
91 RUN_SMATDMATMULT_OPERATION_TEST( CUCb( 128UL, 16UL ), CHDb( 128UL ) );
92 }
93 catch( std::exception& ex ) {
94 std::cerr << "\n\n ERROR DETECTED during sparse matrix/dense matrix multiplication:\n"
95 << ex.what() << "\n";
96 return EXIT_FAILURE;
97 }
98
99 return EXIT_SUCCESS;
100 }
101 //*************************************************************************************************
102