test_patch_westergaard.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	# -- coding: utf-8 --
	# @file
	# @section LICENSE
	#
	# Copyright (©) 2016-19 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/>.

	from __future__ import division, print_function

	from numpy.linalg import norm


	def test_patch_westergaard(patch_westergaard):
	model = patch_westergaard.model
	solution = patch_westergaard.solution
	pressure = patch_westergaard.pressure
	model.getTraction()[...] = pressure[...]
	model.solveNeumann()
	output_displ = model.getDisplacement()

	error = norm(solution - output_displ) / norm(solution)

	assert error < 1e-15 and norm(solution) > 1e-15, \
	"Neumann error = {}".format(error)

	output_displ[:, :] = solution[:, :]
	model.solveDirichlet()
	error = norm(pressure - model.getTraction()) / norm(pressure)

	assert error < 1e-15, "Dirichlet error = {}".format(error)