Running BW1484 Benchmark on TITAN

In [13]:

    

import lxml.etree as ET
import string
xml = ET.parse('/data/work/tme_analysis/early_science/bw1484/bw1484_I.xml')
out = ET.tostring(xml, pretty_print=True, method="html")

root = xml.getroot()

Nr         = 0
mesh       = 0
N          = 0
num_groups = 0
for n in root:
        
    print (n.get('name'))
    
    if n.get('name') == "CORE":
        for m in n:
            d = dict(m.attrib)
            if d.has_key('value'):
                print "%15s %12s" % (d['name'], d['value'])
                
    if n.get('name') == "ASSEMBLIES":
        for m in n.iter():
            d = dict(m.attrib)
            if d['name'] == "num_pins":
                print "%15s %12s" % (d['name'], d['value'])
                Nr = string.atoi(d['value'])
            if d['name'] == "mesh":
                print "%15s %12s" % (d['name'], d['value'])
                
    if n.get('name') == "INSILICO":
        for m in n.iter():
            d = dict(m.attrib)
            if d['name'] == "mesh":
                print "%15s %12s" % (d['name'], d['value'])
                mesh = string.atoi(d['value'])
            if d['name'] == "Pn_order":
                print "%15s %12s" % (d['name'], d['value'])
                N = string.atoi(d['value'])
            if d['name'] == "num_groups":
                print "%15s %12s" % (d['name'], d['value'])
                num_groups = string.atoi(d['value'])


    

    




case_id
STATE
CORE
         apitch       98.148
       assm_map       {FUEL}
         bc_rad       vacuum
      core_size            1
         height       145.00
      lower_mat           al
    lower_thick         5.08
    lower_vfrac          1.0
     rated_flow   773.626492
    rated_power        159.1
          shape          {1}
ASSEMBLIES
           mesh            2
       num_pins           60
INSILICO
       Pn_order            3
           mesh           20
     num_groups           23
EDITS

In [14]:

    

print(out)


    

    




