18380674863

Committed 09 Oct 2025 02:45PM UTC coverage: 75.202% (+0.3%) from 74.862%

Build # 18380674863

Build Type

push

github

Committed by

web-flow

Commit Message

Merge pull request #793 from wdhawkins/remove_cmake_debug_flags

Removing forced use of -O0 for Debug config

Coverage Stats

18014 of 23954 relevant lines covered (75.2%)

56527302.01 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

89.82

/modules/linear_boltzmann_solvers/lbs_problem/lbs_vecops.cc

// SPDX-FileCopyrightText: 2024 The OpenSn Authors <https://open-sn.github.io/opensn/>
// SPDX-License-Identifier: MIT

#include "modules/linear_boltzmann_solvers/lbs_problem/lbs_vecops.h"
#include "modules/linear_boltzmann_solvers/lbs_problem/lbs_problem.h"

namespace opensn
{

template <typename Functor>
int
LBSVecOps::GroupsetScopedCopy(LBSProblem& lbs_problem, int gsi, int gss, Functor func)
{
  CALI_CXX_MARK_SCOPE("LBSVecOps::GroupsetScopedCopy");

  const auto& grid = lbs_problem.GetGrid();
  const auto& cell_transport_views = lbs_problem.GetCellTransportViews();
  auto num_moments = lbs_problem.GetNumMoments();

  int64_t idx = -1;
  for (const auto& cell : grid->local_cells)
  {
    const auto& transport_view = cell_transport_views[cell.local_id];
    for (int i = 0; i < cell.vertex_ids.size(); ++i)
    {
      for (size_t m = 0; m < num_moments; ++m)
      {
        auto mapped_idx = transport_view.MapDOF(i, m, gsi);
        for (int g = 0; g < gss; ++g)
        {
          ++idx;
          func(idx, mapped_idx + g);
        }
      }
    }
  }
  return idx;
}

void
LBSVecOps::SetPhiVectorScalarValues(LBSProblem& lbs_problem, PhiSTLOption phi_opt, double value)
{
  CALI_CXX_MARK_SCOPE("LBSVecOps::SetPhiVectorScalarValues");

  auto& phi = (phi_opt == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal()
                                                 : lbs_problem.GetPhiOldLocal();
  const auto& grid = lbs_problem.GetGrid();
  const auto& groups = lbs_problem.GetGroups();
  const auto& sdm = lbs_problem.GetSpatialDiscretization();
  const auto& unknown_manager = lbs_problem.GetUnknownManager();

  const size_t first_grp = groups.front().id;
  const size_t final_grp = groups.back().id;

  for (const auto& cell : grid->local_cells)
  {
    const auto& cell_mapping = sdm.GetCellMapping(cell);
    const size_t num_nodes = cell_mapping.GetNumNodes();

    for (size_t i = 0; i < num_nodes; ++i)
    {
      const auto dof_map = sdm.MapDOFLocal(cell, i, unknown_manager, 0, 0);
      std::fill(phi.begin() + dof_map + first_grp, phi.begin() + dof_map + final_grp + 1, value);
    }
  }
}

void
LBSVecOps::ScalePhiVector(LBSProblem& lbs_problem, PhiSTLOption phi_opt, double value)
{
  CALI_CXX_MARK_SCOPE("LBSVecOps::ScalePhiVector");

  auto& phi = (phi_opt == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal()
                                                 : lbs_problem.GetPhiOldLocal();
  auto& groupsets = lbs_problem.GetGroupsets();

  Scale(phi, value);
  for (auto& groupset : groupsets)
  {
    if (groupset.angle_agg)
    {
      if (phi_opt == PhiSTLOption::PHI_NEW)
      {
        auto psi = groupset.angle_agg->GetNewDelayedAngularDOFsAsSTLVector();
        Scale(psi, value);
        groupset.angle_agg->SetNewDelayedAngularDOFsFromSTLVector(psi);
      }
      else if (phi_opt == PhiSTLOption::PHI_OLD)
      {
        auto psi = groupset.angle_agg->GetOldDelayedAngularDOFsAsSTLVector();
        Scale(psi, value);
        groupset.angle_agg->SetOldDelayedAngularDOFsFromSTLVector(psi);
      }
    }
  }
}

void
LBSVecOps::SetGSPETScVecFromPrimarySTLvector(LBSProblem& lbs_problem,
                                             const LBSGroupset& groupset,
                                             Vec dest,
                                             PhiSTLOption src)
{
  const auto& src_phi =
    (src == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();
  double* petsc_dest;
  VecGetArray(dest, &petsc_dest);
  int64_t index = GroupsetScopedCopy(lbs_problem,
                                     groupset.groups.front().id,
                                     groupset.groups.size(),
                                     [&](int64_t idx, size_t mapped_idx)
                                     { petsc_dest[idx] = src_phi[mapped_idx]; });
  if (groupset.angle_agg)
  {
    if (src == PhiSTLOption::PHI_NEW)
      groupset.angle_agg->AppendNewDelayedAngularDOFsToArray(index, petsc_dest);
    else if (src == PhiSTLOption::PHI_OLD)
      groupset.angle_agg->AppendOldDelayedAngularDOFsToArray(index, petsc_dest);
  }
  VecRestoreArray(dest, &petsc_dest);
}

void
LBSVecOps::SetPrimarySTLvectorFromGSPETScVec(LBSProblem& lbs_problem,
                                             const LBSGroupset& groupset,
                                             Vec src,
                                             PhiSTLOption dest)
{
  auto& dest_phi =
    (dest == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();
  const double* petsc_src;
  VecGetArrayRead(src, &petsc_src);
  int64_t index = GroupsetScopedCopy(lbs_problem,
                                     groupset.groups.front().id,
                                     groupset.groups.size(),
                                     [&](int64_t idx, size_t mapped_idx)
                                     { dest_phi[mapped_idx] = petsc_src[idx]; });
  if (groupset.angle_agg)
  {
    if (dest == PhiSTLOption::PHI_NEW)
      groupset.angle_agg->SetNewDelayedAngularDOFsFromArray(index, petsc_src);
    else if (dest == PhiSTLOption::PHI_OLD)
      groupset.angle_agg->SetOldDelayedAngularDOFsFromArray(index, petsc_src);
  }
  VecRestoreArrayRead(src, &petsc_src);
}

void
LBSVecOps::SetGroupScopedPETScVecFromPrimarySTLvector(LBSProblem& lbs_problem,
                                                      int first_group_id,
                                                      int last_group_id,
                                                      Vec dest,
                                                      const std::vector<double>& src)
{
  double* petsc_dest;
  VecGetArray(dest, &petsc_dest);
  GroupsetScopedCopy(lbs_problem,
                     first_group_id,
                     last_group_id - first_group_id + 1,
                     [&](int64_t idx, size_t mapped_idx) { petsc_dest[idx] = src[mapped_idx]; });
  VecRestoreArray(dest, &petsc_dest);
}

void
LBSVecOps::SetPrimarySTLvectorFromGroupScopedPETScVec(
  LBSProblem& lbs_problem, int first_group_id, int last_group_id, Vec src, PhiSTLOption dest)
{
  auto& dest_phi =
    (dest == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();
  const double* petsc_src;
  VecGetArrayRead(src, &petsc_src);
  GroupsetScopedCopy(lbs_problem,
                     first_group_id,
                     last_group_id - first_group_id + 1,
                     [&](int64_t idx, size_t mapped_idx)
                     { dest_phi[mapped_idx] = petsc_src[idx]; });
  VecRestoreArrayRead(src, &petsc_src);
}

void
LBSVecOps::GSScopedCopyPrimarySTLvectors(LBSProblem& lbs_problem,
                                         const LBSGroupset& groupset,
                                         const std::vector<double>& src,
                                         std::vector<double>& dest)
{
  GroupsetScopedCopy(lbs_problem,
                     groupset.groups.front().id,
                     groupset.groups.size(),
                     [&](int64_t idx, size_t mapped_idx) { dest[mapped_idx] = src[mapped_idx]; });
}

void
LBSVecOps::GSScopedCopyPrimarySTLvectors(LBSProblem& lbs_problem,
                                         const LBSGroupset& groupset,
                                         PhiSTLOption src,
                                         PhiSTLOption dest)
{
  const auto& src_phi =
    (src == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();
  auto& dest_phi =
    (dest == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();
  GroupsetScopedCopy(lbs_problem,
                     groupset.groups.front().id,
                     groupset.groups.size(),
                     [&](int64_t idx, size_t mapped_idx)
                     { dest_phi[mapped_idx] = src_phi[mapped_idx]; });
  if (groupset.angle_agg)
  {
    if (src == PhiSTLOption::PHI_NEW and dest == PhiSTLOption::PHI_OLD)
      groupset.angle_agg->SetDelayedPsiOld2New();
    else if (src == PhiSTLOption::PHI_OLD and dest == PhiSTLOption::PHI_NEW)
      groupset.angle_agg->SetDelayedPsiNew2Old();
  }
}

void
LBSVecOps::SetMultiGSPETScVecFromPrimarySTLvector(LBSProblem& lbs_problem,
                                                  const std::vector<int>& groupset_ids,
                                                  Vec x,
                                                  PhiSTLOption which_phi)
{
  auto& y = (which_phi == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal()
                                                 : lbs_problem.GetPhiOldLocal();
  auto& groupsets = lbs_problem.GetGroupsets();
  double* x_ref;
  VecGetArray(x, &x_ref);

  for (int gs_id : groupset_ids)
  {
    const auto& groupset = groupsets.at(gs_id);
    int gsi = groupset.groups.front().id;
    int gsf = groupset.groups.back().id;
    int gss = gsf - gsi + 1;

    int64_t index = GroupsetScopedCopy(
      lbs_problem, gsi, gss, [&](int64_t idx, size_t mapped_idx) { x_ref[idx] = y[mapped_idx]; });
    if (groupset.angle_agg)
    {
      if (which_phi == PhiSTLOption::PHI_NEW)
        groupset.angle_agg->AppendNewDelayedAngularDOFsToArray(index, x_ref);
      else if (which_phi == PhiSTLOption::PHI_OLD)
        groupset.angle_agg->AppendOldDelayedAngularDOFsToArray(index, x_ref);
    }
  }

  VecRestoreArray(x, &x_ref);
}

void
LBSVecOps::SetPrimarySTLvectorFromMultiGSPETScVec(LBSProblem& lbs_problem,
                                                  const std::vector<int>& groupset_ids,
                                                  Vec x,
                                                  PhiSTLOption which_phi)
{
  auto& y = (which_phi == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal()
                                                 : lbs_problem.GetPhiOldLocal();
  auto& groupsets = lbs_problem.GetGroupsets();
  const double* x_ref;
  VecGetArrayRead(x, &x_ref);

  for (int gs_id : groupset_ids)
  {
    const auto& groupset = groupsets.at(gs_id);
    int gsi = groupset.groups.front().id;
    int gsf = groupset.groups.back().id;
    int gss = gsf - gsi + 1;

    int64_t index = GroupsetScopedCopy(
      lbs_problem, gsi, gss, [&](int64_t idx, size_t mapped_idx) { y[mapped_idx] = x_ref[idx]; });
    if (groupset.angle_agg)
    {
      if (which_phi == PhiSTLOption::PHI_NEW)
        groupset.angle_agg->SetNewDelayedAngularDOFsFromArray(index, x_ref);
      else if (which_phi == PhiSTLOption::PHI_OLD)
        groupset.angle_agg->SetOldDelayedAngularDOFsFromArray(index, x_ref);
    }
  }

  VecRestoreArrayRead(x, &x_ref);
}

} // namespace opensn

1	// SPDX-FileCopyrightText: 2024 The OpenSn Authors <https://open-sn.github.io/opensn/>
2	// SPDX-License-Identifier: MIT
3
4	#include "modules/linear_boltzmann_solvers/lbs_problem/lbs_vecops.h"
5	#include "modules/linear_boltzmann_solvers/lbs_problem/lbs_problem.h"
6
7	namespace opensn
8	{
9
10	template <typename Functor>
11	int
12	LBSVecOps::GroupsetScopedCopy(LBSProblem& lbs_problem, int gsi, int gss, Functor func)	66,499✔
13	{
14	CALI_CXX_MARK_SCOPE("LBSVecOps::GroupsetScopedCopy");	66,499✔
15
16	const auto& grid = lbs_problem.GetGrid();	66,499✔
17	const auto& cell_transport_views = lbs_problem.GetCellTransportViews();	66,499✔
18	auto num_moments = lbs_problem.GetNumMoments();	66,499✔
19
20	int64_t idx = -1;	66,499✔
21	for (const auto& cell : grid->local_cells)	20,776,153✔
22	{
23	const auto& transport_view = cell_transport_views[cell.local_id];	20,709,654✔
24	for (int i = 0; i < cell.vertex_ids.size(); ++i)	134,988,720✔
25	{
26	for (size_t m = 0; m < num_moments; ++m)	431,690,552✔
27	{
28	auto mapped_idx = transport_view.MapDOF(i, m, gsi);	317,411,486✔
29	for (int g = 0; g < gss; ++g)	2,147,483,647✔
30	{
31	++idx;	2,147,483,647✔
32	func(idx, mapped_idx + g);	2,147,483,647✔
33	}
34	}
35	}
36	}
37	return idx;	132,998✔
38	}	66,499✔
39
40	void
41	LBSVecOps::SetPhiVectorScalarValues(LBSProblem& lbs_problem, PhiSTLOption phi_opt, double value)	59✔
42	{
43	CALI_CXX_MARK_SCOPE("LBSVecOps::SetPhiVectorScalarValues");	59✔
44
45	auto& phi = (phi_opt == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal()	59✔
46	: lbs_problem.GetPhiOldLocal();	59✔
47	const auto& grid = lbs_problem.GetGrid();	59✔
48	const auto& groups = lbs_problem.GetGroups();	59✔
49	const auto& sdm = lbs_problem.GetSpatialDiscretization();	59✔
50	const auto& unknown_manager = lbs_problem.GetUnknownManager();	59✔
51
52	const size_t first_grp = groups.front().id;	59✔
53	const size_t final_grp = groups.back().id;	59✔
54
55	for (const auto& cell : grid->local_cells)	10,439✔
56	{
57	const auto& cell_mapping = sdm.GetCellMapping(cell);	10,380✔
58	const size_t num_nodes = cell_mapping.GetNumNodes();	10,380✔
59
60	for (size_t i = 0; i < num_nodes; ++i)	53,772✔
61	{
62	const auto dof_map = sdm.MapDOFLocal(cell, i, unknown_manager, 0, 0);	43,392✔
63	std::fill(phi.begin() + dof_map + first_grp, phi.begin() + dof_map + final_grp + 1, value);	86,784✔
64	}
65	}
66	}	59✔
67
68	void
69	LBSVecOps::ScalePhiVector(LBSProblem& lbs_problem, PhiSTLOption phi_opt, double value)	8✔
70	{
71	CALI_CXX_MARK_SCOPE("LBSVecOps::ScalePhiVector");	8✔
72
73	auto& phi = (phi_opt == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal()	8✔
74	: lbs_problem.GetPhiOldLocal();	×
75	auto& groupsets = lbs_problem.GetGroupsets();	8✔
76
77	Scale(phi, value);	8✔
78	for (auto& groupset : groupsets)	16✔
79	{
80	if (groupset.angle_agg)	8✔
81	{
82	if (phi_opt == PhiSTLOption::PHI_NEW)	8✔
83	{
84	auto psi = groupset.angle_agg->GetNewDelayedAngularDOFsAsSTLVector();	8✔
85	Scale(psi, value);	8✔
86	groupset.angle_agg->SetNewDelayedAngularDOFsFromSTLVector(psi);	8✔
87	}	8✔
88	else if (phi_opt == PhiSTLOption::PHI_OLD)	×
89	{
90	auto psi = groupset.angle_agg->GetOldDelayedAngularDOFsAsSTLVector();	×
91	Scale(psi, value);	×
92	groupset.angle_agg->SetOldDelayedAngularDOFsFromSTLVector(psi);	×
93	}	×
94	}
95	}
96	}	8✔
97
98	void
99	LBSVecOps::SetGSPETScVecFromPrimarySTLvector(LBSProblem& lbs_problem,	26,267✔
100	const LBSGroupset& groupset,
101	Vec dest,
102	PhiSTLOption src)
103	{
104	const auto& src_phi =	26,267✔
105	(src == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();	26,267✔
106	double* petsc_dest;	26,267✔
107	VecGetArray(dest, &petsc_dest);	26,267✔
108	int64_t index = GroupsetScopedCopy(lbs_problem,	52,534✔
109	groupset.groups.front().id,	26,267✔
110	groupset.groups.size(),	26,267✔
111	[&](int64_t idx, size_t mapped_idx)	2,147,483,647✔
112	{ petsc_dest[idx] = src_phi[mapped_idx]; });	2,147,483,647✔
113	if (groupset.angle_agg)	26,267✔
114	{
115	if (src == PhiSTLOption::PHI_NEW)	26,267✔
116	groupset.angle_agg->AppendNewDelayedAngularDOFsToArray(index, petsc_dest);	25,119✔
117	else if (src == PhiSTLOption::PHI_OLD)	1,148✔
118	groupset.angle_agg->AppendOldDelayedAngularDOFsToArray(index, petsc_dest);	1,148✔
119	}
120	VecRestoreArray(dest, &petsc_dest);	26,267✔
121	}	26,267✔
122
123	void
124	LBSVecOps::SetPrimarySTLvectorFromGSPETScVec(LBSProblem& lbs_problem,	26,267✔
125	const LBSGroupset& groupset,
126	Vec src,
127	PhiSTLOption dest)
128	{
129	auto& dest_phi =	26,267✔
130	(dest == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();	26,267✔
131	const double* petsc_src;	26,267✔
132	VecGetArrayRead(src, &petsc_src);	26,267✔
133	int64_t index = GroupsetScopedCopy(lbs_problem,	52,534✔
134	groupset.groups.front().id,	26,267✔
135	groupset.groups.size(),	26,267✔
136	[&](int64_t idx, size_t mapped_idx)	2,147,483,647✔
137	{ dest_phi[mapped_idx] = petsc_src[idx]; });	2,147,483,647✔
138	if (groupset.angle_agg)	26,267✔
139	{
140	if (dest == PhiSTLOption::PHI_NEW)	26,267✔
141	groupset.angle_agg->SetNewDelayedAngularDOFsFromArray(index, petsc_src);	8,828✔
142	else if (dest == PhiSTLOption::PHI_OLD)	17,439✔
143	groupset.angle_agg->SetOldDelayedAngularDOFsFromArray(index, petsc_src);	17,439✔
144	}
145	VecRestoreArrayRead(src, &petsc_src);	26,267✔
146	}	26,267✔
147
148	void
149	LBSVecOps::SetGroupScopedPETScVecFromPrimarySTLvector(LBSProblem& lbs_problem,	×
150	int first_group_id,
151	int last_group_id,
152	Vec dest,
153	const std::vector<double>& src)
154	{
155	double* petsc_dest;	×
156	VecGetArray(dest, &petsc_dest);	×
157	GroupsetScopedCopy(lbs_problem,	×
158	first_group_id,
159	last_group_id - first_group_id + 1,	×
160	[&](int64_t idx, size_t mapped_idx) { petsc_dest[idx] = src[mapped_idx]; });	×
161	VecRestoreArray(dest, &petsc_dest);	×
162	}	×
163
164	void
165	LBSVecOps::SetPrimarySTLvectorFromGroupScopedPETScVec(	28✔
166	LBSProblem& lbs_problem, int first_group_id, int last_group_id, Vec src, PhiSTLOption dest)
167	{
168	auto& dest_phi =	28✔
169	(dest == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();	28✔
170	const double* petsc_src;	28✔
171	VecGetArrayRead(src, &petsc_src);	28✔
172	GroupsetScopedCopy(lbs_problem,	28✔
173	first_group_id,
174	last_group_id - first_group_id + 1,	28✔
175	[&](int64_t idx, size_t mapped_idx)	2,217,800✔
176	{ dest_phi[mapped_idx] = petsc_src[idx]; });	2,217,800✔
177	VecRestoreArrayRead(src, &petsc_src);	28✔
178	}	28✔
179
180	void
181	LBSVecOps::GSScopedCopyPrimarySTLvectors(LBSProblem& lbs_problem,	47✔
182	const LBSGroupset& groupset,
183	const std::vector<double>& src,
184	std::vector<double>& dest)
185	{
186	GroupsetScopedCopy(lbs_problem,	47✔
187	groupset.groups.front().id,	47✔
188	groupset.groups.size(),	47✔
189	[&](int64_t idx, size_t mapped_idx) { dest[mapped_idx] = src[mapped_idx]; });	759,888✔
190	}	47✔
191
192	void
193	LBSVecOps::GSScopedCopyPrimarySTLvectors(LBSProblem& lbs_problem,	6,190✔
194	const LBSGroupset& groupset,
195	PhiSTLOption src,
196	PhiSTLOption dest)
197	{
198	const auto& src_phi =	6,190✔
199	(src == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();	6,190✔
200	auto& dest_phi =	6,190✔
201	(dest == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal() : lbs_problem.GetPhiOldLocal();	6,190✔
202	GroupsetScopedCopy(lbs_problem,	6,190✔
203	groupset.groups.front().id,	6,190✔
204	groupset.groups.size(),	6,190✔
205	[&](int64_t idx, size_t mapped_idx)	1,396,369,502✔
206	{ dest_phi[mapped_idx] = src_phi[mapped_idx]; });	1,396,369,502✔
207	if (groupset.angle_agg)	6,190✔
208	{
209	if (src == PhiSTLOption::PHI_NEW and dest == PhiSTLOption::PHI_OLD)	6,190✔
210	groupset.angle_agg->SetDelayedPsiOld2New();	6,190✔
211	else if (src == PhiSTLOption::PHI_OLD and dest == PhiSTLOption::PHI_NEW)	×
212	groupset.angle_agg->SetDelayedPsiNew2Old();	×
213	}
214	}	6,190✔
215
216	void
217	LBSVecOps::SetMultiGSPETScVecFromPrimarySTLvector(LBSProblem& lbs_problem,	3,864✔
218	const std::vector<int>& groupset_ids,
219	Vec x,
220	PhiSTLOption which_phi)
221	{
222	auto& y = (which_phi == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal()	3,864✔
223	: lbs_problem.GetPhiOldLocal();	28✔
224	auto& groupsets = lbs_problem.GetGroupsets();	3,864✔
225	double* x_ref;	3,864✔
226	VecGetArray(x, &x_ref);	3,864✔
227
228	for (int gs_id : groupset_ids)	7,728✔
229	{
230	const auto& groupset = groupsets.at(gs_id);	3,864✔
231	int gsi = groupset.groups.front().id;	3,864✔
232	int gsf = groupset.groups.back().id;	3,864✔
233	int gss = gsf - gsi + 1;	3,864✔
234
235	int64_t index = GroupsetScopedCopy(	3,864✔
236	lbs_problem, gsi, gss, [&](int64_t idx, size_t mapped_idx) { x_ref[idx] = y[mapped_idx]; });	372,380,000✔
237	if (groupset.angle_agg)	3,864✔
238	{
239	if (which_phi == PhiSTLOption::PHI_NEW)	3,864✔
240	groupset.angle_agg->AppendNewDelayedAngularDOFsToArray(index, x_ref);	3,836✔
241	else if (which_phi == PhiSTLOption::PHI_OLD)	28✔
242	groupset.angle_agg->AppendOldDelayedAngularDOFsToArray(index, x_ref);	28✔
243	}
244	}
245
246	VecRestoreArray(x, &x_ref);	3,864✔
247	}	3,864✔
248
249	void
250	LBSVecOps::SetPrimarySTLvectorFromMultiGSPETScVec(LBSProblem& lbs_problem,	3,836✔
251	const std::vector<int>& groupset_ids,
252	Vec x,
253	PhiSTLOption which_phi)
254	{
255	auto& y = (which_phi == PhiSTLOption::PHI_NEW) ? lbs_problem.GetPhiNewLocal()	3,836✔
256	: lbs_problem.GetPhiOldLocal();	3,836✔
257	auto& groupsets = lbs_problem.GetGroupsets();	3,836✔
258	const double* x_ref;	3,836✔
259	VecGetArrayRead(x, &x_ref);	3,836✔
260
261	for (int gs_id : groupset_ids)	7,672✔
262	{
263	const auto& groupset = groupsets.at(gs_id);	3,836✔
264	int gsi = groupset.groups.front().id;	3,836✔
265	int gsf = groupset.groups.back().id;	3,836✔
266	int gss = gsf - gsi + 1;	3,836✔
267
268	int64_t index = GroupsetScopedCopy(	3,836✔
269	lbs_problem, gsi, gss, [&](int64_t idx, size_t mapped_idx) { y[mapped_idx] = x_ref[idx]; });	370,162,200✔
270	if (groupset.angle_agg)	3,836✔
271	{
272	if (which_phi == PhiSTLOption::PHI_NEW)	3,836✔
273	groupset.angle_agg->SetNewDelayedAngularDOFsFromArray(index, x_ref);	×
274	else if (which_phi == PhiSTLOption::PHI_OLD)	3,836✔
275	groupset.angle_agg->SetOldDelayedAngularDOFsFromArray(index, x_ref);	3,836✔
276	}
277	}
278
279	VecRestoreArrayRead(x, &x_ref);	3,836✔
280	}	3,836✔
281
282	} // namespace opensn

Open-Sn / opensn / 18380674863

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous