netcdf-mpi

test_voigt.py
View Options

			"""
			@author Lucas Frérot <lucas.frerot@epfl.ch>

			@section LICENSE

			Copyright (©) 2019 EPFL (Ecole Polytechnique Fédérale de
			Lausanne) Laboratory (LSMS - Laboratoire de Simulation en Mécanique des
			Solides)

			Tamaas is free software: you can redistribute it and/or modify it under the
			terms of the GNU Lesser General Public License as published by the Free
			Software Foundation, either version 3 of the License, or (at your option) any
			later version.

			Tamaas 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 Lesser General Public License for more
			details.

			You should have received a copy of the GNU Lesser General Public License
			along with Tamaas. If not, see <http://www.gnu.org/licenses/>.
			"""

			import tamaas as tm
			import numpy as np


			def test_voigt(tamaas_fixture):
			shape = [3, 3, 3]
			stress = np.ones(shape + [9])
			#voigt = np.zeros(shape + [6])
			voigt = tm.to_voigt(stress)
			sol = np.ones(shape + [6])
			sol[:, :, :, 3:] = np.sqrt(2)
			print(voigt)
			print(sol)
			assert np.all(np.abs(sol - voigt) < 1e-15)