working for python 3

Author: tcaduser

2dresult.pdf

@@ -61,7 +61,7 @@ def setup_atox(device, region):
         interface_model(device=device, interface=i, name="iindex", equation="node_index@r1")
       nlist = get_interface_model_values(device=device, interface=i, name="iindex")
       for j in nlist:
-        print j
+        print(j)
         set_node_value(device=device, region=region, index=int(j), name="AtOx", value=1.0)
 
 def setup_dg_variable(device, region, variable):

Electrons_1.5nm.pdf

@@ -31,7 +31,7 @@ def SetUniversalParameters(device, region):
     'k' :          1.3806503e-23,    #, 'J/K'),
     'Permittivity_0' :  8.85e-14     #, 'F/cm^2')
   }
-  for k, v in universal.items():
+  for k, v in list(universal.items()):
     set_parameter(device=device, region=region, name=k, value=v)
 
 def SetOxideParameters(device, region):
@@ -93,7 +93,7 @@ def SetSiliconParameters(device, region):
 #    "T" : 300
   }
 
-  for k, v in par.items():
+  for k, v in list(par.items()):
     set_parameter(device=device, region=region, name=k, value=v)
 
 def CreateQuasiFermiLevels(device, region, electron_model, hole_model, variables):
@@ -176,7 +176,7 @@ def CreateBandEdgeModels(device, region, variables):
     #('EPdiff', 'EVdiff + Lh@n0 - Lh@n1', ('Potential', 'T', 'Lh')),
   )
   for (model, equation, variable_list) in eeq:
-    print model
+    print(model)
     CreateEdgeModel(device, region, model, equation)
     vset = set(variable_list)
     for v in variables:
@@ -581,7 +581,7 @@ def CreateOxidePotentialOnly(device, region, update_type="default"):
     Creates Potential solution variable if not available
   '''
   if not InNodeModelList(device, region, "Potential"):
-    print "Creating Node Solution Potential"
+    print("Creating Node Solution Potential")
     CreateSolution(device, region, "Potential")
 
   efield="(Potential@n0 - Potential@n1)*EdgeInverseLength"

Electrons_3nm.pdf

@@ -28,7 +28,7 @@ def CreateNodeModel(device, region, model, expression):
   '''
   result=node_model(device=device, region=region, name=model, equation=expression)
   if debug:
-    print("NODEMODEL {d} {r} {m} \"{re}\"".format(d=device, r=region, m=model, re=result))
+    print(("NODEMODEL {d} {r} {m} \"{re}\"".format(d=device, r=region, m=model, re=result)))
 
 def CreateNodeModelDerivative(device, region, model, expression, *vars):
   '''
@@ -47,7 +47,7 @@ def CreateContactNodeModel(device, contact, model, expression):
   '''
   result=contact_node_model(device=device, contact=contact, name=model, equation=expression)
   if debug:
-    print("CONTACTNODEMODEL {d} {c} {m} \"{re}\"".format(d=device, c=contact, m=model, re=result))
+    print(("CONTACTNODEMODEL {d} {c} {m} \"{re}\"".format(d=device, c=contact, m=model, re=result)))
 
 
 def CreateContactNodeModelDerivative(device, contact, model, expression, variable):
@@ -65,7 +65,7 @@ def CreateEdgeModel (device, region, model, expression):
   '''
   result=edge_model(device=device, region=region, name=model, equation=expression)
   if debug:
-    print "EDGEMODEL {d} {r} {m} \"{re}\"".format(d=device, r=region, m=model, re=result);
+    print("EDGEMODEL {d} {r} {m} \"{re}\"".format(d=device, r=region, m=model, re=result));
 
 def CreateEdgeModelDerivatives(device, region, model, expression, variable):
   '''
@@ -86,7 +86,7 @@ def CreateContactEdgeModel(device, contact, model, expression):
   '''
   result=contact_edge_model(device=device, contact=contact, name=model, equation=expression)
   if debug:
-    print("CONTACTEDGEMODEL {d} {c} {m} \"{re}\"".format(d=device, c=contact, m=model, re=result))
+    print(("CONTACTEDGEMODEL {d} {c} {m} \"{re}\"".format(d=device, c=contact, m=model, re=result)))
 
 def CreateContactEdgeModelDerivative(device, contact, model, expression, variable):
   '''
@@ -101,7 +101,7 @@ def CreateInterfaceModel(device, interface, model, expression):
   '''
   result=interface_model(device=device, interface=interface, name=model, equation=expression)
   if debug:
-    print("INTERFACEMODEL {d} {i} {m} \"{re}\"".format(d=device, i=interface, m=model, re=result))
+    print(("INTERFACEMODEL {d} {i} {m} \"{re}\"".format(d=device, i=interface, m=model, re=result)))
 
 #def CreateInterfaceModelDerivative(device, interface, model, expression, variable):
 #  '''
