def __init__(self, name=None):
Script.__init__(self, name)
DynamicComponentHarness.__init__(self)
self.client = ''
return
def __init__(self, name="slip_invert"):
Application.__init__(self, name)
self.dataCoords = None
self.dataVals = None
self.dataCov = None
self.dataNames = []
self.numDataPoints = 0
self.numDesignRows = 0
self.numFaultVerts = 0
self.numFaultCells = 0
self.faultVertCoords = None
self.faultCells = None
self.numImpulses = 0
self.impulseInds = None
self.impulseCoords = None
self.summaryHead = 'Penalty-weight\tData-residual\t' + \
'Weighted-data-residual\tPenalty-residual\t' + \
'Weighted-penalty-residual\tTotal-residual\t' + \
'Total-weighted-residual'
self.numSummaryCols = 7
self.design = None
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
# import pdb
# pdb.set_trace()
# File info.
self.vtkInputFile = self.inventory.vtkInputFile
self.vtkOutputFile = self.inventory.vtkOutputFile
# Index information
self.vtkTensorIndex = self.inventory.vtkTensorIndex
self.vtkTensorComponentsOrder = self.inventory.vtkTensorComponentsOrder
# Regional field
s1 = float(self.inventory.regionalSigma1[0])
s2 = float(self.inventory.regionalSigma1[0])
s3 = float(self.inventory.regionalSigma1[0])
sVec = numpy.array([s1, s2, s3], dtype=numpy.float64)
self.regionalSigma = numpy.diag(sVec)
self.regionalAxes = numpy.zeros((3, 3), dtype=numpy.float64)
for i in range(3):
self.regionalAxes[0,
i] = float(self.inventory.regionalSigma1[i + 1])
self.regionalAxes[1,
i] = float(self.inventory.regionalSigma2[i + 1])
self.regionalAxes[2,
i] = float(self.inventory.regionalSigma3[i + 1])
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
self.snapshots = list(map(float, self.snapshots))
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
self.dataDim = self.inventory.dataDim
self.pointsFile = self.inventory.pointsFile
self.segsFile = self.inventory.segsFile
self.slipScale = self.inventory.slipScale
self.spaceDim = 3
self.pointsSpatialDB = self.inventory.pointsSpatialDB
self.upDir = self.inventory.upDir
self.normalDir = self.inventory.normalDir
self.dipSlip = self.inventory.dipSlip
self.dipCutoff = self.inventory.dipCutoff
self.xMinVal = self.inventory.xMin.value
self.xMaxVal = self.inventory.xMax.value
self.yMinVal = self.inventory.yMin.value
self.yMaxVal = self.inventory.yMax.value
self.zMinVal = self.inventory.zMin.value
self.zMaxVal = self.inventory.zMax.value
self.defaultValues = self.inventory.defaultValues
self.upVec = numpy.array([float(self.upDir[0]), float(self.upDir[1]),
float(self.upDir[2])], dtype=numpy.float64)
self.normalVec = numpy.array([float(self.normalDir[0]),
float(self.normalDir[1]),
float(self.normalDir[2])], dtype=numpy.float64)
self.dipCutoffProj = abs(math.sin(self.dipCutoff.value))
return
def __init__(self, name="transform"):
Application.__init__(self, name)
self.numPoints = 0
self.numSegs = 0
self.segInfo = None
self.pointsUTM = []
self.normals = []
return
def __init__(self, name, asCGI=None):
Script.__init__(self, name)
self.stream = None
if asCGI is None:
asCGI = True
self.asCGIScript(asCGI)
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
self.polesInputFile = self.inventory.polesInputFile
self.polesOutputFile = self.inventory.polesOutputFile
self.geogToCart = self.inventory.geogToCart
self.spaceDim = 3
return
def __init__(self, name="euler"):
Application.__init__(self, name)
self.numPoints = 0
self.pointsUTM = []
self.normals = []
self.eulerPole = numpy.array([0.0, 0.0, 0.0], dtype=numpy.float64)
# Note that we use a mean radius since we are doing rotations on a
# spherical Earth.
self.earthRad = 6372795.477598
return
def __init__(self, name="generate_slowslip"):
Application.__init__(self, name)
self.lon = None
self.lat = None
self.z = None
self.grid = None
self.faultSlip = None
return
def _configure(self):
Application._configure(self)
self.modeler = self.inventory.modeler
self.contours = self.inventory.contours
self.extender = self.inventory.extender
self.slabTopFilename = self.inventory.slabTopFilename
self.slabBotFilename = self.inventory.slabBotFilename
self.splayFilename = self.inventory.splayFilename
self.slabThickness = self.inventory.slabThickness
self.slabNormalDir = self.inventory.slabNormalDir
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
self.ucdFaceFile = self.inventory.ucdFaceFile
self.faultIDNum = self.inventory.faultIDNum
self.pointOutputFile = self.inventory.pointOutputFile
self.nodeValuesList = self.inventory.nodeValuesList
self.excludeZeroNormals = self.inventory.excludeZeroNormals
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
self.srcCoordSys = self.inventory.srcCoordSys
self.destCoordSys = self.inventory.destCoordSys
self.dataDim = self.inventory.dataDim
self.bcType = self.inventory.bcType
self.pointsFile = self.inventory.pointsFile
self.pointsSpatialDB = self.inventory.pointsSpatialDB
self.bcScale = self.inventory.bcScale
self.upDir = self.inventory.upDir
self.normalDir = self.inventory.normalDir
self.eulerLat = self.inventory.eulerLat
self.eulerLon = self.inventory.eulerLon
self.eulerRot = self.inventory.eulerRot
self.spaceDim = self.srcCoordSys.spaceDim
self.dipSlip = self.inventory.dipSlip
self.dipCutoff = self.inventory.dipCutoff
self.xMinVal = self.inventory.xMin.value
self.xMaxVal = self.inventory.xMax.value
self.yMinVal = self.inventory.yMin.value
self.yMaxVal = self.inventory.yMax.value
self.zMinVal = self.inventory.zMin.value
self.zMaxVal = self.inventory.zMax.value
self.defaultValues = self.inventory.defaultValues
lat = self.eulerLat.value
lon = self.eulerLon.value
rot = self.eulerRot.value
clat = math.cos(lat)
slat = math.sin(lat)
clon = math.cos(lon)
slon = math.sin(lon)
# Note that the Euler pole already includes the earth radius term.
self.eulerPole[0] = self.earthRad * rot * clat * clon
self.eulerPole[1] = self.earthRad * rot * clat * slon
self.eulerPole[2] = self.earthRad * rot * slat
self.upVec = numpy.array(
[float(self.upDir[0]),
float(self.upDir[1]),
float(self.upDir[2])],
dtype=numpy.float64)
self.normalVec = numpy.array([
float(self.normalDir[0]),
float(self.normalDir[1]),
float(self.normalDir[2])
],
dtype=numpy.float64)
self.dipCutoffProj = abs(math.sin(self.dipCutoff.value))
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
self.psetFile = self.inventory.psetFile
self.ucdFile = self.inventory.ucdFile
self.pointOutputFile = self.inventory.pointOutputFile
self.valuesList = self.inventory.valuesList
self.outputIndex = self.inventory.outputIndex
self.excludeZeroNormals = self.inventory.excludeZeroNormals
return
def __init__(self, name="make_synthetic_gpsdisp"):
Application.__init__(self, name)
self.coords = None
self.stations = None
self.dispRaw = None
self.dispNoise = None
self.numStations = 0
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
self.penaltyWeights = numpy.array(self.penaltyWeightVals, dtype=numpy.float64)
self.numPenaltyWeights = self.penaltyWeights.shape[0]
# Left-lateral and updip components from assumed rake.
self.llComp = math.cos(self.rake.value)
self.udComp = math.sin(self.rake.value)
return
def __init__(self, name="vtkdiff"):
Application.__init__(self, name)
self.vtkInputList = []
self.numVtkInputFiles = 0
self.vtkInputTimes = []
self.timeStampWidth = 0
self.numVertsPerCell = 0
self.numCells = 0
self.numVerts = 0
self.spaceDim = 0
self.cellType = None
self.readMesh = False
return
def __init__(self, name="dem2lines"):
Application.__init__(self, name)
self.numXIn = 0
self.numYIn = 0
self.numZIn = 0
self.numXOut = 0
self.numYOut = 0
self.numZOut = 0
self.xIn = None
self.yIn = None
self.zIn = None
self.xOut = None
self.yOut = None
self.zOut = None
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
# Filenames
self.faultInfoFile = self.inventory.faultInfoFile
self.faultResultsFile = self.inventory.faultResultsFile
self.faultOutputFile = self.inventory.faultOutputFile
# Parameters
self.frictionCoeff = self.inventory.frictionCoeff
self.stressScaleFactor = self.inventory.stressScaleFactor
self.slipScaleFactor = self.inventory.slipScaleFactor
self.shearDirection = self.inventory.shearDirection
return
def __init__(self, name="faultinfo"):
Application.__init__(self, name)
self.numVertsPerCell = 0
self.numCells = 0
self.cellsArray = numpy.array([0])
self.verticesArray = numpy.array([0])
self.normalVec = numpy.array([0])
self.strikeVec = numpy.array([0])
self.dipVec = numpy.array([0])
self.stressVec = numpy.array([0])
self.slipVec = numpy.array([0])
self.numVerts = 0
self.spaceDim = 0
self.cellType = ""
self.readMesh = False
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
import os
import re
# Set up info for input files
totalInputPath = os.path.normpath(os.path.join(os.getcwd(),
self.inventory.vtkInputRoot))
self.vtkInputDir = os.path.dirname(totalInputPath)
baseInputName = os.path.basename(totalInputPath)
baseInputNameLen = len(baseInputName)
if baseInputName.endswith(".vtk"):
baseInputNameStripped = baseInputName[0:baseInputNameLen - 4]
else:
baseInputNameStripped = baseInputName
testFind = re.search('_t[0-9]*$', baseInputNameStripped)
if testFind != None:
timeInd = baseInputNameStripped.rfind(testFind.group(0))
self.vtkInputRoot = baseInputNameStripped[0:timeInd]
else:
self.vtkInputRoot = baseInputNameStripped
# Set up info for output files
totalOutputPath = os.path.normpath(os.path.join(
os.getcwd(), self.inventory.vtkOutputRoot))
self.vtkOutputDir = os.path.dirname(totalOutputPath)
baseOutputName = os.path.basename(totalOutputPath)
baseOutputNameLen = len(baseOutputName)
if baseOutputName.endswith(".vtk"):
baseOutputNameStripped = baseOutputName[0:baseOutputNameLen - 4]
else:
baseOutputNameStripped = baseOutputName
testFind = re.search('_t[0-9]*$', baseOutputNameStripped)
if testFind != None:
timeInd = baseOutputNameStripped.rfind(testFind.group(0))
self.vtkOutputRoot = baseOutputNameStripped[0:timeInd]
else:
self.vtkOutputRoot = baseOutputNameStripped
self.scaleFactor = self.inventory.scaleFactor
return
def __init__(self, name="vtkcff"):
Application.__init__(self, name)
self.vtkInputList = []
self.numVtkInputFiles = 0
self.vtkInputTimes = []
self.timeStampWidth = 0
self.refFileIndex = 0
self.numVertsPerCell = 0
self.numCells = 0
self.cellsArray = numpy.array([0])
self.verticesArray = numpy.array([0])
self.numVerts = 0
self.spaceDim = 0
self.cellType = ""
self.readMesh = False
self.numStressPoints = 0
return
def __init__(self, name="princaxes"):
Application.__init__(self, name)
self.cells = None
self.vertArray = None
self.spaceDim = 0
self.cellType = ""
self.numTensorPoints = 0
self.tensorSorted = None
self.minPrincAxis = None
self.intPrincAxis = None
self.maxPrincAxis = None
self.minEigenvalue = None
self.intEigenvalue = None
self.maxEigenvalue = None
self.regionalField = numpy.zeros((3, 3), dtype=numpy.float64)
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
# Filenames
self.inputDem = self.inventory.inputDem
self.vtkOutputFile = self.inventory.vtkOutputFile
self.masterJournal = self.inventory.masterJournal
self.uLinePrefix = self.inventory.uLinePrefix
self.uLineJournal = self.inventory.uLineJournal
self.vLinePrefix = self.inventory.vLinePrefix
self.vLineJournal = self.inventory.vLineJournal
self.acisFilename = self.inventory.acisFilename
# Parameters
self.xMin = self.inventory.xMin
self.xMax = self.inventory.xMax
self.yMin = self.inventory.yMin
self.yMax = self.inventory.yMax
self.skipInterval = self.inventory.skipInterval
return
def __init__(self, name="geo2cart"):
Application.__init__(self, name)
self.numPoles = 0
self.polesSource = []
return
def __init__(self, name="surfaceapp"):
"""Constructor.
"""
Application.__init__(self, name)
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
import pdb
pdb.set_trace()
# Set up info for input files
totalInputPath = os.path.normpath(
os.path.join(os.getcwd(), self.inventory.vtkInputRoot))
self.vtkInputDir = os.path.dirname(totalInputPath)
baseInputName = os.path.basename(totalInputPath)
baseInputNameLen = len(baseInputName)
if baseInputName.endswith(".vtk"):
baseInputNameStripped = baseInputName[0:baseInputNameLen - 4]
else:
baseInputNameStripped = baseInputName
testFind = re.search('_t[0-9]*$', baseInputNameStripped)
if testFind != None:
timeInd = baseInputNameStripped.rfind(testFind.group(0))
self.vtkInputRoot = baseInputNameStripped[0:timeInd]
else:
self.vtkInputRoot = baseInputNameStripped
# Solution mode info
self.stressRefMode = self.inventory.stressRefMode
self.orientationMode = self.inventory.orientationMode
self.initialStateIndex = self.inventory.initialStateIndex
self.constantStateValues = self.inventory.constantStateValues
self.isotropicPoroelastic = self.inventory.isotropicPoroelastic
# Index information
self.vtkStressIndex = self.inventory.vtkStressIndex
self.vtkStressComponentsOrder = self.inventory.vtkStressComponentsOrder
# Parameters
self.frictionCoeff = self.inventory.frictionCoeff
self.skemptonCoeff = self.inventory.skemptonCoeff
self.cffPlaneNormal = numpy.array([
float(self.inventory.cffPlaneNormal[0]),
float(self.inventory.cffPlaneNormal[1]),
float(self.inventory.cffPlaneNormal[2]),
],
dtype=numpy.float64)
# Set up info for output files
totalOutputPath = os.path.normpath(
os.path.join(os.getcwd(), self.inventory.vtkOutputRoot))
self.vtkOutputDir = os.path.dirname(totalOutputPath)
baseOutputName = os.path.basename(totalOutputPath)
baseOutputNameLen = len(baseOutputName)
if baseOutputName.endswith(".vtk"):
baseOutputNameStripped = baseOutputName[0:baseOutputNameLen - 4]
else:
baseOutputNameStripped = baseOutputName
testFind = re.search('_t[0-9]*$', baseOutputNameStripped)
if testFind != None:
timeInd = baseOutputNameStripped.rfind(testFind.group(0))
self.vtkOutputRoot = baseOutputNameStripped[0:timeInd]
else:
self.vtkOutputRoot = baseOutputNameStripped
return
def _configure(self):
"""Setup members using inventory.
"""
Application._configure(self)
def __init__(self, name="grabpoints"):
Application.__init__(self, name)
self.numPoints = 0
self.indices = []
self.pointCoords = []
return
def __init__(self, name="grabfaces"):
Application.__init__(self, name)
return
