truncate.h

Go to the documentation of this file.

/*
 *  This file is a part of TiledArray.
 *  Copyright (C) 2015  Virginia Tech
 *
 *  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 3 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, see <http://www.gnu.org/licenses/>.
 *
 *  Justus Calvin
 *  Department of Chemistry, Virginia Tech
 *
 *  truncate.h
 *  Apr 15, 2015
 *
 */
 
#ifndef TILEDARRAY_CONVERSIONS_TRUNCATE_H__INCLUDED
#define TILEDARRAY_CONVERSIONS_TRUNCATE_H__INCLUDED
 
#include <TiledArray/conversions/foreach.h>
 
namespace TiledArray {
 
template <typename, typename>
class DistArray;
class DensePolicy;
class SparsePolicy;
 
 
template <typename Tile, typename Policy>
inline std::enable_if_t<is_dense_v<Policy>, void> truncate(
    DistArray<Tile, Policy>& array,
    typename Policy::shape_type::value_type = 0) {}
 
 
template <typename Tile, typename Policy>
inline std::enable_if_t<!is_dense_v<Policy>, void> truncate(
    DistArray<Tile, Policy>& array,
    typename Policy::shape_type::value_type thresh =
        Policy::shape_type::threshold()) {
  TA_ASSERT(thresh >= 0);
  const typename Policy::shape_type::value_type previous_thresh =
      Policy::shape_type::threshold();
  const auto need_to_change_thresh = (thresh != previous_thresh);
  if (need_to_change_thresh)
    array.world().gop.serial_invoke(
        [thresh] { Policy::shape_type::threshold(thresh); });
  typedef typename DistArray<Tile, Policy>::value_type value_type;
  array = foreach (array,
                   [](value_type& result_tile, const value_type& arg_tile) ->
                   typename Policy::shape_type::value_type {
                     using result_type =
                         typename Policy::shape_type::value_type;
                     result_type arg_tile_norm;
                     norm(arg_tile, arg_tile_norm);
                     result_tile = arg_tile;  // Assume this is shallow copy
                     return arg_tile_norm;
                   });
  if (need_to_change_thresh)
    array.world().gop.serial_invoke(
        [previous_thresh] { Policy::shape_type::threshold(previous_thresh); });
}
 
}  // namespace TiledArray
 
#endif  // TILEDARRAY_CONVERSIONS_TRUNCATE_H__INCLUDED