examples/cl_convolution.cpp

   1 /*
   2  * Copyright (c) 2016, 2017 ARM Limited.
   3  *
   4  * SPDX-License-Identifier: MIT
   5  *
   6  * Permission is hereby granted, free of charge, to any person obtaining a copy
   7  * of this software and associated documentation files (the "Software"), to
   8  * deal in the Software without restriction, including without limitation the
   9  * rights to use, copy, modify, merge, publish, distribute, sublicense, and/or
  10  * sell copies of the Software, and to permit persons to whom the Software is
  11  * furnished to do so, subject to the following conditions:
  12  *
  13  * The above copyright notice and this permission notice shall be included in all
  14  * copies or substantial portions of the Software.
  15  *
  16  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  17  * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  18  * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
  19  * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  20  * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
  21  * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
  22  * SOFTWARE.
  23  */
  24 #ifndef ARM_COMPUTE_CL /* Needed by Utils.cpp to handle OpenCL exceptions properly */
  25 #error "This example needs to be built with -DARM_COMPUTE_CL"
  26 #endif /* ARM_COMPUTE_CL */
  27
  28 #include "arm_compute/core/Types.h"
  29 #include "arm_compute/runtime/CL/CLFunctions.h"
  30 #include "arm_compute/runtime/CL/CLScheduler.h"
  31 #include "utils/Utils.h"
  32
  33 using namespace arm_compute;
  34 using namespace utils;
  35
  36 /** Gaussian 3x3 matrix
  37  */
  38 const int16_t gaussian3x3[] =
  39 {
  40     1, 2, 1,
  41     2, 4, 2,
  42     1, 2, 1
  43 };
  44
  45 /** Gaussian 5x5 matrix
  46  */
  47 const int16_t gaussian5x5[] =
  48 {
  49     1, 4, 6, 4, 1,
  50     4, 16, 24, 16, 4,
  51     6, 24, 36, 24, 6,
  52     4, 16, 24, 16, 4,
  53     1, 4, 6, 4, 1
  54 };
  55
  56 void main_cl_convolution(int argc, const char **argv)
  57 {
  58     PPMLoader ppm;
  59     CLImage   src, tmp, dst;
  60
  61     CLScheduler::get().default_init();
  62
  63     if(argc < 2)
  64     {
  65         // Print help
  66         std::cout << "Usage: ./build/cl_convolution [input_image.ppm]\n\n";
  67         std::cout << "No input_image provided, creating a dummy 640x480 image\n";
  68         // Create an empty grayscale 640x480 image
  69         src.allocator()->init(TensorInfo(640, 480, Format::U8));
  70     }
  71     else
  72     {
  73         ppm.open(argv[1]);
  74         ppm.init_image(src, Format::U8);
  75     }
  76
  77     // Configure the temporary and destination images
  78     tmp.allocator()->init(*src.info());
  79     dst.allocator()->init(*src.info());
  80
  81     CLConvolution3x3 conv3x3;
  82     CLConvolution5x5 conv5x5;
  83
  84     // Apply a Gaussian 3x3 filter to the source image followed by a Gaussian 5x5:
  85     conv3x3.configure(&src, &tmp, gaussian3x3, 0 /* Let arm_compute calculate the scale */, BorderMode::UNDEFINED);
  86     conv5x5.configure(&tmp, &dst, gaussian5x5, 0 /* Let arm_compute calculate the scale */, BorderMode::UNDEFINED);
  87
  88     // Allocate all the images
  89     src.allocator()->allocate();
  90     tmp.allocator()->allocate();
  91     dst.allocator()->allocate();
  92     // Fill the input image with the content of the PPM image if a filename was provided:
  93     if(ppm.is_open())
  94     {
  95         ppm.fill_image(src);
  96     }
  97
  98     // Execute the functions:
  99     conv3x3.run();
 100     conv5x5.run();
 101
 102     // Make sure all the OpenCL jobs are done executing:
 103     CLScheduler::get().sync();
 104
 105     // Save the result to file:
 106     if(ppm.is_open())
 107     {
 108         const std::string output_filename = std::string(argv[1]) + "_out.ppm";
 109         save_to_ppm(dst, output_filename); // save_to_ppm maps and unmaps the image to store as PPM
 110     }
 111 }
 112
 113 /** Main program for convolution test
 114  *
 115  * @param[in] argc Number of arguments
 116  * @param[in] argv Arguments ( [optional] Path to PPM image to process )
 117  */
 118 int main(int argc, const char **argv)
 119 {
 120     return utils::run_example(argc, argv, main_cl_convolution);
 121 }