#!/usr/bin/python
import Manifold
import ManifoldVTK
import Utils
import operator
import sys
import numpy
coords = Utils.read_point_file("./datasets/swissroll.txt")
embedded_coords, mapping = Manifold.do_embedding(coords, k=10)
coords = coords.take(mapping, axis=0)
nb_points = len(coords)
# sorting
sorted_index = list(index
for index, item in sorted(enumerate(embedded_coords[:, 0]),
key=operator.itemgetter(1)))
colors = [(255, 255, 255) for i in range(0, nb_points)]
for i in range(0, nb_points):
if i < nb_points / 2:
f = float(i) / float(nb_points / 2)
color = (0, f * 255, 255 - f * 255)
#color = ( 0, 0, 255 )
else:
f = float(i - nb_points / 2) / float(nb_points / 2)
color = (f * 255, 255 - f * 255, 0)
#color = ( 255, 0, 0 )
#!/usr/bin/python
import Manifold
import ManifoldVTK
import Utils
import operator
import sys
import numpy
coords = Utils.read_point_file("./datasets/swissroll.txt")
embedded_coords, mapping = Manifold.do_embedding(coords,k=10)
coords = coords.take(mapping,axis=0)
nb_points = len(coords)
# sorting
sorted_index = list( index for index, item in sorted( enumerate( embedded_coords[:,0] ),
key=operator.itemgetter(1)
) )
colors = [(255,255,255) for i in range(0,nb_points)]
for i in range(0,nb_points):
if i < nb_points / 2:
f = float(i) / float(nb_points / 2)
color = ( 0, f*255, 255 - f*255 )
#color = ( 0, 0, 255 )
else:
f = float(i - nb_points / 2) / float(nb_points / 2)
color = ( f*255, 255 - f*255, 0 )
#color = ( 255, 0, 0 )
#All the well heads are divided into three list for each function
listOfProductionWellHeads = Rojas.listOfWellHeads.getList('Production')
listOfWaterInjectionWellHeads = Rojas.listOfWellHeads.getList('Water Injection')
listOfGasInjectionWellHeads = Rojas.listOfWellHeads.getList('Gas Injection')
#Determines the listOfPoints for every manifold
listOfProductionManifoldPoints = Rojas.betterKmeans(listOfProductionWellHeads,nummberOfManifoldForProductionWH)
listOfWaterInjectionManifoldPoints = Rojas.betterKmeans(listOfWaterInjectionWellHeads,nummberOfManifoldForWaterInjectionWH)
listOfGasInjectionManifoldPoints = Rojas.betterKmeans(listOfGasInjectionWellHeads,nummberOfManifoldForGasInjectionWH)
#Since we have the well heads for every manifold it is posible to INITIALIZE each manifold
listOfProductionManifold=list()
for i in listOfProductionManifoldPoints:
listOfProductionManifold.append(Manifold(i))
listOfWaterInjectionManifold=list()
for i in listOfWaterInjectionManifoldPoints:
listOfWaterInjectionManifold.append(Manifold(i))
listOfGasInjectionManifold=list()
for i in listOfGasInjectionManifoldPoints:
listOfGasInjectionManifold.append(Manifold(i))
#Calculate the total sum of gathering lines:
totalGatheringLength = 0
for i in listOfProductionManifold:
totalGatheringLength += i.lengthOfGathLines
#!/usr/bin/python
import Manifold
import ManifoldPIL
import Utils
from glob import glob
import argparse
parser = argparse.ArgumentParser(usage="usage: %(prog)s directory")
parser.add_argument('image_directory',
metavar='image_directory',
help='directory of images to process')
args = parser.parse_args()
filenames = glob(args.image_directory + '/*')
coords = Utils.read_images(filenames)
embedded_coords, mapping = Manifold.do_embedding(coords, tree='spilltree')
ManifoldPIL.render2D([filenames[i] for i in mapping], embedded_coords)
#!/usr/bin/python
import Manifold
import ManifoldPIL
import Utils
from glob import glob
import argparse
parser = argparse.ArgumentParser(
usage = "usage: %(prog)s directory" )
parser.add_argument(
'image_directory',
metavar='image_directory',
help='directory of images to process' )
args = parser.parse_args()
filenames = glob(args.image_directory + '/*')
coords = Utils.read_images(filenames)
embedded_coords, mapping = Manifold.do_embedding(coords,tree='spilltree')
ManifoldPIL.render2D([ filenames[i] for i in mapping], embedded_coords)
# -*- coding: utf-8 -*-
import sys, os
sys.path.append('PressureDropAndHeatTransfer')
sys.path.append('../REFPROP')
os.environ['RPPREFIX'] = r'C:/Program Files (x86)/REFPROP'
from Manifold import *
if __name__ == '__main__':
path = "../XML/QuarterShell/"
rootXML = "endcapShortyRing.xml"
m = Manifold(path, rootXML, 1.1e-3, 0.01, 2., -40.)
# m.branches[0].setFinalVaporQualityGuess(0.5) # very advanced, don't do this
m.run()
m.plot()
# -*- coding: utf-8 -*-
import sys, os
sys.path.append('PressureDropAndHeatTransfer')
sys.path.append('../REFPROP')
os.environ['RPPREFIX'] = r'C:/Program Files (x86)/REFPROP'
from Manifold import *
if __name__ == '__main__':
path = "../XML/Manifoldv1/"
rootXML = "manifold0.xml"
MF = 5e-3
m = Manifold(path, rootXML, 2.02e-3, 0.01, -35.)
m.run()
m.plot()
"""
Fig. 1
__________
______/ b3 \_____
/ \__________/ \
/ b1 b4 b5 \
______/ \______
b0 \ / b6
\ /
\___________b2__________/
Run b5, guess initial vapor quality, gives vapor quality that we guessed and Temperature
Fig. 2
____b0__________b2__________b4_________b6_____
# -*- coding: utf-8 -*-
import sys, os
sys.path.append('../FluidDynamics')
sys.path.append('../REFPROP')
os.environ['RPPREFIX'] = r'C:/Program Files (x86)/REFPROP'
from Manifold import *
import runinfo
if __name__ == '__main__':
if len(sys.argv) > 1:
runinfo.runname = sys.argv[1]
else:
runinfo.runname = 'test'
path = "../XML/HalfBarrel/"
rootXML = "barrel.xml"
m = Manifold(path, rootXML, 2.*2.304e-3+2.*7.68e-3, 0.0, -1, -40.)
# m = Manifold(path, rootXML, 2.*2.304e-3, 0.0, -1, -40.)
m.run()
m.plot()
