File: itkOpenCLBufferTest.cxx

package info (click to toggle)
elastix 5.2.0-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 42,480 kB
  • sloc: cpp: 68,403; lisp: 4,118; python: 1,013; xml: 182; sh: 177; makefile: 33
file content (116 lines) | stat: -rw-r--r-- 3,667 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
/*=========================================================================
 *
 *  Copyright UMC Utrecht and contributors
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *        http://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/
#include "itkTestHelper.h"
#include "itkOpenCLBufferTest.h"
#include "itkOpenCLBuffer.h"
#include <algorithm>

//------------------------------------------------------------------------------
template <class type>
bool
std_all_of(const std::vector<type> & v, const type value)
{
  for (typename std::vector<type>::const_iterator it = v.begin(); it != v.end(); ++it)
  {
    if (*it != value)
    {
      return false;
    }
  }
  return true;
}


//------------------------------------------------------------------------------
int
main()
{
  itk::OpenCLBuffer bufferNull;

  if (!bufferNull.IsNull())
  {
    return EXIT_FAILURE;
  }

  try
  {
    itk::OpenCLContext::Pointer context = itk::OpenCLContext::GetInstance();
    context->Create(itk::OpenCLContext::DevelopmentSingleMaximumFlopsDevice);
    const std::list<itk::OpenCLDevice> devices = context->GetDevices();

    itk::OpenCLProgram program =
      context->BuildProgramFromSourceCode(devices, itk::OpenCLBufferTestKernel::GetOpenCLSource());

    if (program.IsNull())
    {
      if (context->GetDefaultDevice().HasCompiler())
      {
        itkGenericExceptionMacro("Could not compile the OpenCL test program");
      }
      else
      {
        itkGenericExceptionMacro("OpenCL implementation does not have a compiler");
      }
    }

    // Local variables
    const std::size_t     bufferSize = 16;
    const std::size_t     bufferSizeBytes = sizeof(float) * bufferSize;
    const float           value = 5.1f;
    const itk::OpenCLSize workSize(bufferSize);

    // Tests the OpenCLBuffer
    itk::OpenCLBuffer deviceBuffer = context->CreateBufferDevice(itk::OpenCLMemoryObject::WriteOnly, bufferSizeBytes);
    ITK_OPENCL_COMPARE(deviceBuffer.IsNull(), false);
    std::cout << deviceBuffer << std::endl;
    itk::OpenCLKernel setFloatKernel = program.CreateKernel("SetFloat");
    ITK_OPENCL_COMPARE(setFloatKernel.IsNull(), false);
    setFloatKernel(deviceBuffer, value);

    // Compare results test1
    std::vector<float> hostBuffer(bufferSize);
    deviceBuffer.Read(&hostBuffer[0], sizeof(float));
    ITK_OPENCL_COMPARE(hostBuffer[0], value);
    ITK_OPENCL_COMPARE(hostBuffer[1], 0.0f);

    // Compare results test2
    setFloatKernel.SetGlobalWorkSize(workSize);
    itk::OpenCLEvent event = setFloatKernel.LaunchKernel();
    ITK_OPENCL_COMPARE(event.IsNull(), false);
    event.WaitForFinished();

    deviceBuffer.Read(&hostBuffer[0], bufferSizeBytes);

    const bool result = std_all_of<float>(hostBuffer, value);

    // Only in the latest compilers
    // const bool result = std::all_of( hostBuffer.begin(), hostBuffer.end(),
    //  [ value ]( float i ){
    //    return i == value;
    //  } );

    ITK_OPENCL_COMPARE(result, true);
  }
  catch (itk::ExceptionObject & e)
  {
    std::cerr << "Caught ITK exception: " << e << std::endl;
    return EXIT_FAILURE;
  }

  return EXIT_SUCCESS;
}