Coverage for tests / tests_library / tests_model_components / test_velocity_law_hht.py: 100%

1# SPDX-FileCopyrightText: Copyright INRIA 

2# 

3# SPDX-License-Identifier: LGPL-3.0-only 

4# 

5# Copyright INRIA 

6# 

7# This file is part of PhysioBlocks, a library mostly developed by the 

8# [Ananke project-team](https://team.inria.fr/ananke) at INRIA. 

9# 

10# Authors: 

11# - Colin Drieu 

12# - Dominique Chapelle 

13# - François Kimmig 

14# - Philippe Moireau 

15# 

16# PhysioBlocks is free software: you can redistribute it and/or modify it under the 

17# terms of the GNU Lesser General Public License as published by the Free Software 

18# Foundation, version 3 of the License. 

19# 

20# PhysioBlocks is distributed in the hope that it will be useful, but WITHOUT ANY 

21# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A 

22# PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details. 

23# 

24# You should have received a copy of the GNU Lesser General Public License along with 

25# PhysioBlocks. If not, see <https://www.gnu.org/licenses/>. 

26 

27import numpy as np 

28import pytest 

29 

30from physioblocks.computing.quantities import Quantity 

31from physioblocks.library.model_components.velocity_law import ( 

32 VelocityLawHHTModelComponent, 

33) 

34from physioblocks.simulation.state import State 

35from physioblocks.simulation.time_manager import Time 

36from physioblocks.utils.gradient_test_utils import gradient_test_from_model 

37 

38 

39@pytest.fixture 

40def ref_block() -> VelocityLawHHTModelComponent: 

41 return VelocityLawHHTModelComponent( 

42 disp=Quantity(0.1), 

43 vel=Quantity(0.2), 

44 accel=Quantity(0.3), 

45 scheme_ts_hht=Quantity(0.4), 

46 time=Time(0.0), 

47 ) 

48 

49 

50class TestVelocityLawHHTModelComponent: 

51 def test_check_gradient(self, ref_block: VelocityLawHHTModelComponent): 

52 state = State() 

53 state["disp"] = ref_block.disp 

54 state["vel"] = ref_block.vel 

55 state["accel"] = ref_block.accel 

56 ref_block.time.update(0.001) 

57 

58 magnitudes = np.array([1.0, 1.0, 1.0]) 

59 

60 assert gradient_test_from_model(ref_block, state, magnitudes)