cs_array_2dspan.h

Go to the documentation of this file.

#ifndef __CS_ARRAY_2DSPAN_H__
#define __CS_ARRAY_2DSPAN_H__
 
/*============================================================================
 * Templated 2D array object
 *============================================================================*/
 
/*
  This file is part of code_saturne, a general-purpose CFD tool.
 
  Copyright (C) 1998-2025 EDF S.A.
 
  This program is free software; you can redistribute it and/or modify it under
  the terms of the GNU General Public License as published by the Free Software
  Foundation; either version 2 of the License, or (at your option) any later
  version.
 
  This program is distributed in the hope that it will be useful, but WITHOUT
  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
  details.
 
  You should have received a copy of the GNU General Public License along with
  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
  Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
 
/*----------------------------------------------------------------------------*/
 
/*----------------------------------------------------------------------------
 *  Local headers
 *----------------------------------------------------------------------------*/
 
#include "bft/bft_error.h"
 
#include "base/cs_defs.h"
 
#include "base/cs_array.h"
#include "base/cs_mem.h"
 
#if defined(__cplusplus)
 
/*=============================================================================
 * BUILTIN TEST
 *============================================================================*/
 
#ifndef __has_builtin
#define __has_builtin(x) 0
#endif
 
/*=============================================================================
 * Public C++ class template
 *============================================================================*/
 