<?xml-stylesheet version="1.0" type="text/xsl" href="PL9.xsl"?>
<ParameterList name="CASEID">
  <Parameter name="case_id" type="string" value="BW 1484 Core I"></Parameter>
  <ParameterList name="STATE">
    <Parameter name="boron" type="double" value="0"></Parameter>
    <Parameter name="feedback" type="string" value="off"></Parameter>
    <Parameter name="modden" type="double" value="0.997992"></Parameter>
    <Parameter name="power" type="double" value="0.0"></Parameter>
    <Parameter name="sym" type="string" value="qtr"></Parameter>
    <Parameter name="tfuel" type="double" value="294.15"></Parameter>
    <Parameter name="tinlet" type="double" value="21"></Parameter>
  </ParameterList>
  <ParameterList name="CORE">
    <ParameterList name="Materials">
      <ParameterList name="Material_al">
        <Parameter name="density" type="double" value="2.688"></Parameter>
        <Parameter name="key_name" type="string" value="al"></Parameter>
        <Parameter name="mat_fracs" type="Array(double)" value="{1}"></Parameter>
        <Parameter name="mat_names" type="Array(string)" value="{al6061}"></Parameter>
      </ParameterList>
      <ParameterList name="Material_b4c">
        <Parameter name="density" type="double" value="0.879"></Parameter>
        <Parameter name="key_name" type="string" value="b4c"></Parameter>
        <Parameter name="mat_fracs" type="Array(double)" value="{1}"></Parameter>
        <Parameter name="mat_names" type="Array(string)" value="{b4c}"></Parameter>
      </ParameterList>
      <ParameterList name="Material_he">
        <Parameter name="density" type="double" value="0.0001786"></Parameter>
        <Parameter name="key_name" type="string" value="he"></Parameter>
        <Parameter name="mat_fracs" type="Array(double)" value="{1}"></Parameter>
        <Parameter name="mat_names" type="Array(string)" value="{he}"></Parameter>
      </ParameterList>
    </ParameterList>
    <Parameter name="apitch" type="double" value="98.148"></Parameter>
    <Parameter name="assm_map" type="Array(string)" value="{FUEL}"></Parameter>
    <Parameter name="bc_rad" type="string" value="vacuum"></Parameter>
    <Parameter name="core_size" type="int" value="1"></Parameter>
    <Parameter name="height" type="double" value="145.00"></Parameter>
    <Parameter name="lower_mat" type="string" value="al"></Parameter>
    <Parameter name="lower_thick" type="double" value="5.08"></Parameter>
    <Parameter name="lower_vfrac" type="double" value="1.0"></Parameter>
    <Parameter name="rated_flow" type="double" value="773.626492"></Parameter>
    <Parameter name="rated_power" type="double" value="159.1"></Parameter>
    <Parameter name="shape" type="Array(int)" value="{1}"></Parameter>
  </ParameterList>
  <ParameterList name="ASSEMBLIES">
    <ParameterList name="Assembly_FUEL">
      <ParameterList name="CellMaps">
        <ParameterList name="CellMap_FUEL">
          <Parameter name="cell_map" type="Array(string)" value="{3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,1,1,1,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3}"></Parameter>
          <Parameter name="label" type="string" value="FUEL"></Parameter>
        </ParameterList>
        <ParameterList name="CellMap_PLUG">
          <Parameter name="cell_map" type="Array(string)" value="{3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3}"></Parameter>
          <Parameter name="label" type="string" value="PLUG"></Parameter>
        </ParameterList>
      </ParameterList>
      <ParameterList name="Cells">
        <ParameterList name="Cell_1">
          <Parameter name="label" type="string" value="1"></Parameter>
          <Parameter name="mats" type="Array(string)" value="{U246,he,al}"></Parameter>
          <Parameter name="num_rings" type="int" value="3"></Parameter>
          <Parameter name="radii" type="Array(double)" value="{0.515,0.522,0.603}"></Parameter>
          <Parameter name="type" type="string" value="fuel"></Parameter>
        </ParameterList>
        <ParameterList name="Cell_2">
          <Parameter name="label" type="string" value="2"></Parameter>
          <Parameter name="mats" type="Array(string)" value="{al}"></Parameter>
          <Parameter name="num_rings" type="int" value="1"></Parameter>
          <Parameter name="radii" type="Array(double)" value="{0.603}"></Parameter>
          <Parameter name="type" type="string" value="other"></Parameter>
        </ParameterList>
        <ParameterList name="Cell_3">
          <Parameter name="label" type="string" value="3"></Parameter>
          <Parameter name="mats" type="Array(string)" value="{mod}"></Parameter>
          <Parameter name="num_rings" type="int" value="1"></Parameter>
          <Parameter name="radii" type="Array(double)" value="{0.8}"></Parameter>
          <Parameter name="mesh" type="int" value="2"></Parameter>
          <Parameter name="type" type="string" value="other"></Parameter>
        </ParameterList>
      </ParameterList>
      <ParameterList name="Fuels">
        <ParameterList name="Fuel_U246">
          <Parameter name="density" type="double" value="10.220"></Parameter>
          <Parameter name="enrichments" type="Array(double)" value="{2.459}"></Parameter>
          <Parameter name="fuel_names" type="Array(string)" value="{u-235}"></Parameter>
          <Parameter name="key_name" type="string" value="U246"></Parameter>
          <Parameter name="thden" type="double" value="94.5"></Parameter>
        </ParameterList>
      </ParameterList>
      <ParameterList name="SpacerGrids">
        <ParameterList name="Grid_END">
          <Parameter name="height" type="double" value="2.54"></Parameter>
          <Parameter name="label" type="string" value="END"></Parameter>
          <Parameter name="mass" type="double" value="1638"></Parameter>
          <Parameter name="material" type="string" value="al"></Parameter>
        </ParameterList>
      </ParameterList>
      <Parameter name="axial_elevations" type="Array(double)" value="{0.0,0.3175,145.0}"></Parameter>
      <Parameter name="axial_labels" type="Array(string)" value="{PLUG,FUEL}"></Parameter>
      <Parameter name="grid_elev" type="Array(double)" value="{1.5875}"></Parameter>
      <Parameter name="grid_map" type="Array(string)" value="{END}"></Parameter>
      <Parameter name="label" type="string" value="FUEL"></Parameter>
      <Parameter name="num_pins" type="int" value="60"></Parameter>
      <Parameter name="ppitch" type="double" value="1.6358"></Parameter>
      <Parameter name="title" type="string" value="BW 1484 Core"></Parameter>
    </ParameterList>
  </ParameterList>
  <ParameterList name="INSILICO">
    <Parameter name="Pn_correction" type="bool" value="false"></Parameter>
    <Parameter name="Pn_order" type="int" value="3"></Parameter>
    <Parameter name="cell_homogenize" type="bool" value="false"></Parameter>
    <Parameter name="dimension" type="int" value="3"></Parameter>
    <Parameter name="do_transport" type="bool" value="false"></Parameter>
    <ParameterList name="eigenvalue_db">
      <Parameter name="L2_tolerance" type="double" value="1.0e-3"></Parameter>
    </ParameterList>
    <Parameter name="eq_set" type="string" value="ld"></Parameter>
    <Parameter name="mat_library" type="string" value="bw1484_comp.sh5"></Parameter>
    <Parameter name="max_delta_z" type="double" value="5.0"></Parameter>
    <Parameter name="mesh" type="int" value="20"></Parameter>
    <Parameter name="new_grp_bounds" type="Array(double)" value="{8.2085e+05,1.1109e+05,5.5308e+03,1.8644e+02,3.7612e+01,3.5379e+01,2.7697e+01,2.1684e+01,2.0397e+01,1.5968e+01,7.1500e+00,6.7000e+00,6.3000e+00,1.0970e+00,1.0450e+00,9.5000e-01,3.5000e-01,2.0600e-01,1.0700e-01,5.8000e-02,2.5000e-02,1.0000e-02,1.0000e-05}"></Parameter>
    <Parameter name="num_blocks_i" type="int" value="1"></Parameter>
    <Parameter name="num_blocks_j" type="int" value="1"></Parameter>
    <Parameter name="num_groups" type="int" value="23"></Parameter>
    <Parameter name="num_sets" type="int" value="1"></Parameter>
    <Parameter name="num_z_blocks" type="int" value="1"></Parameter>
    <Parameter name="problem_name" type="string" value="out_bw1484_I_sn"></Parameter>
    <ParameterList name="quadrature_db">
      <Parameter name="azimuthals_octant" type="int" value="12"></Parameter>
      <Parameter name="polars_octant" type="int" value="6"></Parameter>
      <Parameter name="quad_type" type="string" value="qr"></Parameter>
    </ParameterList>
    <Parameter name="tolerance" type="double" value="1e-3"></Parameter>
    <Parameter name="transport" type="string" value="deterministic"></Parameter>
    <ParameterList name="upscatter_db">
      <Parameter name="tolerance" type="double" value="1.0e-8"></Parameter>
    </ParameterList>
    <Parameter name="xs_library" type="string" value="lib252_hetbondoneabs-noabssigp"></Parameter>
  </ParameterList>
  <ParameterList name="EDITS">
    <Parameter name="axial_edit_bounds" type="Array(double)" value="{0.3175,11.3175,22.3175,33.3175,44.3175,55.3175,66.3175,77.3175,88.3175,99.3175,110.3175,121.3175,132.3175,145.0}"></Parameter>
  </ParameterList>
</ParameterList>

In [15]:

    

# number of axial cells with 1 cm spacing
Nz = 154

# total mesh axially (20 cells per pin in core, 2 cells per pin in water)
Nr = 552


    

In [7]:

    

bw = Image(filename='/data/work/tme_analysis/early_science/bw1484/out.png')


    

In [8]:

    

display(bw)


    

In [4]:

    

%loadpy /data/denovo/Exnihilo/packages/Denovo/tools/kbaprofiler.py


    

In [19]:

    

###############################################################################
## kbaprofiler.py
## Thomas M. Evans
## Wed Mar 10 09:01:32 2010
## $Id: template.py,v 1.3 2008/01/02 17:18:47 9te Exp $
###############################################################################
## Copyright (C) 2008 Oak Ridge National Laboratory, UT-Battelle, LLC.
###############################################################################

import sys, os

##---------------------------------------------------------------------------##
## PROBLEM VARIABLES : SETUP BY USER
##---------------------------------------------------------------------------##

# Number of cells in each dimension
Ni = Nr
Nj = Nr
Nk = Nz

# Pn order
Pn = N

# Number of groups in memory (determined by loading the cross section set)
Ng = num_groups

# Specify a range of solved groups by setting range to 1 and then specifying
# the first and last group in range [g_first, g_last] (inclusive), the groups
# are defined over the range [0, Ng)
Range   = 0    # set to 1 if a solved range is specified
g_first = 0    # first group in range [0, Ng)
g_last  = 26   # last  group in range [0, Ng)

