Skip to content

Commit

Permalink
reset 2x2 file
Browse files Browse the repository at this point in the history
  • Loading branch information
blegouix committed Jun 26, 2024
1 parent 2597448 commit 93a2b2d
Showing 1 changed file with 2 additions and 30 deletions.
32 changes: 2 additions & 30 deletions include/ddc/kernels/splines/splines_linear_problem_2x2_blocks.hpp
Original file line number Diff line number Diff line change
Expand Up @@ -50,7 +50,7 @@ class SplinesLinearProblem2x2Blocks : public SplinesLinearProblem<ExecSpace>
* @brief SplinesLinearProblem2x2Blocks constructor.
*
* @param mat_size The size of one of the dimensions of the square matrix.
* @param top_left_block A pointer toward the top-left SplinesLinearProblem. `setup_solver` must not have been called on it.
* @param top_left_block A pointer toward the top-left SplinesLinearProblem. `setup_solver` must not have been called on `q`.
*/
explicit SplinesLinearProblem2x2Blocks(
std::size_t const mat_size,
Expand All @@ -74,34 +74,6 @@ class SplinesLinearProblem2x2Blocks : public SplinesLinearProblem<ExecSpace>
Kokkos::deep_copy(m_bottom_left_block.h_view, 0.);
}

protected:
/**
* @brief SplinesLinearProblem2x2Blocks constructor.
*
* @param mat_size The size of one of the dimensions of the square matrix.
* @param q A pointer toward the top-left SplinesLinearProblem.
*/
explicit SplinesLinearProblem2x2Blocks(
std::size_t const mat_size,
std::unique_ptr<SplinesLinearProblem<ExecSpace>> top_left_block,
std::size_t const lambda_size1,
std::size_t const lambda_size2)
: SplinesLinearProblem<ExecSpace>(mat_size)
, m_top_left_block(std::move(top_left_block))
, m_top_right_block(
"top_right_block",
m_top_left_block->size(),
mat_size - m_top_left_block->size())
, m_bottom_left_block("bottom_left_block", lambda_size1, lambda_size2)
, m_bottom_right_block(
new SplinesLinearProblemDense<ExecSpace>(mat_size - m_top_left_block->size()))
{
assert(m_top_left_block->size() <= mat_size);

Kokkos::deep_copy(m_top_right_block.h_view, 0.);
Kokkos::deep_copy(m_bottom_left_block.h_view, 0.);
}

double get_element(std::size_t const i, std::size_t const j) const override
{
assert(i < size());
Expand Down Expand Up @@ -201,7 +173,7 @@ class SplinesLinearProblem2x2Blocks : public SplinesLinearProblem<ExecSpace>
* @param LinOp
* @param transpose
*/
virtual void gemv_minus1_1(
void gemv_minus1_1(
MultiRHS const x,
MultiRHS const y,
Kokkos::View<double**, Kokkos::LayoutRight, typename ExecSpace::memory_space> const
Expand Down

0 comments on commit 93a2b2d

Please sign in to comment.