#  Copyright (C) 2012,2013
#      Max Planck Institute for Polymer Research
#  Copyright (C) 2008,2009,2010,2011
#      Max-Planck-Institute for Polymer Research & Fraunhofer SCAI
#
#  This file is part of ESPResSo++.
#
#  ESPResSo++ 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.
#
#  ESPResSo++ 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/>.


r"""
*******************************
espressopp.FixedTripleAngleList
*******************************


.. function:: espressopp.FixedTripleAngleList(storage)

                :param storage:
                :type storage:

.. function:: espressopp.FixedTripleAngleList.add(pid1, pid2, pid3)

                :param pid1:
                :param pid2:
                :param pid3:
                :type pid1:
                :type pid2:
                :type pid3:
                :rtype:

.. function:: espressopp.FixedTripleAngleList.addTriples(triplelist)

                :param triplelist:
                :type triplelist:
                :rtype:

.. function:: espressopp.FixedTripleAngleList.getAngle(pid1, pid2, pid3)

                :param pid1:
                :param pid2:
                :param pid3:
                :type pid1:
                :type pid2:
                :type pid3:
                :rtype:

.. function:: espressopp.FixedTripleAngleList.getTriples()

                :rtype:

.. function:: espressopp.FixedTripleAngleList.getTriplesAngles()

                :rtype:

.. function:: espressopp.FixedTripleAngleList.size()

                :rtype:
"""
from espressopp import pmi
import _espressopp
#import espressopp
from espressopp.esutil import cxxinit

class FixedTripleAngleListLocal(_espressopp.FixedTripleAngleList):


    def __init__(self, storage):

        #if pmi.workerIsActive():
        if not (pmi._PMIComm and pmi._PMIComm.isActive()) or pmi._MPIcomm.rank in pmi._PMIComm.getMPIcpugroup():
            cxxinit(self, _espressopp.FixedTripleAngleList, storage)

    def add(self, pid1, pid2, pid3):

        if pmi.workerIsActive():
            return self.cxxclass.add(self, pid1, pid2, pid3)

    def addTriples(self, triplelist):
        """
        Each processor takes the broadcasted triplelist and
        adds those triples whose first particle is owned by
        this processor.
        """
        if pmi.workerIsActive():
            for triple in triplelist:
                pid1, pid2, pid3 = triple
                self.cxxclass.add(self, pid1, pid2, pid3)

    def size(self):

        if pmi.workerIsActive():
            return self.cxxclass.size(self)

    def getTriples(self):

        if pmi.workerIsActive():
            triples = self.cxxclass.getTriples(self)
            return triples

    'returns the list of (pid1, pid2, pid3, angle(123))'
    def getTriplesAngles(self):

        if pmi.workerIsActive():
            triples_angles = self.cxxclass.getTriplesAngles(self)
            return triples_angles

    def getAngle(self, pid1, pid2, pid3):
        if pmi.workerIsActive():
            return self.cxxclass.getAngle(self, pid1, pid2, pid3)

if pmi.isController:
    class FixedTripleAngleList(metaclass=pmi.Proxy):
        pmiproxydefs = dict(
            cls = 'espressopp.FixedTripleAngleListLocal',
            localcall = [ "add" ],
            pmicall = [ "addTriples" ],
            pmiinvoke = ["getTriples", "getTriplesAngles", "size"]
        )

    def getAngle(self, pid1, pid2, pid3 ):
        angles = pmi.invoke(self.pmiobject, 'getAngle', pid1, pid2, pid3 )
        for i in angles:
            if( i != -1 ):
                return i