# Quadrature Setup (choose Level-symmetric or Product)
Levelsym = 0
Product  = 1

# if level-symmetric, set the Sn order
Sn = 8

# if product set the number of polar per octant and the number of azimuthals
# per polar
Npolar  = 6
Nazimth = 12

# set the equation set by placing setting it equal to 1
sc  = 0  # step-characteristics
tld = 0  # tri-linear discontinuous
ld  = 1  # linear-discontinuous
twd = 0  # theta-weighted-diamond
wdd = 0  # weighted-diamond-difference (use for regular weighted diamond and
         # diamond difference with flux-fixup)

# Domain Filters
Min_cores        = 2000   # minimum cores
Max_cores        = 25000  # maximum cores
Core_multiple    = 16     # restricted to multiples of this core count
Min_Cells_block  = 2000   # minimum cells per block
Min_Cells_domain = 1      # minimum cells per domain
Max_MB_domain    = 1000.0 # maximum state size per domain

# Ignore uniform block sizes (domain/block/efficiencies will be approximate)
Ignore_uniformity = 0 # set to 1 to ignore uniformity

# Minimum desired parallel efficiency
Min_eff = 0.80

##---------------------------------------------------------------------------##
## NO NEED FOR EDITING BELOW HERE
##---------------------------------------------------------------------------##
## Computer memory conversions

# Byte-sizes of words
DOUBLE = 8
INT    = 4

# Conversions to memory units
KB = 2.0**10
MB = 2.0**20
GB = 2.0**30
TB = 2.0**40

##---------------------------------------------------------------------------##
## Problem sizes

# angles per octant
M = 0
if Levelsym == 1:
    M = Sn * (Sn + 2) / 8
elif Product == 1:
    M = Npolar * Nazimth

# total angles
Na = M * 8

# total cells
Nc = Ni * Nj * Nk

# unknowns per cell
Nu = 1
if tld == 1: Nu = 8
if ld  == 1: Nu = 4

# number of moments
Nm = (Pn + 1) * (Pn + 1)

# number of solved groups
Ng_s = Ng
if Range == 1:
    Ng_s = 1 + g_last - g_first

# Global size of problem state
N_unknowns = Nu * Nc * Ng * Nm
B_unknowns = N_unknowns * DOUBLE / GB

# Number of unknowns per spatial cell
Nc_unknowns = Nu * Ng * Nm
Bc_unknowns = Nc_unknowns * DOUBLE / KB

# Global size of solved unknowns (degrees-of-freedom)
N_solved = Nu * Nc * Ng_s * Nm * Na

##---------------------------------------------------------------------------##
## Print Data Tables

