#!/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()