Homec4science
Diffusion tamaas (features/netcdf-mpi) Browse

tamaas/tests/test_patch_plasticity.py
998dbf9a8e3afeatures/netcdf-mpi

Last Change
Raw File

test_patch_plasticity.py
View Options

# -*- coding: utf-8 -*-
# @file
# @section LICENSE
#
# Copyright (©) 2016-20 EPFL (École Polytechnique Fédérale de Lausanne),
# Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU Affero General Public License as published
# by the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Affero General Public License for more details.
#
# You should have received a copy of the GNU Affero General Public License
# along with this program. If not, see <https://www.gnu.org/licenses/>.
import pytest
from numpy.linalg import norm
import tamaas
from tamaas.nonlinear_solvers import \
DFSANESolver as PySolver, \
DFSANECXXSolver as CXXSolver, \
NewtonKrylovSolver
@pytest.fixture(params=[CXXSolver, NewtonKrylovSolver, PySolver])
def solvers(request):
return request.param
def test_patch_plasticity(patch_isotropic_plasticity, solvers):
"Test analyitical solution of 1d plasticity"
tamaas.set_log_level(tamaas.LogLevel.info)
Solver = solvers
model = patch_isotropic_plasticity.model
residual = patch_isotropic_plasticity.residual
applied_pressure = 0.1
solver = Solver(residual)
solver.tolerance = 1e-15
pressure = model['traction'][..., 2]
pressure[:] = applied_pressure
solver.solve()
solver.updateState()
solution, normal = patch_isotropic_plasticity.solution(applied_pressure)
for key in solution:
error = norm(model[key] - solution[key]) / normal[key]
assert error < 2e-15
c4science · Help