print "********************************************************************"
print "Denovo KBA Problem Profiling"
print "===================================================================="
print "Problem Setup"
print "--------------------------------------------------------------------"
print
print "Number of cells in I           = %6i" % (Ni)
print "Number of cells in J           = %6i" % (Nj)
print "Number of cells in K           = %6i" % (Nk)
print "Number of spatial moments/cell = %6i" % (Nu)
print "Number of angular moments      = %6i" % (Nm)
print "Number of angles               = %6i" % (Na)
print "Number of angles per octant    = %6i" % (M)
print "Number of loaded energy groups = %6i" % (Ng)
print "Number of solved energy groups = %6i" % (Ng_s)
print
print "--------------------------------------------------------------------"
print "Global Sizes"
print "--------------------------------------------------------------------"
print
print "Number of total cells   = %16i" % (Nc)
print "Size of problem state   = %16i (words)" % (N_unknowns)
print "Size of problem state   = %16.5e (GB)" % (B_unknowns)
print
print "Number of unknowns/cell = %16i (words)" % (Nc_unknowns)
print "Memory per cell         = %16.5e (KB)" % (Bc_unknowns)
print
print "Solved DoF              = %16.5e (total unknowns)" % (N_solved)
print
print "--------------------------------------------------------------------"
print "Parallel Domain Analysis"
print "--------------------------------------------------------------------"
print

# Pi = number of partitions in I
# Pj = number of partitions in J

print "%7s %4s %4s %4s %5s %6s %8s %10s %12s" % \
    ("Domains", "Pi", "Pj", "Bk", "Eff", "Cells", "Cells", "MB/Domain", "Space-Angle")
print "%35s %8s %23s" % ("/Block", "/Domain", "/Domain")
print "...................................................................."

# maximum attainable efficiency
max_eff = 0.0

# loop over potential partitions in i, j
for Pj in xrange(1, Nj):
    for Pi in xrange(1, Ni):

        # check for equal-sized blocks
        if Ignore_uniformity != 1:
            if Ni % Pi != 0: continue
            if Nj % Pj != 0: continue

        # check if Pi * Pj > Min_cores
        N_domains = Pi * Pj
        if (N_domains % Core_multiple != 0): continue
        if (N_domains < Min_cores):          continue
        if (N_domains > Max_cores):          continue

        # number of cells in (I,J) in 1 partition
        Ib = Ni / float(Pi)
        Jb = Nj / float(Pj)

        # loop over z-blocks
        for Bk in xrange(1, Nk):

            # check for equal-sized blocks
            if Nk % Bk != 0: continue

            # number of cells in K-block
            Kb = Nk / Bk

            # cells per block
            Nb = Ib * Jb * Kb

            if Nb < Min_Cells_block: continue

            # cells per domain
            Dc = Ib * Jb * Nk

            if Dc < Min_Cells_domain: continue

            # memory per domain in MB
            Dm = Dc * Bc_unknowns * KB / MB

            if Dm > Max_MB_domain : continue

            # number of space-angle unknowns per block
            Usa = Nu * Nm * Na * Nb

            # calculate hypothetical efficiency
            eff = float(2 * M * Bk) / float(2 * M * Bk + Pi + Pj - 2)

            # calculate maximum efficiency
            max_eff = max(eff, max_eff)

            # print out
            if eff > Min_eff:
                print "%7i %4i %4i %4i %5.2f %6i %8i %10.3e %12i" % \
                    (N_domains, Pi, Pj, Bk, eff, Nb, Dc, Dm, Usa)

print "...................................................................."
print
print "Maximum efficiency attainable : %10.6f" % (max_eff)
print
print "********************************************************************"

###############################################################################
##                            end of kbaprofiler.py
###############################################################################


    

    




********************************************************************
Denovo KBA Problem Profiling
====================================================================
Problem Setup
--------------------------------------------------------------------

Number of cells in I           =    552
Number of cells in J           =    552
Number of cells in K           =    154
Number of spatial moments/cell =      4
Number of angular moments      =     16
Number of angles               =    576
Number of angles per octant    =     72
Number of loaded energy groups =     23
Number of solved energy groups =     23

--------------------------------------------------------------------
Global Sizes
--------------------------------------------------------------------

Number of total cells   =         46924416
Size of problem state   =      69072740352 (words)
Size of problem state   =      5.14632e+02 (GB)

Number of unknowns/cell =             1472 (words)
Memory per cell         =      1.15000e+01 (KB)

Solved DoF              =      3.97859e+13 (total unknowns)

--------------------------------------------------------------------
Parallel Domain Analysis
--------------------------------------------------------------------

Domains   Pi   Pj   Bk   Eff  Cells    Cells  MB/Domain  Space-Angle
                             /Block  /Domain                 /Domain
....................................................................
   2208  184   12    7  0.84   3036    21252  2.387e+02    111919104
   2208   92   24    7  0.90   3036    21252  2.387e+02    111919104
   3312  138   24    7  0.86   2024    14168  1.591e+02     74612736
   2208   24   92    7  0.90   3036    21252  2.387e+02    111919104
   3312   24  138    7  0.86   2024    14168  1.591e+02     74612736
   2208   12  184    7  0.84   3036    21252  2.387e+02    111919104
....................................................................

Maximum efficiency attainable :   0.898396

********************************************************************

In [14]:

    

Nb = 2208
Px = 92
Py = 24
Bz = 7
Ng = 252


    

In [15]:

    

num_sets = Ng / 6
num_sets


    

    
Out[15]:





42

In [16]:

    

domains = num_sets * Nb


    

In [17]:

    

domains


    

    
Out[17]:





92736

In [18]:

    

num_nodes = domains / 16
print(num_nodes)


    

    




5796

In [19]:

    

domains % 16


    

    
Out[19]:





0

In [20]:

    

Nz % Bz


    

    
Out[20]:





0

In [40]:

    

xml = ET.parse('/data/work/tme_analysis/early_science/bw1484/bw1484_I252G.xml')
out = ET.tostring(xml, pretty_print=True, method="html")

root = xml.getroot()

for n in root[4]:
    print n.items()


    

    




[('name', 'Pn_correction'), ('type', 'bool'), ('value', 'false')]
[('name', 'Pn_order'), ('type', 'int'), ('value', '3')]
[('name', 'cell_homogenize'), ('type', 'bool'), ('value', 'false')]
[('name', 'dimension'), ('type', 'int'), ('value', '3')]
[('name', 'do_transport'), ('type', 'bool'), ('value', 'true')]
[('name', 'eigenvalue_db')]
[('name', 'eq_set'), ('type', 'string'), ('value', 'ld')]
[('name', 'mat_library'), ('type', 'string'), ('value', 'bw1484_comp.sh5')]
[('name', 'max_delta_z'), ('type', 'double'), ('value', '1.0')]
[('name', 'mesh'), ('type', 'int'), ('value', '20')]
[('name', 'num_blocks_i'), ('type', 'int'), ('value', '92')]
[('name', 'num_blocks_j'), ('type', 'int'), ('value', '24')]
[('name', 'num_groups'), ('type', 'int'), ('value', '252')]
[('name', 'num_sets'), ('type', 'int'), ('value', '42')]
[('name', 'num_z_blocks'), ('type', 'int'), ('value', '7')]
[('name', 'problem_name'), ('type', 'string'), ('value', 'out_bw1484_I252G_sn')]
[('name', 'quadrature_db')]
[('name', 'tolerance'), ('type', 'double'), ('value', '1e-3')]
[('name', 'transport'), ('type', 'string'), ('value', 'deterministic')]
[('name', 'upscatter_db')]
[('name', 'xs_library'), ('type', 'string'), ('value', 'lib252_hetbondoneabs-noabssigp')]

In [23]:

    

num_nodes = domains / 8
print(num_nodes)


    

    




11592

In [25]:

    

domains


    

    
Out[25]:





92736

In [26]:

    

num_nodes * 8


    

    
Out[26]:





92736