Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
GILLES Sebastien
MoReFEM
Commits
c2055d13
Commit
c2055d13
authored
Apr 13, 2019
by
GILLES Sebastien
Browse files
#859 Remove debug lines.
parent
f8cb21ac
Changes
1
Hide whitespace changes
Inline
Side-by-side
Sources/FiniteElement/FiniteElementSpace/FEltSpace.cpp
View file @
c2055d13
...
...
@@ -194,15 +194,9 @@ namespace MoReFEM
{
assert
(
!
(
!
geometric_type_ptr
));
if
(
domain
.
GetUniqueId
()
==
3
)
std
::
cout
<<
"FELT 0"
<<
std
::
endl
;
if
(
!
domain
.
DoRefGeomEltMatchCriteria
(
*
geometric_type_ptr
))
continue
;
if
(
domain
.
GetUniqueId
()
==
3
)
std
::
cout
<<
"FELT 10"
<<
std
::
endl
;
auto
&&
ref_felt_space_ptr
=
std
::
make_unique
<
const
Internal
::
RefFEltNS
::
RefLocalFEltSpace
>
(
geometric_type_ptr
,
unknown_storage
,
...
...
@@ -220,14 +214,9 @@ namespace MoReFEM
static_cast
<
std
::
size_t
>
(
geometric_elt_range
.
second
-
geometric_elt_range
.
first
);
local_felt_space_list
.
reserve
(
Ngeometric_elt_in_range
);
if
(
domain
.
GetUniqueId
()
==
3
)
std
::
cout
<<
"FELT 20"
<<
std
::
endl
;
for
(
auto
it
=
geometric_elt_range
.
first
;
it
!=
geometric_elt_range
.
second
;
++
it
)
{
if
(
domain
.
GetUniqueId
()
==
3
)
std
::
cout
<<
"FELT 30"
<<
std
::
endl
;
const
auto
&
geometric_elt_ptr
=
*
it
;
assert
(
!
(
!
geometric_elt_ptr
));
const
auto
&
geometric_elt
=
*
geometric_elt_ptr
;
...
...
@@ -235,22 +224,15 @@ namespace MoReFEM
if
(
domain
.
IsGeometricEltInside
(
geometric_elt
))
{
if
(
domain
.
GetUniqueId
()
==
3
)
std
::
cout
<<
"FELT 35"
<<
std
::
endl
;
auto
&&
ptr
=
new
LocalFEltSpace
(
*
ref_felt_space_ptr
,
geometric_elt_ptr
);
local_felt_space_list
.
insert
({
geometric_elt
.
GetIndex
(),
LocalFEltSpace
::
shared_ptr
(
ptr
)});
}
}
if
(
domain
.
GetUniqueId
()
==
3
)
std
::
cout
<<
"FELT 40"
<<
std
::
endl
;
assert
(
Ngeometric_elt_in_range
>=
local_felt_space_list
.
size
());
// equal if no domain restriction applied.
if
(
local_felt_space_list
.
empty
())
// might happen if for instance RefGeomElt not considered in the domain.
continue
;
if
(
domain
.
GetUniqueId
()
==
3
)
std
::
cout
<<
"FELT 50"
<<
std
::
endl
;
felt_list_per_type
.
push_back
(
std
::
make_pair
(
std
::
move
(
ref_felt_space_ptr
),
local_felt_space_list
));
}
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment