def test_ReadEnsightGeo(self):
parts = [1, 2]
partRec = []
THIS_FOLDER = os.path.dirname(os.path.abspath(__file__))
my_file = os.path.join(THIS_FOLDER, 'testEnsightReader.geo')
mesh = EnsightReader2.readEnsightGeo(my_file, parts, partRec)
# mesh = EnsightReader2.readEnsightGeo('MMPTestCase_v1.geo', parts, partRec)
timeUnits = PQ.PhysicalUnit('s', 1., [0, 0, 1, 0, 0, 0, 0, 0, 0])
temperatureUnits = PQ.PhysicalUnit('K', 1.,
[0, 0, 0, 0, 1, 0, 0, 0, 0])
# f = EnsightReader2.readEnsightField('fld_TEMPERATURE.escl', parts, partRec, 1, FieldID.FID_Temperature, mesh, temperatureUnits, 0)
self.assertTrue(1 == 1)
def test_cmp(self):
timeUnits = PQ.PhysicalUnit('s', 1., [0,0,1,0,0,0,0,0,0])
timeUnits2 = PQ.PhysicalUnit('s', 2., [0, 0, 1, 0, 0, 0, 0, 0, 0])
timeUnits3 = PQ.PhysicalUnit('s', -1., [0, 0, 1, 0, 0, 0, 0, 0, 0])
timeUnits4 = PQ.PhysicalUnit('s', 1., [0, 0, 1, 0, 0, 0, 0, 0, 0])
temperatureUnits = PQ.PhysicalUnit('K', 1., [0,0,0,0,1,0,0,0,0])
try:
a = timeUnits.__cmp__(temperatureUnits)
except TypeError:
pass
except Exception as e:
self.fail('Unexpected exception raised:', e)
else:
self.fail('Exception not raised')
self.assertTrue(timeUnits.__cmp__(timeUnits2) < 0)
self.assertTrue(timeUnits.__cmp__(timeUnits3) > 0)
self.assertTrue(timeUnits.__cmp__(timeUnits4) == 0)
self.assertTrue(timeUnits.__cmp__(5) == -1)
#!/usr/bin/env python
from __future__ import print_function
from builtins import str
import sys
sys.path.append('../../..')
from mupif import *
import logging
log = logging.getLogger()
import mupif.Physics.PhysicalQuantities as PQ
timeUnits = PQ.PhysicalUnit('s', 1., [0,0,1,0,0,0,0,0,0])
class application1(Application.Application):
"""
Simple application that generates a property with a value equal to actual time
"""
def __init__(self, file):
super(application1,self).__init__(self,file)
def getProperty(self, propID, time, objectID=0):
if (propID == PropertyID.PID_Velocity):
return Property.ConstantProperty(self.value, PropertyID.PID_Velocity, ValueType.Scalar, 'm/s', time, 0)
else:
raise APIError.APIError ('Unknown property ID')
def solveStep(self, tstep, stageID=0, runInBackground=False):
time = tstep.getTime().inUnitsOf(timeUnits).getValue()
self.value=1.0*time
def getCriticalTimeStep(self):
import sys
sys.path.append('../../..')
from mupif import *
import mupif.Physics.PhysicalQuantities as PQ
timeUnits = PQ.PhysicalUnit('s', 1., [0, 0, 1, 0, 0, 0, 0, 0, 0])
temperatureUnits = PQ.PhysicalUnit('K', 1., [0, 0, 0, 0, 1, 0, 0, 0, 0])
class Celsian(Application.Application):
def __init__(self, file):
super(Celsian, self).__init__(file)
self.mesh = None
def getField(self, fieldID, time):
parts = [1, 2]
partRec = []
if (self.mesh == None):
self.mesh = EnsightReader2.readEnsightGeo(
'paraview/MMPTestCase_v1.geo', parts, partRec)
f = EnsightReader2.readEnsightField('paraview/fld_TEMPERATURE.escl',
parts, partRec, 1,
FieldID.FID_Temperature, self.mesh,
temperatureUnits, time)
return f
def solveStep(self, tstep, stageID=0, runInBackground=False):
time = tstep.getTime()
self.value = 1.0 * time
#!/usr/bin/env python
from builtins import range
#
# This example requires pyvtk module, install it using
# pip install pyvtk
# Note: pyvtk is only available in Python 2.x (Sept, 2014)
#
import sys
import logging
sys.path.append('../../..')
from mupif import *
import mupif.Physics.PhysicalQuantities as PQ
log = logging.getLogger()
temperatureUnit = PQ.PhysicalUnit('K', 1., [0, 0, 0, 0, 1, 0, 0, 0, 0])
def main():
# generate field and corresponding mesh
mesh = Mesh.UnstructuredMesh()
vertices = []
values1 = []
values2 = []
num = 0
for i in range(40):
for j in range(15):
vertices.append(
Vertex.Vertex((i * 15) + j, (i * 15) + j + 1,
(float(i), float(j), 0.0)))
values1.append((num, ))
import sys
sys.path.append('../../..') # for mupif import
import os
import platform
from subprocess import *
import socket
import time as timeTime
import pyvtk
from mupif import *
import logging
import mupif.Physics.PhysicalQuantities as PQ
dummyUnits = PQ.PhysicalUnit('dummy', 1., [0, 0, 0, 0, 0, 0, 0, 0, 0])
timeUnits = PQ.PhysicalUnit('s', 1., [0, 0, 1, 0, 0, 0, 0, 0, 0])
import micressConfig as mConf
#Micress import
#import Property
import MICPropertyID
import MICFieldID
VtkReader2.pyvtk_monkeypatch()
class micress(Application.Application):
"""
MICRESS application interface class
