Python set_atypeの例

コード例 #1

ファイル: IBM-Unsteady-3D-Cantilever_MPMCoupling.py プロジェクト: btanasoi/fvm

#!/usr/bin/env python

import sys
import pdb

import fvm
import fvm.fvmbaseExt as fvmbaseExt
fvm.set_atype('double')
import fvm.importers as importers
import fvm.exporters_atyped_double as exporters
from mpi4py import MPI


import math
from FluidStructure import *

atype = 'double'
#atype = 'tangent'

from FluentCase import FluentCase

#fvmbaseExt.enableDebug("cdtor")

fileBase = None
numIterationsPerStep = 20
#mpm number of steps shoul numTimeSteps-1
numTimeSteps = 2 
fileBase = "/home/yildirim/memosa/src/fvm/test/"

vFile = open(fileBase + "velocity-fullbeam-test.out","w")
pFile = open(fileBase + "pIntegral-fullbeam-new.out","w")

コード例 #2

ファイル: parallel_cylinder2D.py プロジェクト: wuhao1938/fvm

#!/usr/bin/env python
### import  modules ###
"""
mpirun -np   1 python ./parallel_cylinder.py  ../test/uniform-cart_5K.cas ../test/circle_66.cas --type=quad
"""
import pdb
import sys
from math import *
import fvm
fvm.set_atype('double')
#import tecplotExporter
import fvm.fvmbaseExt as fvmbaseExt
import fvm.importers as importers
import fvm.models_atyped_double as models
import fvm.exporters_atyped_double as exporters
import fvm.fvmparallel as fvmparallel
from FluentCase import FluentCase
from fvm.fvmbaseExt import VecD3
import time
from mpi4py import MPI
import optparse
from Tools import *
from ComputeForce import *
from TimeStep import *
from MeshSetup import *
from tecplotParallelFlowField import *
from mpi4py import MPI

etype = {'tri': 1, 'quad': 2, 'tetra': 3, 'hexa': 4}

コード例 #3

ファイル: FvmTestFlowModel.py プロジェクト: btanasoi/fvm

if len(sys.argv) < 2:
    usage()

fileBase = sys.argv[1]
if len(sys.argv) == 3:
    outfile = sys.argv[2]
else:
    outfile = fileBase+".dat"

(options, args) = parser.parse_args()

if options.tangent:
    atype = 'tangent'
else:
    atype = 'double'
fvm.set_atype(atype)
exporters = fvm.exporters

#fvmbaseExt.enableDebug("cdtor")

print "reading",fileBase+".cas"
reader = FluentCase(fileBase+".cas")
reader.read();
meshes = reader.getMeshList()

t0 = time.time()
geomFields =  fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields,meshes)
metricsCalculator.init()

if atype == 'tangent':

コード例 #4

ファイル: testCant_Torder.py プロジェクト: drm42/fvm-drm

#!/usr/bin/env python

import sys

sys.setdlopenflags(0x100 | 0x2)

import fvm
import fvm.fvmbaseExt as fvmbaseExt
import fvm.importers as importers

fvm.set_atype("double")

import math

if fvm.atype == "double":
    import fvm.models_atyped_double as models
    import fvm.exporters_atyped_double as exporters
elif fvm.atype == "tangent":
    import fvm.models_atyped_tangent_double as models
    import fvm.exporters_atyped_tangent_double as exporters

from FluentCase import FluentCase
from optparse import OptionParser

# fvmbaseExt.enableDebug("cdtor")

fileBase = None
fileBase = "/scratch/prism/shankha/memosa/src/fvm/test/CANT-TORDER/Co_0.05/cbeam7"
# fileBase = "/home/sm/app-memosa/src/fvm/test/cav32"
# fileBase = "/home/sm/a/data/wj"

コード例 #5

ファイル: FvmTestFlowModel.py プロジェクト: wuhao1938/fvm

if len(sys.argv) < 2:
    usage()

fileBase = sys.argv[1]
if len(sys.argv) == 3:
    outfile = sys.argv[2]
else:
    outfile = fileBase + ".dat"

(options, args) = parser.parse_args()

if options.tangent:
    atype = 'tangent'
else:
    atype = 'double'
fvm.set_atype(atype)
exporters = fvm.exporters

#fvmbaseExt.enableDebug("cdtor")

print "reading", fileBase + ".cas"
reader = FluentCase(fileBase + ".cas")
reader.read()
meshes = reader.getMeshList()

t0 = time.time()
geomFields = fvm.models.GeomFields('geom')
metricsCalculator = fvm.models.MeshMetricsCalculatorA(geomFields, meshes)
metricsCalculator.init()

if atype == 'tangent':