namespace cs {
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
template <class T>
class array_2dspan {
public:
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  array_2dspan() :
    _dim1(0),
    _dim2(0),
    _size(0),
    _is_owner(true),
    _full_array(nullptr),
    _mode(cs_alloc_mode)
  {
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  array_2dspan
  (
    cs_lnum_t dim1, 
    cs_lnum_t dim2, 
#if (defined(__GNUC__) || defined(__clang__)) && \
  __has_builtin(__builtin_LINE) && \
 __has_builtin(__builtin_FILE)
    const char *file_name   = __builtin_FILE(), 
    const int   line_number = __builtin_LINE()  
#else
    const char *file_name   = __FILE__, 
    const int   line_number = __LINE__  
#endif
  )
  :
    _dim1(dim1),
    _dim2(dim2),
    _size(dim1*dim2),
    _is_owner(true),
    _mode(cs_alloc_mode)
  {
    allocate_(file_name, line_number);
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  array_2dspan
  (
    cs_lnum_t       dim1,       
    cs_lnum_t       dim2,       
    cs_alloc_mode_t alloc_mode, 
#if (defined(__GNUC__) || defined(__clang__)) && \
  __has_builtin(__builtin_LINE) && \
  __has_builtin(__builtin_FILE)
    const char *file_name   = __builtin_FILE(), 
    const int   line_number = __builtin_LINE()  
#else
    const char *file_name   = __FILE__, 
    const int   line_number = __LINE__  
#endif
  )
  :
    _dim1(dim1),
    _dim2(dim2),
    _size(dim1 * dim2),
    _is_owner(true),
    _mode(alloc_mode)
  {
    allocate_(file_name, line_number);
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  array_2dspan
  (
    cs_lnum_t       dim1,                      
    cs_lnum_t       dim2,                      
    T*              data_array,                
    cs_alloc_mode_t alloc_mode = cs_alloc_mode 
  )
  :
    _dim1(dim1),
    _dim2(dim2),
    _size(dim1 * dim2),
    _is_owner(false),
    _mode(alloc_mode),
    _full_array(data_array)
  {
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  array_2dspan
  (
    array_2dspan& other,              
    bool          shallow_copy=false, 
#if (defined(__GNUC__) || defined(__clang__)) && \
  __has_builtin(__builtin_LINE) && \
  __has_builtin(__builtin_FILE)
    const char *file_name   = __builtin_FILE(), 
    const int   line_number = __builtin_LINE()  
#else
    const char *file_name   = __FILE__, 
    const int   line_number = __LINE__  
#endif
  )
  {
    set_size_(other._dim1, other._dim2);
    _mode = other._mode;
 
    /* If shallow copy new instance is not owner. Otherwise same ownership
     * as original instance since we copy it.
     */
    _is_owner = (shallow_copy) ? false : other._is_owner;
 
    if (_is_owner) {
      allocate_(file_name, line_number);
      cs_array_copy<T>(_size, other._full_array, _full_array);
    }
    else {
      _full_array = other._full_array;
    }
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  array_2dspan
  (
    array_2dspan&& other 
  )
    : array_2dspan()
  {
    swap(*this, other);
  }
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  ~array_2dspan()
  {
    clear();
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  friend void
  swap
  (
    array_2dspan& first, 
    array_2dspan& second 
  )
  {
    using std::swap;
    /* Swap the different members between the two references. */
    swap(first._dim1, second._dim1);
    swap(first._dim2, second._dim2);
    swap(first._size, second._size);
    swap(first._is_owner, second._is_owner);
    swap(first._mode, second._mode);
    swap(first._full_array, second._full_array);
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  array_2dspan& operator=(array_2dspan other)
  {
    swap(*this, other);
 
    return *this;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void
  clear()
  {
    if (_is_owner) {
      CS_FREE(_full_array);
    }
    else {
      _full_array = nullptr;
    }
 
    set_size_(0,0);
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void
  empty()
  {
    set_size_(0, 0);
    _is_owner = false;
    _full_array = nullptr;
    _mode = cs_alloc_mode;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void
  point_to
  (
    array_2dspan& other 
  )
  {
    clear();
    *this = other.view();
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void set_to_val
  (
    T val 
  )
  {
    cs_arrays_set_value<T,1>(_size, val, _full_array);
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void resize
  (
    cs_lnum_t       dim1, 
    cs_lnum_t       dim2, 
#if (defined(__GNUC__) || defined(__clang__)) && \
  __has_builtin(__builtin_LINE) && \
  __has_builtin(__builtin_FILE)
    const char *file_name   = __builtin_FILE(), 
    const int   line_number = __builtin_LINE()  
#else
    const char *file_name   = __FILE__, 
    const int   line_number = __LINE__  
#endif
  )
  {
    assert(dim1 >= 0 && dim2 >= 0);
 
    /* If same dimensions, nothing to do ... */
    if (dim1 == _dim1 && dim2 == _dim2)
      return;
 
    if (_is_owner) {
      clear();
      set_size_(dim1, dim2);
      allocate_(file_name, line_number);
    }
 
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void resize
  (
    cs_lnum_t       dim1,         
    cs_lnum_t       dim2,         
    bool            copy_data,    
    cs_lnum_t       size_to_keep, 
#if (defined(__GNUC__) || defined(__clang__)) && \
  __has_builtin(__builtin_LINE) && \
  __has_builtin(__builtin_FILE)
    const char *file_name   = __builtin_FILE(), 
    const int   line_number = __builtin_LINE()  
#else
    const char *file_name   = __FILE__, 
    const int   line_number = __LINE__  
#endif
  )
  {
    assert(dim1 >= 0 && dim2 >= 0);
    assert(size_to_keep <= dim2 && size_to_keep <= _dim2);
 
    /* If same dimensions, nothing to do ... */
    if (dim1 == _dim1 && dim2 == _dim2)
      return;
 
    if (copy_data && dim1 < _dim1)
      bft_error(__FILE__, __LINE__, 0,
                "%s: Data cannot be saved when new dim1 is smaller than previous.\n",
                __func__);
 
    if (_is_owner) {
      if (copy_data) {
        /* If we change dim1 -> Realloc is sufficient */
        if (_dim1 != dim1) {
          set_size_(dim1, dim2);
          reallocate_(file_name, line_number);
        }
        else {
          /* Temporary copy */
          T *tmp = nullptr;
          CS_MALLOC_HD(tmp, _size, T, _mode);
 
          cs_array_copy<T>(_size, _full_array, tmp);
 
          cs_lnum_t old_dim2 = _dim2;
 
          resize(dim1, dim2, file_name, line_number);
 
          /* We loop on "_dim1" since dim1 = _dim1 */
          for (cs_lnum_t i = 0; i < _dim1; i++) {
            cs_array_copy<T>(size_to_keep,
                             tmp + i*old_dim2,
                             _full_array + i*_dim2);
          }
 
          CS_FREE(tmp);
        }
      }
      else {
        resize(dim1, dim2, file_name, line_number);
      }
    }
 
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  array_2dspan view()
  {
    /* Use the shallow copy constructor to return a temporary instance */
    return array_2dspan(*this, true);
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void set_alloc_mode
  (
    cs_alloc_mode_t mode 
  )
  {
    if (_mode != mode) {
      if (_size != 0)
        clear();
 
      _mode = mode;
    }
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  T *vals()
  {
    return _full_array;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  const T *vals() const
  {
    return _full_array;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  T* operator[]
  (
    int i 
  )
  {
    return _full_array + i * _dim2;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  const T* operator[]
  (
    int i 
  ) const
  {
    return _full_array + i * _dim2;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  cs_lnum_t dim1()
  {
    return _dim1;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  cs_lnum_t dim2()
  {
    return _dim2;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  cs_lnum_t size()
  {
    return _size;
  }
 
private:
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void
  set_size_
  (
    cs_lnum_t dim1,
    cs_lnum_t dim2
  )
  {
    _dim1 = dim1;
    _dim2 = dim2;
    _size = dim1 * dim2;
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void
  allocate_
  (
    const char *file_name,  
    const int   line_number 
  )
  {
    const char *_ptr_name = "_full_array";
    _full_array = static_cast<T *>(cs_mem_malloc_hd(_mode,
                                                    _size,
                                                    sizeof(T),
                                                    _ptr_name,
                                                    file_name,
                                                    line_number));
  }
 
  /*--------------------------------------------------------------------------*/
  /*--------------------------------------------------------------------------*/
 
  CS_F_HOST_DEVICE
  void
  reallocate_
  (
    const char *file_name,  
    const int   line_number 
  )
  {
    /* If not owner no-op */
    if (_is_owner) {
      const char *_ptr_name = "_full_array";
      _full_array = static_cast<T *>(cs_mem_realloc_hd(_full_array,
                                                       _mode,
                                                       _size,
                                                       sizeof(T),
                                                       _ptr_name,
                                                       file_name,
                                                       line_number));
    }
  };
 
  /*--------------------------------------------------------------------------*/
  /* Private members */
  /*--------------------------------------------------------------------------*/
 
  cs_lnum_t       _dim1;       
  cs_lnum_t       _dim2;       
  cs_lnum_t       _size;       
  bool            _is_owner;   
  T*              _full_array; 
  cs_alloc_mode_t _mode;       
};
 
}
 
/*----------------------------------------------------------------------------*/
 
#endif // __cplusplus
 
/*----------------------------------------------------------------------------*/
 
#endif // __CS_ARRAY_2DSPAN_H__