Skip to main content

Python interface to Piccolo and Picalor simulation kernel

Project description

Piccolo5/Piccolo6/Picalor6 and Ganessa_SIM API python calling interface

What is it?

A python package enabling users to call a set of Picwin32.dll or Ganessa_SIM.dll API functions and subroutines within python scripts. Initially provided for python 2.7 32 bits only. 1.9.x versions served as python 3.5 / 3.6 pre-release.

Starting with version 2.0.0, the wheel package is provided for 32 bits version for python 2.7 / 3.5 / 3.6 / 3.7 (Piccolo kernel is a 32 bits application) and 64 bits version for python 3.6 / 3.7 (Ganessa_SIM is available as 32 and 64 bits).

Starting with version 2.1.0, the wheel package is provided for 32 bits version for python 2.7 / 3.7 / 3.8 (Piccolo kernel is a 32 bits application) and 64 bits version for python 3.7 / 3.8 (Ganessa_SIM is available as 32 and 64 bits).

Installation

Windows only

pip install ganessa

Requirements (32 bits)

  1. python 2.7, 3.5, 3.6, 3.7, 3.8
  2. python requirements: numpy 1.11 or above (python 2.7) / numpy 1.13.1 or above (python 3.5 - 3.6) / numpy 1.15.1 (python 3.7 / 3.8)
  3. Piccolo or Picalor kernel library (picwin32.dll), starting from version 5 (141203) for python 2.7 / version 6 for python 3.x
  4. valid Piccolo or Picalor or Ganessa_SIM license

Requirements (64 bits)

  1. python 3.7 or 3.8
  2. numpy 1.13.1 or above (python 3.6) / numpy 1.15.1 (python 3.7 / 3.8)
  3. Ganessa_SIM_x64.dll (64 bit dll)
  4. valid Ganessa_SIM license

This tool expects Picwin32.dll to be in the PATH or in one of the following folders:

[ C: or D: ] / ['Program Files/Adelior/Piccolo5_'  or
               'Program Files (x86)/Gfi Progiciels/Piccolo6_'  or
               'Program Files (x86)/Gfi Progiciels/Picalor6_']
               'Program Files (x86)/Safege/Ganessa_']
             + ['fr' or 'uk' or 'esp' or 'eng'] + ['' or '_ck']

Or one of Ganessa_SIM.dll or Ganessa_TH.dll or Ganessa_SIM_x64.dll in:

[ C: or D: ] / ['Program Files (x86)/Safege/Ganessa_']
             + ['fr' or 'uk' or 'esp' or 'eng']

With custom installations, PICCOLO_DIR or GANESSA_DIR environment variables can be set to define the path.

Content

The package provides:
  1. ‘sim’ package:
    • a few basic functions for reading or loading a model, running hydraulic and WQ simulations.
    • ‘getter’ functions for individual objects and attributes, time series vectors, attribute vector of all object.
    • iterators over links, nodes, tanks, demands, and tables, or over Piccolo command language selections.
      • a Graph class built from the current model graph, with propagation functions.
  2. ‘th’ package: same functions for Picalor except running extended period simulations and time series getters.
  3. ‘util’ package: miscellaneous functions.
  4. ‘OpenFileMMI’ provides classes for opening dialog frame for a .dat/.bin model file, and output (result) file. It should be imported after ‘sim’ or ‘th’.
  5. ‘plot’ provides easy multiple simulation vs reference graphs plot settings - requires matplotlib.
  6. ‘sort’ provides a heapsort based on heapq.
  7. ‘midfile’ provides minimal mif/mid functions similar to shp/dbf shapefile handler (pyshp package).
  8. ‘geojsonfile’ provides minimal geojson functions similar to shp/dbf shapefile handler (pyshp package).
  9. ‘epanet’ provides epanet2 python API for win32 (thanks to Assela Pathirana - mailto: assela@pathirana.net) similar to EpanetTools-0.4.0; with epanet.getlinknodes corrected (from 0.9.2).
  10. ‘parallel’ provides a simple parallel simulations handling framework based on multiprocessing, returning either elementary or reduced results.

Model objects and parameters can be modified using Piccolo command language (see cmd, cmdfile and execute).

Recently introduced functions and classes may not be functional with old versions of Piccolo / Ganessa_SIM kernel.

Documentation

Documentation is provided as pyGanessa.html in the intallation folder.

History of the README document

Created 2013-07-04:
  • Revised 2015-05-03: since 2014-12-04 Picwin32.dll is compatible with this API.
  • Revised 2016-07-07: provided as .rst.
  • Revised 2017-08-08: install using pip; Piccolo/Ganessa dll folder search order.
  • Revised 2017-09-12: split sim into core, core_sim, core_th.
  • Revised 2017-11-13: added sort, midfile, epanet modules.
  • Revised 2017-11-30: added parallel.
  • Revised 2018-03-29: minor changes / extension to python 3.5-3.6.
  • Revised 2018-06-08: OpenFileMMI and epanet details, added example.
  • Revised 2018-08-16: 64 bits compatibility (Ganessa_SIM and python 3.7 / 3.8 only).
  • Revised 2018-08-27: numpy requirements.
  • Revised 2019-09-13: minor changes.
  • Revised 2019-10-09: extension to python 3.8.
  • Revised 2019-10-15: README rst fix.
  • Revised 2020-02-05: added sim.Graph; README rst fix.
  • Revised 2020-11-20: README rst fix.

Example

# -*- coding: utf-8 -*-
'''Test program for ganessa
requires a Piccolo model 'model.dat' in the current folder'''
from __future__ import (print_function, unicode_literals)
import os
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

import ganessa.sim as pic

# This code sample runs under python 2.7, 3.5, 3.6, 3.7
# but requires ganessa 2.1.1 or higher

if __name__ == '__main__':
    wkfold = os.path.dirname(__file__)
    if not wkfold:
        wkfold = os.getcwd()
        print('Current folder is:', wkfold)
    os.chdir(wkfold)
    # will create work files here
    pic.init(wkfold)
    # loads the model - change the name here...
    pic.cmdfile('model.dat')
    # select small pipe (diameter below 155)
    # using iterator and D getter, by id
    pipefilter = lambda d: (d > 0 and d < 155)
    pipes = [p for p in pic.Links() if pipefilter(pic.linkattr(p, 'D'))]
    # using iterator and D getter, by index (less concise but slightly faster)
    pipes = [pic.getid(pic.LINK, k) for k in range(1, pic.nbobjects(pic.LINK)+1)
                                    if pipefilter(pic.nlinkattr(k, 'D'))]
    # Most Piccoloic, using (french) command language selection
    pic.setlang('French')
    pipes = list(pic.Selected('TUYAU (D < 155) FIN', return_type=False))
    # minimal setter functions for links, nodes, tanks attributes
    pid = pipes[len(pipes)//2]
    diam = pic.linkattr(pid, 'D')
    pic.setlinkattr(pid, 'D', diam*2)
    # language-independant commands
    KWEND = pic.getkeyword(pic.DICT.END)
    pic.gencmdw(pic.M.MOD, pic.M.MOD.MULT, pic.M.NONE, 'D 2')
    pic.gencmd(pic.M.LNK, pic.M.LNK.NONE, pic.M.NONE, pid + ' ' + KWEND, 1)
    # Use command language to set / modify / delete objects
    pic.execute('MODIF', 'DETR ARC {} FIN'.format(pid))
    pic.execute('MODIF', 'FUSION NOEU (AA = 1) FIN',
                         'FUSION NOEU (AA = 1) FIN')
    # or change parameters
    pic.execute('SIMUL', 'CRITX 0.05', 'MXITER 75')
    # get demand profile TS
    nb = pic.tsdemandlen('DEFAUT', 'HAUT')
    print('Profile length for area HAUT:', nb)
    res = pic.tsdemand('DEFAUT', 'HAUT')
    # Unlike epanet, first run simulation then browse/ query results
    pic.full_solveH(silent=True, retry=True)
    print('Simulation done.')
    pic.browseH(3600*8.5)
    fig = plt.figure(figsize=(16-1, 9-1))
    # plot 3D elevation and hydraulic head
    ax = fig.gca(projection='3d')
    for a in pic.Links():
        x, y, z, v, nb = pic.linkXYZV(a, 'CH')
        if nb > 0:
            ax.plot(x, y, z, 'r-')
            ax.plot(x, y, v, 'b-')
    plt.title('Plan de charge a t= 08:30')
    plt.show()
    # print min/max pressure over the EPS
    attr = 'P'
    minpax_p = list(pic.getMinMax(pic.NODE, attr))
    for nid, vmin, vmax, vmoy in minpax_p[:10]:
        print(nid, attr+'min:', vmin, attr+'max:', vmax)
    print('\t...')
    for nid, vmin, vmax, vmoy in minpax_p[-10:]:
        print(nid, attr+'min:', vmin, attr+'max:', vmax)
    # get current volume simulation (internal variable steps) and measurement TS
    for tid in pic.Tanks():
        t1, v1, nb1 = pic.tsval(pic.TANK, tid, 'VC')
        t2, v2, nb2 = pic.msval(pic.TANK, tid, 'VC')
    # get current volume simulation results at measurement sampling
    for tid in pic.Tanks():
        t3, v3, nb3 = pic.tsvalbymts(pic.TANK, tid, 'VC')
    # export as Epanet .inp
    pic.exportEpanet('model.inp')
    pic.execute('XXMEM ARC : D')
    pic.savemodel('model_xx.dat', extra_data=[('ARC (D > 149) FIN', 'XX')])
    # WQ source data optional argument is a specific node - not used here
    print('\nWQ source data (input forcing and TS at nodes):')
    for n, cd, attr, t, c, nb in pic.WQSources():
        print(n, cd, attr, nb, t, c)
    print('Done.')
       # use reset() for switching model, close() to terminate
    pic.close()

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for ganessa, version 2.1.6
Filename, size File type Python version Upload date Hashes
Filename, size ganessa-2.1.6-cp27-cp27m-win32.whl (3.3 MB) File type Wheel Python version cp27 Upload date Hashes View
Filename, size ganessa-2.1.6-cp37-cp37m-win32.whl (3.2 MB) File type Wheel Python version cp37 Upload date Hashes View
Filename, size ganessa-2.1.6-cp37-cp37m-win_amd64.whl (817.7 kB) File type Wheel Python version cp37 Upload date Hashes View
Filename, size ganessa-2.1.6-cp38-cp38-win32.whl (3.2 MB) File type Wheel Python version cp38 Upload date Hashes View
Filename, size ganessa-2.1.6-cp38-cp38-win_amd64.whl (878.8 kB) File type Wheel Python version cp38 Upload date Hashes View
Filename, size ganessa-2.1.6-cp39-cp39-win32.whl (2.9 MB) File type Wheel Python version cp39 Upload date Hashes View
Filename, size ganessa-2.1.6-cp39-cp39-win_amd64.whl (896.4 kB) File type Wheel Python version cp39 Upload date Hashes View

Supported by

Pingdom Pingdom Monitoring Google Google Object Storage and Download Analytics Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN DigiCert DigiCert EV certificate StatusPage StatusPage Status page