Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F85654436
test_patch_westergaard.py
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Mon, Sep 30, 19:15
Size
1 KB
Mime Type
text/x-python
Expires
Wed, Oct 2, 19:15 (2 d)
Engine
blob
Format
Raw Data
Handle
21172825
Attached To
rTAMAAS tamaas
test_patch_westergaard.py
View Options
#!/usr/bin/env python
# coding: utf-8
# -----------------------------------------------------------------------------
# @author Lucas Frérot <lucas.frerot@epfl.ch>
#
# @section LICENSE
#
# Copyright (©) 2016 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/>.
# -----------------------------------------------------------------------------
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
)
Event Timeline
Log In to Comment