Commit f40ec4cd authored by COULAUD Olivier's avatar COULAUD Olivier
Browse files

Particle re-distribution to balance cells or particles works.

parent 7b9b7677
......@@ -65,13 +65,12 @@ namespace scalfmm::interpolation
static constexpr std::size_t kn = matrix_kernel_type::kn;
static constexpr std::size_t km = matrix_kernel_type::km;
using base_type = interpolator<chebyshev_interpolator<value_type, dimension, matrix_kernel_type>>;
using complex_type = std::complex<value_type>;
using k_tensor_type = xt::xarray<complex_type>;
using k_tensor_type = xt::xarray<value_type>;
using interaction_matrix_type = xt::xtensor_fixed<k_tensor_type, xt::xshape<kn, km>>;
using base_type::base_type;
using buffer_inner_type = xt::xarray<std::complex<value_type>>;
using buffer_shape_type = xt::xshape<kn>;
using buffer_type = xt::xtensor_fixed<buffer_inner_type, buffer_shape_type>;
// using buffer_inner_type = xt::xarray<std::complex<value_type>>;
// using buffer_shape_type = xt::xshape<kn>;
// using buffer_type = xt::xtensor_fixed<buffer_inner_type, buffer_shape_type>;
chebyshev_interpolator() = default;
chebyshev_interpolator(chebyshev_interpolator const&) = default;
......@@ -215,6 +214,9 @@ namespace scalfmm::interpolation
// },
// target_expansion);
}
private:
/// Store Tn(x_i)
};
template<typename ValueType, std::size_t Dimension, typename FarFieldMatrixKernel>
......@@ -224,12 +226,12 @@ namespace scalfmm::interpolation
static constexpr std::size_t dimension = Dimension;
using far_field_matrix_kernel_type = FarFieldMatrixKernel;
using base_type = interpolator<chebyshev_interpolator<value_type, dimension, far_field_matrix_kernel_type>>;
using k_tensor_type = xt::xarray<std::complex<value_type>>;
using buffer_inner_type = xt::xarray<std::complex<value_type>>;
using k_tensor_type = xt::xarray<svalue_type>;
// using buffer_inner_type = xt::xarray<std::complex<value_type>>;
using interaction_matrix_type =
xt::xtensor_fixed<k_tensor_type,
xt::xshape<far_field_matrix_kernel_type::kn, far_field_matrix_kernel_type::km>>;
using buffer_type = xt::xtensor_fixed<buffer_inner_type, xt::xshape<far_field_matrix_kernel_type::kn>>;
// using buffer_type = xt::xtensor_fixed<buffer_inner_type, xt::xshape<far_field_matrix_kernel_type::kn>>;
};
} // namespace scalfmm::interpolation
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment