Skip to content
Snippets Groups Projects
Commit a7a6c502 authored by Tom Fischer's avatar Tom Fischer
Browse files

[PL] std::fabs -> std::abs.

parent 9ccd4f8c
No related branches found
No related tags found
No related merge requests found
...@@ -64,7 +64,7 @@ bool Output::shallDoOutput(int timestep, double const t) ...@@ -64,7 +64,7 @@ bool Output::shallDoOutput(int timestep, double const t)
auto const fixed_output_time = std::lower_bound( auto const fixed_output_time = std::lower_bound(
cbegin(_fixed_output_times), cend(_fixed_output_times), t); cbegin(_fixed_output_times), cend(_fixed_output_times), t);
if ((fixed_output_time != cend(_fixed_output_times)) && if ((fixed_output_time != cend(_fixed_output_times)) &&
(std::fabs(*fixed_output_time - t) < (std::abs(*fixed_output_time - t) <
std::numeric_limits<double>::epsilon())) std::numeric_limits<double>::epsilon()))
{ {
return true; return true;
......
...@@ -319,7 +319,7 @@ void PhaseFieldProcess<DisplacementDim>::postNonLinearSolverConcreteProcess( ...@@ -319,7 +319,7 @@ void PhaseFieldProcess<DisplacementDim>::postNonLinearSolverConcreteProcess(
_process_data.pressure = _process_data.pressure =
_process_data.injected_volume / _process_data.crack_volume; _process_data.injected_volume / _process_data.crack_volume;
_process_data.pressure_error = _process_data.pressure_error =
std::fabs(_process_data.pressure_old - _process_data.pressure) / std::abs(_process_data.pressure_old - _process_data.pressure) /
_process_data.pressure; _process_data.pressure;
INFO("Internal pressure: {:g} and Pressure error: {:.4e}", INFO("Internal pressure: {:g} and Pressure error: {:.4e}",
_process_data.pressure, _process_data.pressure_error); _process_data.pressure, _process_data.pressure_error);
......
...@@ -458,7 +458,7 @@ double TimeLoop::computeTimeStepping(const double prev_dt, double& t, ...@@ -458,7 +458,7 @@ double TimeLoop::computeTimeStepping(const double prev_dt, double& t,
dt = NumLib::possiblyClampDtToNextFixedTime(t, dt, dt = NumLib::possiblyClampDtToNextFixedTime(t, dt,
_output->getFixedOutputTimes()); _output->getFixedOutputTimes());
// Check whether the time stepping is stabilized // Check whether the time stepping is stabilized
if (std::fabs(dt - prev_dt) < std::numeric_limits<double>::epsilon()) if (std::abs(dt - prev_dt) < std::numeric_limits<double>::epsilon())
{ {
if (_last_step_rejected) if (_last_step_rejected)
{ {
...@@ -656,7 +656,7 @@ bool TimeLoop::loop() ...@@ -656,7 +656,7 @@ bool TimeLoop::loop()
} }
dt = computeTimeStepping(prev_dt, t, accepted_steps, rejected_steps); dt = computeTimeStepping(prev_dt, t, accepted_steps, rejected_steps);
if (std::fabs(t - _end_time) < std::numeric_limits<double>::epsilon() || if (std::abs(t - _end_time) < std::numeric_limits<double>::epsilon() ||
t + dt > _end_time) t + dt > _end_time)
{ {
break; break;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment