Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
MoReFEM
CoreLibrary
MoReFEM
Commits
f7f0cd35
Commit
f7f0cd35
authored
Nov 04, 2015
by
GILLES Sebastien
Browse files
#723
Remove dev prints.
parent
5e47c222
Changes
2
Hide whitespace changes
Inline
Side-by-side
Sources/ModelInstances/FSI_EI/Newton/ModelForward.hxx
View file @
f7f0cd35
...
...
@@ -83,17 +83,6 @@ namespace HappyHeart
__FILE__
,
__LINE__
);
dh_displacement_
->
ZeroEntries
(
__FILE__
,
__LINE__
);
{
const
auto
&
foo
=
GetSolidDisplacementEvolutionInSoF
();
const
PetscReal
evaluation_state_min
=
foo
.
Min
(
__FILE__
,
__LINE__
).
second
;
const
PetscReal
evaluation_state_max
=
foo
.
Max
(
__FILE__
,
__LINE__
).
second
;
std
::
cout
<<
"EVO EXTREMA ARE "
<<
evaluation_state_min
<<
" and "
<<
evaluation_state_max
<<
" on "
<<
foo
.
GetProgramWiseSize
(
__FILE__
,
__LINE__
)
<<
" elements."
<<
std
::
endl
;
// OK
}
dH_solver_
->
SolveLinear
(
test
,
test
,
GetSolidDisplacementEvolutionInSoF
(),
*
dh_displacement_
,
// \todo #723 Inverted for tests
...
...
@@ -108,8 +97,6 @@ namespace HappyHeart
std
::
cout
<<
"AFTER GMRES EXTREMA ARE "
<<
evaluation_state_min
<<
" and "
<<
evaluation_state_max
<<
" on "
<<
foo
.
GetProgramWiseSize
(
__FILE__
,
__LINE__
)
<<
" elements."
<<
std
::
endl
;
// OK
}
// dH(*dh_displacement_);
}
...
...
@@ -122,18 +109,6 @@ namespace HappyHeart
input_displ_vector
.
SetFromPetscVec
(
petsc_input_vector
,
__FILE__
,
__LINE__
);
{
const
PetscReal
evaluation_state_min
=
input_displ_vector
.
Min
(
__FILE__
,
__LINE__
).
second
;
const
PetscReal
evaluation_state_max
=
input_displ_vector
.
Max
(
__FILE__
,
__LINE__
).
second
;
std
::
cout
<<
"Input DISPL EXTREMA ARE "
<<
evaluation_state_min
<<
" and "
<<
evaluation_state_max
<<
" on "
<<
input_displ_vector
.
GetProgramWiseSize
(
__FILE__
,
__LINE__
)
<<
" elements."
<<
std
::
endl
<<
std
::
endl
;
}
// Update Dirichlet boundary condition.
{
auto
&
displ_vector_on_interface
=
*
work_vector_on_interface_1_
;
...
...
@@ -151,28 +126,8 @@ namespace HappyHeart
auto
&
fluid_varf
=
GetNonCstImplicitStepFluidVariationalFormulation
();
{
const
auto
&
fluid_sol
=
fluid_varf
.
GetSystemSolution
(
fluid_varf
.
GetNumberingSubset
());
const
PetscReal
evaluation_state_min
=
fluid_sol
.
Min
(
__FILE__
,
__LINE__
).
second
;
const
PetscReal
evaluation_state_max
=
fluid_sol
.
Max
(
__FILE__
,
__LINE__
).
second
;
std
::
cout
<<
"
\n
delta[Vel,p] EXTREMA ARE "
<<
evaluation_state_min
<<
" and "
<<
evaluation_state_max
<<
" on "
<<
fluid_sol
.
GetProgramWiseSize
(
__FILE__
,
__LINE__
)
<<
" elements."
<<
std
::
endl
<<
std
::
endl
;
// OK Within solve!
}
const
auto
&
differential_fluid_residual
=
fluid_varf
.
ComputeResidual
(
-
1.
);
// \todo #723 Hack!
{
const
PetscReal
evaluation_state_min
=
differential_fluid_residual
.
Min
(
__FILE__
,
__LINE__
).
second
;
const
PetscReal
evaluation_state_max
=
differential_fluid_residual
.
Max
(
__FILE__
,
__LINE__
).
second
;
std
::
cout
<<
"differential_fluid_residual EXTREMA ARE "
<<
evaluation_state_min
<<
" and "
<<
evaluation_state_max
<<
" on "
<<
differential_fluid_residual
.
GetProgramWiseSize
(
__FILE__
,
__LINE__
)
<<
" elements."
<<
std
::
endl
<<
std
::
endl
;
// opposite of the Freefem output
}
const
auto
&
solid_varf
=
SolidVariationalFormulationPolicyT
::
GetVariationalFormulation
();
const
auto
&
solid_numbering_subset
=
solid_varf
.
GetNumberingSubset
();
...
...
@@ -184,19 +139,8 @@ namespace HappyHeart
solid_residual
,
__FILE__
,
__LINE__
);
{
const
PetscReal
evaluation_state_min
=
solid_residual
.
Min
(
__FILE__
,
__LINE__
).
second
;
const
PetscReal
evaluation_state_max
=
solid_residual
.
Max
(
__FILE__
,
__LINE__
).
second
;
std
::
cout
<<
"delta coupling before Dirichlet EXTREMA ARE "
<<
evaluation_state_min
<<
" and "
<<
evaluation_state_max
<<
" on "
<<
solid_residual
.
GetProgramWiseSize
(
__FILE__
,
__LINE__
)
<<
" elements."
<<
std
::
endl
<<
std
::
endl
;
// opposite of the Freefem output
}
const
auto
&
solid_tangent
=
solid_varf
.
GetSystemMatrix
(
solid_numbering_subset
,
solid_numbering_subset
);
// solid_tangent.View(mpi,__FILE__, __LINE__);
const
auto
&
god_of_dof
=
parent
::
GetGodOfDof
(
EnumUnderlyingType
(
MeshIndex
::
mesh
));
auto
&
dirichlet_bc
=
...
...
@@ -209,28 +153,11 @@ namespace HappyHeart
solution
,
__FILE__
,
__LINE__
);
{
const
PetscReal
evaluation_state_min
=
solution
.
Min
(
__FILE__
,
__LINE__
).
second
;
const
PetscReal
evaluation_state_max
=
solution
.
Max
(
__FILE__
,
__LINE__
).
second
;
std
::
cout
<<
"diff solid solution EXTREMA ARE "
<<
evaluation_state_min
<<
" and "
<<
evaluation_state_max
<<
" on "
<<
solution
.
GetProgramWiseSize
(
__FILE__
,
__LINE__
)
<<
" elements."
<<
std
::
endl
<<
std
::
endl
;
// opposite of the Freefem output
}
Wrappers
::
Petsc
::
AXPY
(
-
1.
,
solution
,
input_displ_vector
,
__FILE__
,
__LINE__
);
{
const
PetscReal
evaluation_state_min
=
input_displ_vector
.
Min
(
__FILE__
,
__LINE__
).
second
;
const
PetscReal
evaluation_state_max
=
input_displ_vector
.
Max
(
__FILE__
,
__LINE__
).
second
;
std
::
cout
<<
"End of dH extrema ARE "
<<
evaluation_state_min
<<
" and "
<<
evaluation_state_max
<<
" on "
<<
input_displ_vector
.
GetProgramWiseSize
(
__FILE__
,
__LINE__
)
<<
" elements."
<<
std
::
endl
;
// opposite of the Freefem output
}
// We can't take directly the internal pointer of a wrapped vector which pattern would be similar to input_displ_vector
// : residual computation is wrong if we do.
auto
error_code
=
VecCopy
(
input_displ_vector
.
Internal
(),
petsc_output_vector
);
...
...
@@ -360,9 +287,6 @@ namespace HappyHeart
implicit_step_fluid_formulation
.
WriteSolution
(
time_manager
,
numbering_subset
);
#ifndef NDEBUG
std
::
cout
<<
"NON DIFF FLUID VARF -> "
<<
implicit_step_fluid_formulation
.
GetSystemRhs
(
numbering_subset
).
GetProgramWiseSize
(
__FILE__
,
__LINE__
);
Wrappers
::
Petsc
::
AssertNumericValues
(
implicit_step_fluid_formulation
.
GetSystemRhs
(
numbering_subset
),
__FILE__
,
__LINE__
);
Wrappers
::
Petsc
::
AssertNumericValues
(
implicit_step_fluid_formulation
.
GetSystemSolution
(
numbering_subset
),
__FILE__
,
__LINE__
);
#endif // NDEBUG
...
...
Sources/ThirdParty/Wrappers/Petsc/Matrix/ShellMatrix.hxx
View file @
f7f0cd35
...
...
@@ -60,7 +60,6 @@ namespace HappyHeart
if
(
error_code
)
throw
ExceptionNS
::
Exception
(
error_code
,
"MatShellSetOperation"
,
invoking_file
,
invoking_line
);
}
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new 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