Vous avez reçu un message "Your GitLab account has been locked ..." ? Pas d'inquiétude : lisez cet article https://docs.gricad-pages.univ-grenoble-alpes.fr/help/unlock/

Commit 2c36367f authored by Olga Stamati's avatar Olga Stamati
Browse files

[skip-ci] closes #189

parent 840068fc
Pipeline #52295 passed with stages
in 24 minutes and 14 seconds
......@@ -122,7 +122,10 @@ if not args.Q8:
fileNameBase = args.OUT_DIR+"/"+args.PREFIX+"-strain-Geers"
mode = "Geers"
else:
fileNameBase = args.OUT_DIR+"/"+args.PREFIX+"-strain-Q8"
if twoD:
fileNameBase = args.OUT_DIR+"/"+args.PREFIX+"-strain-Q4"
else:
fileNameBase = args.OUT_DIR+"/"+args.PREFIX+"-strain-Q8"
mode = "Q8"
# Save strain fields
......
"""
"""
Library of SPAM functions for dealing with fields of Phi or fields of F
Copyright (C) 2020 SPAM Contributors
......@@ -418,7 +418,7 @@ def decomposeFfield(Ffield, components, twoD=False):
output[component] = numpy.zeros(fieldRavelLength)
if component == 'r' or component == 'z':
output[component] = numpy.zeros((fieldRavelLength, 3))
if component == 'U' or component == 'E' or component == 'e':
if component == 'U' or component == 'e':
output[component] = numpy.zeros((fieldRavelLength, 3, 3))
# Iterate through flat field of Fs
......@@ -436,7 +436,7 @@ def decomposeFfield(Ffield, components, twoD=False):
if component == 'r' or component == 'z':
output[component] = numpy.array(output[component]).reshape(Ffield.shape[0:-1])
if component == 'U' or component == 'E' or component == 'e':
if component == 'U' or component == 'e':
output[component] = numpy.array(output[component]).reshape(Ffield.shape)
return output
......@@ -483,7 +483,7 @@ def decomposePhiField(PhiField, components, twoD=False):
output[component] = numpy.zeros(fieldRavelLength)
if component == 't' or component == 'r' or component == 'z':
output[component] = numpy.zeros((fieldRavelLength, 3))
if component == 'U' or component == 'E' or component == 'e':
if component == 'U' or component == 'e':
output[component] = numpy.zeros((fieldRavelLength, 3, 3))
# Iterate through flat field of Phis
......@@ -501,7 +501,7 @@ def decomposePhiField(PhiField, components, twoD=False):
if component == 't' or component == 'r' or component == 'z':
output[component] = numpy.array(output[component]).reshape(*PhiField.shape[0:-2], 3)
if component == 'U' or component == 'E' or component == 'e':
if component == 'U' or component == 'e':
output[component] = numpy.array(output[component]).reshape(*PhiField.shape[0:-2], 3, 3)
return output
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment