ALLEN INSTITUTE

# Source code for bmtk.simulator.bionet.default_setters.synaptic_weights

    # Copyright 2017. Allen Institute. All rights reserved
    #
    # Redistribution and use in source and binary forms, with or without modification, are permitted provided that the
    # following conditions are met:
    #
    # 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following
    # disclaimer.
    #
    # 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following
    # disclaimer in the documentation and/or other materials provided with the distribution.
    #
    # 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote
    # products derived from this software without specific prior written permission.
    #
    # THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES,
    # INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
    # DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
    # SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
    # SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
    # WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
    # OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
    #
    import math

    from bmtk.simulator.bionet.pyfunction_cache import add_weight_function

    def default_weight_fnc(edge_props, src_props, trg_props):
        return edge_props['syn_weight']

    def wmax(edge_props, src_props, trg_props):
        return edge_props["syn_weight"]

    def gaussianLL(edge_props, src_props, trg_props):
        src_tuning = src_props['tuning_angle']
        tar_tuning = trg_props['tuning_angle']

        w0 = edge_props["syn_weight"]
        sigma = edge_props["weight_sigma"]

        delta_tuning = abs(abs(abs(180.0 - abs(float(tar_tuning) - float(src_tuning)) % 360.0) - 90.0) - 90.0)
        weight = w0 * math.exp(-(delta_tuning / sigma) ** 2)

        return weight

    add_weight_function(wmax, 'wmax', overwrite=False)
    add_weight_function(gaussianLL, 'gaussianLL', overwrite=False)
    add_weight_function(default_weight_fnc, 'default_weight_fnc', overwrite=False)
    add_weight_function(default_weight_fnc, 'set_syn_weight', overwrite=False)