# $Id: 3298c1e on Wed Mar 27 15:18:09 2013 +0100 by Gianluca De Santis $
##
##  This file is part of pyFormex 0.9.1  (Tue Oct 15 21:05:25 CEST 2013)
##  pyFormex is a tool for generating, manipulating and transforming 3D
##  geometrical models by sequences of mathematical operations.
##  Home page: http://pyformex.org
##  Project page:  http://savannah.nongnu.org/projects/pyformex/
##  Copyright 2004-2013 (C) Benedict Verhegghe (benedict.verhegghe@ugent.be)
##  Distributed under the GNU General Public License version 3 or later.
##
##  This program is free software: you can redistribute it and/or modify
##  it under the terms of the GNU General Public License as published by
##  the Free Software Foundation, either version 3 of the License, or
##  (at your option) any later version.
##
##  This program is distributed in the hope that it will be useful,
##  but WITHOUT ANY WARRANTY; without even the implied warranty of
##  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
##  GNU General Public License for more details.
##
##  You should have received a copy of the GNU General Public License
##  along with this program.  If not, see http://www.gnu.org/licenses/.
##
"""General framework for attributing properties to geometrical elements.

Properties can really be just about any Python object.
Properties can be attributed to a set of geometrical elements.
"""
from __future__ import print_function

import pyformex as pf
from flatkeydb import FlatDB
from mydict import Dict,CDict
from arraytools import *

#################################################################
# This first part still needs to be changed.
# It should probably be moved to a separate module

_matDB = None
_secDB = None

def setMaterialDB(mat):
    """Set the global materials database.

    If mat is a MaterialDB, it will be used as the global MaterialDB.
    Else, a new global MaterialDB will be created, initialized from
    the argument mat.
    """
    global _matDB
    if isinstance(mat,MaterialDB) or mat is None:
        _matDB = mat
    else:
        _matDB = MaterialDB(mat)

def setSectionDB(sec):
    """Set the global sections database.

    If sec is a SectionDB, it will be used as the global SectionDB.
    Else, a new global SectionDB will be created, initialized from
    the argument sec.
    """
    global _secDB
    if isinstance(sec,SectionDB) or sec is None:
        _secDB = sec
    else:
        _secDB = SectionDB(sec)


class Database(Dict):
    """A class for storing properties in a database."""

    def __init__(self,data={}):
        """Initialize a database.

        The database can be initialized with a dict.
        """
        Dict.__init__(self,data)


    def readDatabase(self,filename,*args,**kargs):
        """Import all records from a database file.

        For now, it can only read databases using flatkeydb.
        args and kargs can be used to specify arguments for the
        FlatDB constructor.
        """
        mat = FlatDB(*args,**kargs)
        mat.readFile(filename)
        for k,v in mat.iteritems():
            self[k] = Dict(v)


class MaterialDB(Database):
    """A class for storing material properties."""

    def __init__(self,data={}):
        """Initialize a materials database.

        If data is a dict, it contains the database.
        If data is a string, it specifies a filename where the
        database can be read. If it is an empty string, the
        configured database name will be used.
        """
        Database.__init__(self,{})
        if data == '':
            data = pf.cfg.get('prop/matdb',{})
        if type(data) == str:
            self.readDatabase(data,['name'],beginrec='material',endrec='endmaterial')
        elif type(data) == dict:
            self.update(data)
        else:
            raise ValueError,"Expected a filename or a dict."


class SectionDB(Database):
    """A class for storing section properties."""

    def __init__(self,data={}):
        """Initialize a section database.

        If data is a dict, it contains the database.
        If data is a string, it specifies a filename where the
        database can be read. If it is an empty string, the
        configured database name will be used.
        """
        Database.__init__(self,{})
        if data == '':
            data = pf.cfg.get('prop/secdb',{})
        if type(data) == str:
            self.readDatabase(data,['name'],beginrec='section',endrec='endsection')
        elif type(data) == dict:
            self.update(data)
        else:
            raise ValueError,"Expected a filename or a dict."


class ElemSection(CDict):
    """Properties related to the section of an element.

    An element section property can hold the following sub-properties:

    section
      the geometric properties of the section. This can be a dict
      or a string. If it is a string, its value is looked up in the global
      section database. The section dict should at least have a key
      'sectiontype', defining the type of section.

      Currently the following sectiontype values are known by module
      :mod:`fe_abq` for export to Abaqus/Calculix:

        - 'solid'   : a solid 2D or 3D section,
        - 'circ'    : a plain circular section,
        - 'rect'    : a plain rectangular section,
        - 'pipe'    : a hollow circular section,
        - 'box'     : a hollow rectangular section,
        - 'I'       : an I-beam,
        - 'general' : anything else (automatically set if not specified).
        - 'rigid'   : a rigid body

      .. note: Currently only 'solid', 'general' and 'rigid' are allowed.

      The other possible (useful) keys in the section dict depend on the
      sectiontype. Again for :mod:`fe_abq`:

        - for sectiontype 'solid' : thickness
        - the sectiontype 'general': cross_section, moment_inertia_11,
          moment_inertia_12, moment_inertia_22, torsional_constant
        - for sectiontype 'circ': radius
        - for sectiontype 'rigid': refnode, density, thickness

    material
      the element material. This can be a dict or a string.
      Currently known keys to fe_abq.py are: young_modulus,
      shear_modulus, density, poisson_ratio . (see fmtMaterial in fe_abq)
      It should not be specified for rigid sections.

    orientation
      - a Dict, or
      - a list of 3 direction cosines of the first beam section axis.

    """
    def __init__(self,section=None,material=None,orientation=None,**kargs):
        """Create a new element section property. Empty by default."""
        if _matDB is None:
            setMaterialDB({})
        if _secDB is None:
            setSectionDB({})
        CDict.__init__(self,kargs)
        self.addMaterial(material)
        self.addSection(section)
        if orientation is not None:
            self.orientation = orientation



    def addSection(self, section):
        """Create or replace the section properties of the element.

        If 'section' is a dict, it will be added to the global SectionDB.
        If 'section' is a string, this string will be used as a key to
        search in the global SectionDB.
        """
        if isinstance(section, str):
            if section in _secDB:
                self.section = _secDB[section]
            else:
                pf.warning("Section '%s' is not in the database" % section)
        elif isinstance(section,dict):
            # WE COULD ADD AUTOMATIC CALCULATION OF SECTION PROPERTIES
            #self.computeSection(section)
            #print(section)
            _secDB[section['name']] = CDict(section)
            self.section = _secDB[section['name']]
        elif section==None:
            self.section = section
        else:
            raise ValueError,"Expected a string or a dict"


    def computeSection(self,section):
        """Compute the section characteristics of specific sections."""
        if 'sectiontype' not in section:
            return
        if section['sectiontype'] == 'circ':
            r = section['radius']
            A = pi * r**2
            I = pi * r**4 / 4
            section.update({'cross_section':A,
                            'moment_inertia_11':I,
                            'moment_inertia_22':I,
                            'moment_inertia_12':0.0,
                            'torsional_constant':2*I,
                            })
        else:
            raise ValueError,"Invalid sectiontype"


    def addMaterial(self, material):
        """Create or replace the material properties of the element.

        If the argument is a dict, it will be added to the global MaterialDB.
        If the argument is a string, this string will be used as a key to
        search in the global MaterialDB.
        """
        if isinstance(material, str) :
            if material in _matDB:
                self.material = _matDB[material]
            else:
                pf.warning("Material '%s'  is not in the database" % material)
        elif isinstance(material, dict):
            _matDB[material['name']] = CDict(material)
            self.material = _matDB[material['name']]
        elif material==None:
            self.material=material
        else:
            raise ValueError,"Expected a string or a dict"


class ElemLoad(CDict):
    """Distributed loading on an element."""

    def __init__(self,label=None,value=None,dir=None):
        """Create a new element load. Empty by default.

        An element load can hold the following sub-properties:

        - label: the distributed load type label
        - value: the magnitude of the distibuted load
        - dir: vector specifying the direction of the load
        """
        if label == 'GRAV':
            if dir is None:
                dir = [0, 0 ,-1]
            if value is None:
                value = 9.81
        Dict.__init__(self,{'label':label,'value':value,'dir':dir})


class EdgeLoad(CDict):
    """Distributed loading on an element edge."""

    def __init__(self,edge=-1,label=None,value=None):
        """Create a new element edge load. Empty by default.

        An element edgeload can hold the following sub-properties:
        - edge: the element edge number
        - label: the distributed load type label ('x','y','z').
        - value: the magnitude of the distibuted load.
        """
        Dict.__init__(self,{'edge':edge,'label':label,'value':value})


############## Basic property data classes ########################

class CoordSystem(object):
    """A class for storing coordinate systems."""

    valid_csys = 'RSC'

    def __init__(self,csys,cdata):
        """Create a new coordinate system.

        csys is one of 'Rectangular', 'Spherical', 'Cylindrical'. Case is
          ignored and the first letter suffices.
        cdata is a list of 6 coordinates specifying the two points that
          determine the coordinate transformation
        """
        try:
            csys = csys[0].upper()
            if not csys in CoordSystem.valid_csys:
                raise
            cdata = asarray(cdata).reshape(2,3)
        except:
            raise ValueError,"Invalid initialization data for CoordSystem"
        self.sys = csys
        self.data = cdata


class Amplitude(object):
    """A class for storing an amplitude.

    The amplitude is a list of tuples (time,value).

    `atime` (amplitude time) can be either STEP TIME (default in Abaqus)
    or TOTAL TIME

    `smoothing` (optional) is a float (from 0. to 0.5, suggested value 0.05)
    representing the fraction of the time interval before and after each time
    point during which the piecewise linear time variation will be replaced by
    a smooth quadratic time variation (avoiding infinite accelerations).
    Smoothing should be used in combination with TABULAR (set 0.05 as default
    value?)
    """

    def __init__(self,data,definition='TABULAR',atime='STEP TIME', smoothing=None):
        """Create a new amplitude."""
        if definition in [ 'TABULAR', 'SMOOTH STEP' ]:
            if atime in [ 'STEP TIME', 'TOTAL TIME' ]:
                self.data = checkArray(data,(-1,2),'f','i')
                self.type = definition
                self.atime = atime
                if definition == 'TABULAR':
                    if smoothing is not None :
                        self.type += ', SMOOTHING=%s'%smoothing
        else:
            raise ValueError,"Expected definition = 'TABULAR' or 'SMOOTH STEP'"


###################################################
############ Utility routines #####################


def checkIdValue(values):
    """Check that a variable is a list of (id,value) tuples

    `id` should be convertible to an int, value to a float.
    If ok, return the values as a list of (int,float) tuples.
    """
    try:
        l = [ len(v) for v in values ]
        if min(l) == 2 and max(l) == 2:
            return [ (int(i),float(v)) for i,v in values ]
    except:
        raise ValueError,"Expected a list of (int,float) tuples"


def checkArrayOrIdValue(values):
    """Check that a variable is an list of values or (id,value) tuples

    This convenience function checks that the argument is either:

    - a list of 6 float values (or convertible to it), or
    - a list of (id,value) tuples where id is convertible to an int,
      value to a float.

    If ok, return the values as a list of (int,float) tuples.
    """
    ##print("VALUES IN: %s" % values)
    try:
        v = checkArray(values,(6,),'f','i')
        w = where(v != 0.0)[0]
        values = [ (i,v[i]) for i in w ]
    except:
        values = checkIdValue(values)
    ##print("VALUES OUT: %s" % values)
    return values


def checkString(a,valid):
    """Check that a string a has one of the valid values.

    This is case insensitive, and returns the upper case string if valid.
    Else, an error is raised.
    """
    try:
        a = a.upper()
        if a in valid:
            return a
    except:
        print("Expected one of %s, got: %s" % (valid,a))
    raise ValueError


# Create automatic names for node and element sets

def autoName(base,*args):
    return (base + '_%s' * len(args)) % args

# The following are not used by the PropertyDB class,
# but may be convenient for the user in applications

def Nset(*args):
    return autoName('Nset',*args)

def Eset(*args):
    return autoName('Eset',*args)

#############################################################
##################### Properties Database ###################


def FindListItem(l,p):
    """Find the item p in the list l.

    If p is an item in the list (not a copy of it!), this returns
    its position. Else, -1 is returned.

    Matches are found with a 'is' function, not an '=='.
    Only the first match will be reported.
    """
    for i,j in enumerate(l):
        if j is p:
            return i
    return -1


def RemoveListItem(l,p):
    """Remove the item p from the list l.

    If p is an item in the list (not a copy of it!), it is removed from
    the list.
    Matches are found with a 'is' comparison. This is different from the
    normal Python list.remove() method, which uses '=='.
    As a result, we can find complex objects which do not allow '==',
    such as ndarrays.
    """
    i = FindListItem(l,p)
    if i >= 0:
        del l[i]


class PropertyDB(Dict):
    """A database class for all properties.

    This class collects all properties that can be set on a
    geometrical model.

    This should allow for storing:

    - materials
    - sections
    - any properties
    - node properties
    - elem properties
    - model properties (current unused: use unnamed properties)

    Materials and sections use their own database for storing. They can be
    specified on creating the property database. If not specified, default
    ones are created from the files distributed with pyFormex.
    """

    bound_strings = [ 'XSYMM', 'YSYMM', 'ZSYMM', 'ENCASTRE', 'PINNED' ]

    def __init__(self,mat='',sec=''):
        """Create a new properties database."""
        setMaterialDB(mat)
        setSectionDB(sec)
        Dict.__init__(self)
        self.prop = []
        self.nprop = []
        self.eprop = []
        #self.mprop = []

    @staticmethod
    def matDB():
        return _matDB

    @staticmethod
    def secDB():
        return _secDB

    @classmethod
    def autoName(clas,kind,*args):
        return autoName((kind+'set').capitalize(),*args)

    def setMaterialDB(self,aDict):
        """Set the materials database to an external source"""
        setMaterialDB(aDict)

    def setSectionDB(self,aDict):
        """Set the sections database to an external source"""
        setSectionDB(aDict)

    def print(self):
        """Print the property database"""
        print("General properties")
        for p in self.getProp(''):
            print(p)
        print("Node properties")
        for p in self.getProp('n'):
            print(p)
        print("Element properties")
        for p in self.getProp('e'):
            print(p)


    def Prop(self,kind='',tag=None,set=None,name=None,**kargs):
        """Create a new property, empty by default.

        A property can hold almost anything, just like any Dict type.
        It has however four predefined keys that should not be used for
        anything else than explained hereafter:

        - nr: a unique id, that never should be set/changed by the user.
        - tag: an identification tag used to group properties
        - name: the name to be used for this set. Default is to use an
          automatically generated name.
        - set: identifies the geometrical elements for which the defined
          properties will hold. This can be either:

          - a single number,
          - a list of numbers,
          - the name of an already defined set,
          - a list of such names.

        Besides these, any other fields may be defined and will be added
        without checking.
        """
        d = CDict()
        # update with kargs first, to make sure tag,set and nr are sane
        d.update(dict(**kargs))

        prop = getattr(self,kind+'prop')
        d.nr = len(prop)
        if tag is not None:
            d.tag = str(tag)
        if name is None and 'setname' in kargs:
            # allow for backwards compatibility
            pf.utils.warn("!! 'setname' is deprecated, please use 'name'")
            name = setname
        if name is None and type(set) is str:
            ### convenience to allow set='name' as alias for name='name'
            ### to reuse already defined set
            name,set = set,name
        if name is None:
            name = self.autoName(kind,d.nr)
        elif type(name) is not str:
            raise ValueError,"Property name should be a string"
        d.name = name
        if set is not None:
            if type(set) is int or type(set) is str:
                set = [ set ]
            d.set = unique(set)

        prop.append(d)
        return d


    # This should maybe change to operate on the property keys
    # and finally return the selected keys or properties?

    def getProp(self,kind='',rec=None,tag=None,attr=[],noattr=[],delete=False):
        """Return all properties of type kind matching tag and having attr.

        kind is either '', 'n', 'e' or 'm'
        If rec is given, it is a list of record numbers or a single number.
        If a tag or a list of tags is given, only the properties having a
        matching tag attribute are returned.

        attr and noattr are lists of attributes. Only the properties having
        all the attributes in attr and none of the properties in noattr are
        returned.
        Attributes whose value is None are treated as non-existing.

        If delete==True, the returned properties are removed from the database.
        """
        prop = getattr(self,kind+'prop')
        if rec is not None:
            if type(rec) != list:
                rec = [ rec ]
            rec = [ i for i in rec if i < len(prop) ]
            prop = [ prop[i] for i in rec ]
        if tag is not None:
            if type(tag) != list:
                tag = [ tag ]
            tag = map(str,tag)   # tags are always converted to strings!
            prop = [ p for p in prop if 'tag' in p and p['tag'] in tag ]
        for a in attr:
            prop = [ p for p in prop if a in p and p[a] is not None ]
        for a in noattr:
            prop = [ p for p in prop if a not in p or p[a] is None ]
        if delete:
            self._delete(prop,kind=kind)
        return prop


    def _delete(self,plist,kind=''):
        """Delete the specified properties from the database.

        plist is a list of property records (not copies of them!),
        such as returned by getProp.
        The kind parameter can specify a specific property database.
        """
        prop = getattr(self,kind+'prop')
        if not type(plist) == list:
            pdel = [ plist ]
        for p in plist:
            RemoveListItem(prop,p)
        self._sanitize(kind)


    def _sanitize(self,kind):
        """Sanitize the record numbers after deletion"""
        prop = getattr(self,kind+'prop')
        for i,p in enumerate(prop):
            p.nr = i


    def delProp(self,kind='',rec=None,tag=None,attr=[]):
        """Delete properties.

        This is equivalent to getProp() but the returned properties
        are removed from the database.
        """
        return self.getProp(kind=kind,rec=rec,tag=tag,attr=attr,delete=True)


    def nodeProp(self,prop=None,set=None,name=None,tag=None,cload=None,bound=None,displ=None,veloc=None,accel=None,csys=None,ampl=None,**kargs):
        """Create a new node property, empty by default.

        A node property can contain any combination of the following fields:

        - tag: an identification tag used to group properties (this is e.g.
          used to flag Step, increment, load case, ...)
        - set: a single number or a list of numbers identifying the node(s)
          for which this property will be set, or a set name
          If None, the property will hold for all nodes.
        - cload: a concentrated load: a list of 6 float values
          [FX,FY,FZ,MX,MY,MZ] or a list of (dofid,value) tuples.
        - displ,veloc,accel: prescribed displacement, velocity or
          acceleration: a list of 6 float values
          [UX,UY,UZ,RX,RY,RZ] or  a list of tuples (dofid,value)
        - bound: a boundary condition: a str, a list of 6 codes (0/1), or
          a list of tuples (what??)
        - csys: a CoordSystem
        - ampl: the name of an Amplitude
        """
        try:
            d = kargs
            if cload is not None:
                d['cload'] = checkArrayOrIdValue(cload)
            if displ is not None:
                d['displ'] = checkArrayOrIdValue(displ)
            if veloc is not None:
                d['veloc'] = checkArrayOrIdValue(veloc)
            if accel is not None:
                d['accel'] = checkArrayOrIdValue(accel)
            if bound is not None:
                if type(bound) == str:
                    d['bound'] = checkString(bound,self.bound_strings)
                elif type(bound) == list:
                    if type(bound[0]) != tuple:
                        d['bound'] = checkArray1D(bound,kind='i',size=6)
                    else:
                        d['bound'] = bound # unchecked
            if csys is not None:
                if isinstance(csys,CoordSystem):
                    d['csys'] = csys
                else:
                    raise ValueError,"Invalid Coordinate System"

            # Currently unchecked!
            if ampl is not None:
                d['ampl'] = ampl
            return self.Prop(kind='n',prop=prop,tag=tag,set=set,name=name,**d)
        except:
            print("tag=%s,set=%s,name=%s,cload=%s,bound=%s,displ=%s,csys=%s" % (tag,set,name,cload,bound,displ,csys))
            raise ValueError,"Invalid Node Property"


    def elemProp(self,prop=None,grp=None,set=None,name=None,tag=None,section=None,eltype=None,dload=None,eload=None,ampl=None,**kargs):
        """Create a new element property, empty by default.

        An elem property can contain any combination of the following fields:

        - tag: an identification tag used to group properties (this is e.g.
          used to flag Step, increment, load case, ...)
        - set: a single number or a list of numbers identifying the element(s)
          for which this property will be set, or a set name
          If None, the property will hold for all elements.
        - grp: an elements group number (default None). If specified, the
          element numbers given in set are local to the specified group.
          If not, elements are global and should match the global numbering
          according to the order in which element groups will be specified
          in the Model.
        - eltype: the element type (currently in Abaqus terms).
        - section: an ElemSection specifying the element section properties.
        - dload: an ElemLoad specifying a distributed load on the element.
        - ampl: the name of an Amplitude
        """
        try:
            d = {}
            if eltype is not None:
                d['eltype'] = eltype.upper()
            if section is not None:
                d['section'] = section
            if dload is not None:
                d['dload'] = dload
            if eload is not None:
                d['eload'] = eload
            # Currently unchecked!
            if ampl is not None:
                d['ampl'] = ampl
            d.update(kargs)


            return self.Prop(kind='e',prop=prop,tag=tag,set=set,name=name,**d)
        except:
            raise ValueError,"Invalid Elem Property\n  tag=%s,set=%s,name=%s,eltype=%s,section=%s,dload=%s,eload=%s" % (tag,set,name,eltype,section,dload,eload)


##################################### Test ###########################

if __name__ == "script" or  __name__ == "draw":


    if pf.GUI:
        chdir(__file__)
    print(os.getcwd())

    P = PropertyDB()

    Stick = P.Prop(color='green',name='Stick',weight=25,comment='This could be anything: a gum, a frog, a usb-stick,...')
    print(Stick)

    author = P.Prop(tag='author',alias='Alfred E Neuman',address=CDict({'street':'Krijgslaan', 'city':'Gent','country':'Belgium'}))

    print(P.getProp(tag='author')[0])

    Stick.weight=30
    Stick.length=10
    print(Stick)

    print(author.street)
    author.street='Voskenslaan'
    print(author.street)
    print(author.address.street)
    author.address.street = 'Wiemersdreef'
    print(author.address.street)

    author = P.Prop(tag='author',name='John Doe',address={'city': 'London', 'street': 'Downing Street 10', 'country': 'United Kingdom'})
    print(author)

    for p in P.getProp(rec=[0,2]):
        print(p.name)

    for p in P.getProp(tag=['author']):
        print(p.name)

    for p in P.getProp(attr=['name']):
        print(p.nr)


    P.Prop(set=[0,1,3],name='green_elements',color='green')
    P.Prop(name='green_elements',transparent=True)
    a = P.Prop(set=[0,2,4,6],thickness=3.2)
    P.Prop(name=a.name,material='steel')

    for p in P.getProp(attr=['name']):
        print(p)

    P.Prop(set='green_elements',transparent=False)
    for p in P.getProp(attr=['name']):
        if p.name == 'green_elements':
            print(p.nr,p.transparent)

    print("before")
    for p in P.getProp():
        print(p)

    P.getProp(attr=['transparent'],delete=True)
    P.delProp(attr=['color'])
    pl = P.getProp(rec=[5,6])
    print(pl)
    P._delete(pl)

    print("after")
    for p in P.getProp():
        print(p)

    #exit()

    times = arange(10)
    values = square(times)
    amp = Amplitude(column_stack([times,values]))
    P.Prop(amplitude=amp,name='amp1')

    P1 = [ 1.0,1.0,1.0, 0.0,0.0,0.0 ]
    P2 = [ 0.0 ] * 3 + [ 1.0 ] * 3
    B1 = [ 1 ] + [ 0 ] * 5
    CYL = CoordSystem('cylindrical',[0,0,0,0,0,1])
    # node property on single node
    P.nodeProp(1,cload=[5,0,-75,0,0,0])
    # node property on nodes 2 and 3
    P.nodeProp(set=[2,3],bound='pinned')
    # node property on ALL nodes
    P.nodeProp(cload=P1,bound=B1,csys=CYL,ampl='amp1')
    # node property whose set will be reused
    nset1 = P.nodeProp(tag='step1',set=[2,3,4],cload=P1).nr
    # node properties with an already named set
    P.nodeProp(tag='step2',set=Nset(nset1),cload=P2)

    print('nodeproperties')
    print(P.nprop)

    print('all nodeproperties')
    print(P.getProp('n'))

    print("properties 0 and 2")
    for p in P.getProp('n',rec=[0,2]):
        print(p)

    print("tags 1 and step1")
    for p in P.getProp('n',tag=[1,'step1']):
        print(p)

    print("cload attributes")
    for p in P.getProp('n',attr=['cload']):
        print(p)

    vert = ElemSection('IPEA100', 'steel')
    hor = ElemSection({'name':'IPEM800','A':951247,'I':CDict({'Ix':1542,'Iy':6251,'Ixy':352})}, {'name':'S400','E':210,'fy':400})
    circ = ElemSection({'name':'circle','radius':10,'sectiontype':'circ'},'steel')

    print("Materials")
    for m in Mat:
        print(Mat[m])

    print("Sections")
    for s in Sec:
        print(Sec[s])


    q1 = ElemLoad('PZ',2.5)
    q2 = ElemLoad('PY',3.14)


    top = P.elemProp(set=[0,1,2],eltype='B22',section=hor,dload=q1)
    column = P.elemProp(eltype='B22',section=vert)
    diagonal = P.elemProp(eltype='B22',section=hor)
    bottom = P.elemProp(section=hor,dload=q2,ampl='amp1')


    print('elemproperties')
    for p in P.eprop:
        print(p)

    print("section properties")
    for p in P.getProp('e',attr=['section']):
        print(p.nr)

    P.Prop(set='cylinder',name='cylsurf',surftype='element',label='SNEG')

    print(P.getProp(attr=['surftype']))

# End
