summaryrefslogtreecommitdiff
path: root/src/main/python/utils/Bin_Phase_env.py
diff options
context:
space:
mode:
authorbrenda-br2023-03-04 11:32:15 +0530
committerbrenda-br2023-03-04 11:32:15 +0530
commit3cbdd4238867bc860282f7cf702d73b5be6e3f86 (patch)
tree7a9323daf6923801a61bb10b5aa35ad4795c52b1 /src/main/python/utils/Bin_Phase_env.py
parent7cf8eca28f9ca77e6195aa43e3bd9ae319518b66 (diff)
downloadChemical-Simulator-GUI-3cbdd4238867bc860282f7cf702d73b5be6e3f86.tar.gz
Chemical-Simulator-GUI-3cbdd4238867bc860282f7cf702d73b5be6e3f86.tar.bz2
Chemical-Simulator-GUI-3cbdd4238867bc860282f7cf702d73b5be6e3f86.zip
Revert "Restructuring Finalized for App Bundling"
This reverts commit 7af3526e105cc330422f8742ec5edec1c4a0a98f.
Diffstat (limited to 'src/main/python/utils/Bin_Phase_env.py')
-rw-r--r--src/main/python/utils/Bin_Phase_env.py241
1 files changed, 241 insertions, 0 deletions
diff --git a/src/main/python/utils/Bin_Phase_env.py b/src/main/python/utils/Bin_Phase_env.py
new file mode 100644
index 0000000..dd6f938
--- /dev/null
+++ b/src/main/python/utils/Bin_Phase_env.py
@@ -0,0 +1,241 @@
+import sys
+import pandas as pd
+import numpy as np
+import os
+import csv
+from subprocess import Popen, PIPE
+
+from PyQt5.QtCore import *
+from PyQt5.QtWidgets import *
+from PyQt5.QtGui import *
+import PyQt5.QtGui as QtGui
+import PyQt5.QtCore as QtCore
+import PyQt5.QtWidgets as QtWidgets
+from PyQt5.uic import loadUiType
+
+
+import pyqtgraph as pg
+import pyqtgraph.exporters
+
+current = os.path.dirname(os.path.realpath(__file__))
+parent = os.path.dirname(current)
+parentPath = os.path.dirname(parent)
+sys.path.append(parentPath)
+
+ui_dialog,_ = loadUiType(parentPath+'/ui/utils/Binary_Phase_Env.ui')
+pg.CONFIG_OPTIONS['crashWarning'] = False
+
+class BinPhaseEnv(QWidget,ui_dialog):
+ def __init__(self,comp):
+ QWidget.__init__(self)
+ self.setupUi(self)
+
+ self.comp = comp
+ self.compunds = self.comp.get_compounds()
+
+ for i in self.compunds:
+ self.comboBox.addItem(str(i))
+
+ for i in range(len(self.compunds)):
+ if i!=0:
+ self.comboBox_2.addItem(str(self.compunds[i]))
+
+ self.comboBox_2.addItem(str(self.compunds[0]))
+
+ self.lines = [line.rstrip('\n') for line in open(parentPath+'/python/utils/thermopackage.txt')]
+ for j in self.lines:
+ self.comboBox_3.addItem(str(j))
+
+ self.radioButton.setChecked(False)
+ self.radioButton_2.setChecked(False)
+
+ self.button_handler()
+ self.counter = 1
+
+ def button_handler(self):
+ self.radioButton.clicked.connect(lambda:self.T_xy())
+ self.radioButton_2.clicked.connect(lambda:self.P_xy())
+ self.pushButton_2.clicked.connect(lambda:self.plot())
+
+ def T_xy(self):
+ self.radioButton.setChecked(True)
+ self.radioButton_2.setChecked(False)
+ for i in reversed(range(self.formLayout.count())):
+ self.formLayout.itemAt(i).widget().setParent(None)
+
+ self.first = QLineEdit()
+ self.type = "P"
+ self.other = "T"
+ self.otherunit = "(K)"
+ self.formLayout.addRow(QLabel("P(Pa)"),self.first)
+ self.points = QLineEdit()
+ self.points.setText("40")
+ self.formLayout.addRow(QLabel("Number of data points"),self.points)
+
+ def P_xy(self):
+ self.radioButton_2.setChecked(True)
+ self.radioButton.setChecked(False)
+ for i in reversed(range(self.formLayout.count())):
+ self.formLayout.itemAt(i).widget().setParent(None)
+
+ self.first = QLineEdit()
+ self.type = "T"
+ self.other = "P"
+ self.otherunit = "(Pa)"
+ self.points = QLineEdit()
+ self.points.setText("40")
+ self.formLayout.addRow(QLabel("T(K)"),self.first)
+ self.formLayout.addRow(QLabel("Number of data points"),self.points)
+
+ def get_omc_path(self):
+ try:
+ self.omhome = os.environ.get('OPENMODELICAHOME')
+ if self.omhome is None:
+ self.omhome = os.path.split(os.path.split(os.path.realpath(spawn.find_executable("omc")))[0])[0]
+ elif os.path.exists('/opt/local/bin/omc'):
+ self.omhome = '/opt/local'
+ elif os.path.exists('/usr/bin/omc'):
+ self.omhome = '/usr'
+ return os.path.join(self.omhome, 'bin', 'omc')
+ except BaseException:
+ #print("The OpenModelica compiler is missing in the System path please install it" )
+ raise
+
+ def plot(self):
+ try:
+ val = int(self.first.text(),10)
+ except:
+ val = 0
+ try:
+ data_points = int(self.points.text(),10)
+ except:
+ data_points = 0
+
+ self.curr_path = os.getcwd()
+ self.sim_dir_path = os.path.join(self.curr_path,'./../Simulator')
+ self.Graphmo_path = os.path.join(self.sim_dir_path,'Graph.mo')
+ self.plot_mos_path = os.path.join(self.sim_dir_path,'PlotGraph.mos')
+
+ self.data = []
+
+ self.comp1 = self.comboBox.currentText()
+ self.comp2 = self.comboBox_2.currentText()
+ self.comp_1 = self.comboBox.currentText().split('(')[0]
+ self.comp_2 = self.comboBox_2.currentText().split('(')[0]
+
+ self.thermoPack = self.comboBox_3.currentText()
+
+ self.data.append("model Graph\n")
+ self.data.append("import data = Simulator.Files.ChemsepDatabase;\n")
+ self.data.append("parameter data."+self.comp_1+" comp1;\n")
+ self.data.append("parameter data."+self.comp_2+" comp2;\n")
+ self.data.append("extends BinaryEnvelopes."+self.thermoPack+"(Nc = 2, data_points = "+str(data_points)+ ", comp = { comp1, comp2 }, "+self.type+" = fill( "+str(val)+", "+str(data_points)+"));\n")
+ self.data.append("end Graph;")
+
+ with open(self.Graphmo_path, 'w') as txtfile:
+ for d in self.data:
+ txtfile.write(str(d))
+
+ with open(self.plot_mos_path, 'w') as mosFile:
+ mosFile.write("loadModel(Modelica);\n")
+ mosFile.write("loadFile(\"Simulator/package.mo\");\n")
+ mosFile.write("loadFile(\"BinaryEnvelopes.mo\");\n")
+ mosFile.write("loadFile(\"Graph.mo\");\n")
+ mosFile.write("simulate(Graph, outputFormat=\"csv\", stopTime=1.0, numberOfIntervals=1);\n")
+
+ self.resdata = []
+ self.omc_path = self.get_omc_path()
+ simpath = self.plot_mos_path
+ os.chdir(self.sim_dir_path)
+
+ process = Popen([self.omc_path, '-s',simpath], stdout=PIPE, stderr=PIPE)
+ self.stdout, self.stderr = process.communicate()
+
+ os.chdir(self.curr_path)
+
+ csvpath = os.path.join(self.sim_dir_path,'Graph_res.csv')
+
+ self.datay = []
+ self.datax1 = []
+ self.datax2 = []
+ self.rows = []
+
+ with open (csvpath,'r') as resultFile:
+ self.resdata = []
+ csvreader = csv.reader(resultFile,delimiter=',')
+ for row in csvreader:
+ self.resdata.append(row)
+ self.rows.append(row)
+
+ if self.type=='T':
+ for k in range(len(self.rows[0])):
+ if self.rows[0][k][0]=='P':
+ self.datay.append(float(self.rows[1][k]))
+ length = len(self.rows[0][k])
+ if self.rows[0][k][0]=='x' and self.rows[0][k][length-2]=='1':
+ self.datax1.append(float(self.rows[1][k]))
+ if self.rows[0][k][0]=='y' and self.rows[0][k][length-2]=='1':
+ self.datax2.append(float(self.rows[1][k]))
+ else:
+ for k in range(len(self.rows[0])):
+ if self.rows[0][k][0]=='T':
+ self.datay.append(float(self.rows[1][k]))
+ length = len(self.rows[0][k])
+ if self.rows[0][k][0]=='x' and self.rows[0][k][length-2]=='1':
+ self.datax1.append(float(self.rows[1][k]))
+
+ if self.rows[0][k][0]=='y' and self.rows[0][k][length-2]=='1':
+ self.datax2.append(float(self.rows[1][k]))
+
+ plt = pg.PlotWidget()
+ plt.showGrid(x=True,y=True)
+ plt.addLegend()
+ plt.setXRange(0,1)
+
+ c1 = plt.plot(self.datax1, self.datay,pen=pg.mkPen('b',width = 1), name='dew points')
+ c2 = plt.plot(self.datax2, self.datay,pen=pg.mkPen('r',width = 1), name='bubble points')
+ view_box = plt.plotItem.vb
+ self.tool_tip = ""
+
+ def press_event(evt):
+ a = 10
+ pos = evt
+ mousepoint = view_box.mapSceneToView(pos)
+ roi = pg.ROI(pos)
+ find_color = plt.mapToGlobal(pos.toPoint())
+
+ screen = QGuiApplication.primaryScreen()
+ image = screen.grabWindow(QApplication.desktop().winId()).toImage()
+ colour = QtGui.QColor(image.pixel(find_color.x(),find_color.y()))
+
+ if colour.red()==255 or colour.blue()==255:
+ self.lineEdit_x.setText(str(round(mousepoint.x(),3)))
+ self.lineEdit_y.setText(str(round(mousepoint.y(),3)))
+ self.tool_tip = str(round(mousepoint.x(),3)) + ", " + str(round(mousepoint.y(),3))
+ QApplication.setOverrideCursor(QCursor(QtCore.Qt.CrossCursor))
+ else:
+ self.lineEdit_x.setText("")
+ self.lineEdit_y.setText("")
+ self.tool_tip = ""
+ QApplication.setOverrideCursor(QCursor(QtCore.Qt.ArrowCursor))
+
+ def entered(items):
+ for i in items:
+ if i.__class__.__name__ =="LegendItem":
+ self.lineEdit_x.setText("")
+ self.lineEdit_y.setText("")
+ QApplication.setOverrideCursor(QCursor(QtCore.Qt.ArrowCursor))
+ else:
+ i.setToolTip(self.tool_tip)
+
+ plt.scene().sigMouseMoved.connect(press_event)
+ plt.scene().sigMouseHover.connect(entered)
+
+ plt.setLabel('left',self.other+self.otherunit,units = '')
+ plt.setLabel('bottom',self.comp1+'(mol. frac.)',units = '')
+
+ self.new_tab = plt
+ self.new_tab.setObjectName("Plot "+str(self.counter))
+
+ self.tabWidget.addTab(self.new_tab,"Plot "+str(self.counter))
+ self.counter+=1