@@ -138,19 +138,19 @@ def EnsureEdgeFromNodeModelExists(device, region, nodemodel):
     Checks if the edge models exists
   '''
   if not InNodeModelList(device, region, nodemodel):
-    raise("{} must exist" % (nodemodel))
+    raise "{} must exist"
 
   emlist = get_edge_model_list(device=device, region=region)
   emtest = ("{0}@n0".format(nodemodel) and "{0}@n1".format(nodemodel))
   if not emtest:
     if debug:
-      print "INFO: Creating ${0}@n0 and ${0}@n1".format(nodemodel)
+      print("INFO: Creating ${0}@n0 and ${0}@n1".format(nodemodel))
     edge_from_node_model(device=device, region=region, node_model=nodemodel)
 
 def CreateElementModel2d(device, region, model, expression):
   result=element_model(device=device, region=region, name=model, equation=expression)
   if debug:
-    print("ELEMENTMODEL {d} {r} {m} \"{re}\"".format(d=device, r=region, m=model, re=result))
+    print(("ELEMENTMODEL {d} {r} {m} \"{re}\"".format(d=device, r=region, m=model, re=result)))
 
 
 def CreateElementModelDerivative2d(device, region, model_name, expression, *args):

Holes_1.5nm.pdf

@@ -28,7 +28,7 @@ def rampparam(device, region, param, stop, step_size, min_step, solver_params=No
     step_sign=-1
   last_param=start_param
   while(abs(last_param - stop) > min_step):
-    print("last end %e %e") % (last_param, stop)
+    print("last end %e %e" % (last_param, stop))
     next_param=last_param + step_sign * step_size
     if next_param < stop:
       next_step_sign=1
@@ -37,17 +37,17 @@ def rampparam(device, region, param, stop, step_size, min_step, solver_params=No
 
     if next_step_sign != step_sign:
       next_param=stop
-      print "setting to last param %e" % (stop)
-      print "setting next param %e" % (next_param)
+      print("setting to last param %e" % (stop))
+      print("setting next param %e" % (next_param))
     set_parameter(device=device, region=region, name=param, value=next_param)
     try:
       solve(**solver_params)
     except devsim.error as msg:
-      if msg[0].find("Convergence failure") != 0:
+      if str(msg).find("Convergence failure") != 0:
         raise
       set_parameter(device=device, region=region, name=param, value=last_param)
       step_size *= 0.5
-      print "setting new step size %e" % (step_size)
+      print("setting new step size %e" % (step_size))
       #raise NameError("STOP")
       if step_size < min_step:
         raise RuntimeError("Min step size too small")
@@ -67,7 +67,7 @@ def rampbias(device, contact, end_bias, step_size, min_step, max_iter, rel_error
     step_sign=-1
   last_bias=start_bias
   while(abs(last_bias - end_bias) > min_step):
-    print("last end %e %e") % (last_bias, end_bias)
+    print(("last end %e %e") % (last_bias, end_bias))
     next_bias=last_bias + step_sign * step_size
     if next_bias < end_bias:
       next_step_sign=1
@@ -76,8 +76,8 @@ def rampbias(device, contact, end_bias, step_size, min_step, max_iter, rel_error
 
     if next_step_sign != step_sign:
       next_bias=end_bias
-      print "setting to last bias %e" % (end_bias)
-      print "setting next bias %e" % (next_bias)
+      print("setting to last bias %e" % (end_bias))
+      print("setting next bias %e" % (next_bias))
     set_parameter(device=device, name=GetContactBiasName(contact), value=next_bias)
     try:
       solve(type="dc", absolute_error=abs_error, relative_error=rel_error, maximum_iterations=max_iter)
@@ -86,11 +86,11 @@ def rampbias(device, contact, end_bias, step_size, min_step, max_iter, rel_error
         raise
       set_parameter(device=device, name=GetContactBiasName(contact), value=last_bias)
       step_size *= 0.5
-      print "setting new step size %e" % (step_size)
+      print("setting new step size %e" % (step_size))
       if step_size < min_step:
-        raise "Min step size too small"
+        raise RuntimeError("Min step size too small")
       continue
-    print "Succeeded"
+    print("Succeeded")
     last_bias=next_bias
     callback(device)
 

Holes_3nm.pdf

@@ -1,5 +1,5 @@
 #!/bin/bash
-# use python 2.7 environment
+# use python 3.7 environment
 set -e
 for i in '1.5' '3'; do \
   for j in 'Electrons' 'Holes'; do \

dg_common.py

@@ -207,11 +207,11 @@ def setup_potential_only():
   carrier_var = sys.argv[2]
   pdfname = sys.argv[3]
 except Exception as e:
-  print '''args tox carrier pdfname
+  print('''args tox carrier pdfname
 tox:     oxide thickness (nm)
 carrier: Electrons or Holes
 pdfname: pdf file name
-'''
+''')
   raise
 
 
@@ -349,7 +349,7 @@ def mycb():
   set_node_values(device=device, region=region_ox, name="Potential", init_from="zero")
 
   set_parameter(name=GetContactBiasName("top"), value=0.0)
-  print d
+  print(d)
   set_parameter(device=device, region=region_si, name="Gamman", value=10)
   set_parameter(device=device, region=region_si, name="Gammanox", value=1)
   set_parameter(device=device, region=region_si, name="Gammap", value=0.1)

dg_physics.py

@@ -345,7 +345,7 @@ def mycb():
   set_node_values(device=device, region=region_ox, name="Potential", init_from="zero")
 
   set_parameter(name=GetContactBiasName("top"), value=0.0)
-  print d
+  print(d)
   set_parameter(device=device, region=region_si, name="Gamman", value=1.0)
   set_parameter(device=device, region=region_si, name="Gammanox", value=0.1)
   set_parameter(device=device, region=region_si, name="Gammap", value=0.1)