.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/thirdorderblwdf.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_thirdorderblwdf.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_thirdorderblwdf.py:


=============================
Third-order Bireciprocal LWDF
=============================

Small bireciprocal lattice wave digital filter.

.. GENERATED FROM PYTHON SOURCE LINES 8-29

.. code-block:: Python


    import numpy as np
    from mplsignal.freq_plots import freqz_fir

    from b_asic.core_operations import Addition, SymmetricTwoportAdaptor
    from b_asic.schedule import Schedule
    from b_asic.signal_flow_graph import SFG
    from b_asic.signal_generator import Impulse
    from b_asic.simulation import Simulation
    from b_asic.special_operations import Delay, Input, Output

    in0 = Input("x")
    D0 = Delay(in0)
    D1 = Delay()
    D2 = Delay(D1)
    s = SymmetricTwoportAdaptor(-0.375, in0, D2)
    D1 <<= s.output(1)
    a = s.output(0) + D0
    out0 = Output(a, "y")

    sfg = SFG(inputs=[in0], outputs=[out0], name="Third-order BLWDF")







.. GENERATED FROM PYTHON SOURCE LINES 30-31

The SFG looks like

.. GENERATED FROM PYTHON SOURCE LINES 31-33

.. code-block:: Python

    sfg






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <?xml version="1.0" encoding="UTF-8" standalone="no"?>
    <!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"
     "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">
    <!-- Generated by graphviz version 2.43.0 (0)
     -->
    <!-- Title: %3 Pages: 1 -->
    <svg width="431pt" height="122pt"
     viewBox="0.00 0.00 431.09 122.00" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
    <g id="graph0" class="graph" transform="scale(1 1) rotate(0) translate(4 118)">
    <title>%3</title>
    <polygon fill="white" stroke="transparent" points="-4,4 -4,-118 427.09,-118 427.09,4 -4,4"/>
    <!-- in0 -->
    <g id="node1" class="node">
    <title>in0</title>
    <polygon fill="none" stroke="black" points="42,-81 0,-81 0,-55 42,-55 54,-68 42,-81"/>
    <text text-anchor="middle" x="27" y="-71.8" font-family="Times,serif" font-size="14.00">x</text>
    <text text-anchor="middle" x="27" y="-56.8" font-family="Times,serif" font-size="14.00">(in0)</text>
    </g>
    <!-- in0.0 -->
    <g id="node2" class="node">
    <title>in0.0</title>
    <ellipse fill="black" stroke="black" cx="91.8" cy="-68" rx="1.8" ry="1.8"/>
    </g>
    <!-- in0&#45;&gt;in0.0 -->
    <g id="edge2" class="edge">
    <title>in0&#45;&gt;in0.0</title>
    <path fill="none" stroke="black" d="M54.27,-68C69.33,-68 86,-68 89.94,-68"/>
    </g>
    <!-- t0 -->
    <g id="node3" class="node">
    <title>t0</title>
    <polygon fill="none" stroke="black" points="195.05,-113 159.05,-113 159.05,-77 195.05,-77 195.05,-113"/>
    <text text-anchor="middle" x="177.05" y="-91.3" font-family="Times,serif" font-size="14.00">t0</text>
    </g>
    <!-- in0.0&#45;&gt;t0 -->
    <g id="edge1" class="edge">
    <title>in0.0&#45;&gt;t0</title>
    <path fill="none" stroke="black" d="M93.73,-68.3C99.16,-70.06 127.14,-79.14 148.99,-86.23"/>
    <polygon fill="black" stroke="black" points="148.21,-89.65 158.8,-89.41 150.36,-82.99 148.21,-89.65"/>
    </g>
    <!-- sym2p0 -->
    <g id="node4" class="node">
    <title>sym2p0</title>
    <ellipse fill="none" stroke="black" cx="177.05" cy="-41" rx="47.39" ry="18"/>
    <text text-anchor="middle" x="177.05" y="-37.3" font-family="Times,serif" font-size="14.00">sym2p0</text>
    </g>
    <!-- in0.0&#45;&gt;sym2p0 -->
    <g id="edge3" class="edge">
    <title>in0.0&#45;&gt;sym2p0</title>
    <path fill="none" stroke="black" d="M93.73,-67.7C97.68,-66.42 113.51,-61.28 130.26,-55.85"/>
    <polygon fill="black" stroke="black" points="131.71,-59.06 140.14,-52.65 129.55,-52.4 131.71,-59.06"/>
    <text text-anchor="middle" x="135.14" y="-56.45" font-family="Times,serif" font-size="14.00">0</text>
    </g>
    <!-- add0 -->
    <g id="node6" class="node">
    <title>add0</title>
    <ellipse fill="none" stroke="black" cx="294.29" cy="-95" rx="33.6" ry="18"/>
    <text text-anchor="middle" x="294.29" y="-91.3" font-family="Times,serif" font-size="14.00">add0</text>
    </g>
    <!-- t0&#45;&gt;add0 -->
    <g id="edge5" class="edge">
    <title>t0&#45;&gt;add0</title>
    <path fill="none" stroke="black" d="M195.05,-95C209.59,-95 231.16,-95 250.38,-95"/>
    <polygon fill="black" stroke="black" points="250.46,-98.5 260.46,-95 250.46,-91.5 250.46,-98.5"/>
    <text text-anchor="middle" x="255.46" y="-83.8" font-family="Times,serif" font-size="14.00">1</text>
    </g>
    <!-- sym2p0&#45;&gt;add0 -->
    <g id="edge6" class="edge">
    <title>sym2p0&#45;&gt;add0</title>
    <path fill="none" stroke="black" d="M207.72,-54.91C223.44,-62.27 242.83,-71.36 259.17,-79.01"/>
    <polygon fill="black" stroke="black" points="257.85,-82.26 268.39,-83.33 260.82,-75.92 257.85,-82.26"/>
    <text text-anchor="middle" x="273.39" y="-72.13" font-family="Times,serif" font-size="14.00">0</text>
    <text text-anchor="middle" x="202.72" y="-58.71" font-family="Times,serif" font-size="14.00">0</text>
    </g>
    <!-- t2 -->
    <g id="node8" class="node">
    <title>t2</title>
    <polygon fill="none" stroke="black" points="312.29,-59 276.29,-59 276.29,-23 312.29,-23 312.29,-59"/>
    <text text-anchor="middle" x="294.29" y="-37.3" font-family="Times,serif" font-size="14.00">t2</text>
    </g>
    <!-- sym2p0&#45;&gt;t2 -->
    <g id="edge8" class="edge">
    <title>sym2p0&#45;&gt;t2</title>
    <path fill="none" stroke="black" d="M224.59,-41C238.7,-41 253.72,-41 266.12,-41"/>
    <polygon fill="black" stroke="black" points="266.24,-44.5 276.24,-41 266.24,-37.5 266.24,-44.5"/>
    <text text-anchor="middle" x="229.59" y="-29.8" font-family="Times,serif" font-size="14.00">1</text>
    </g>
    <!-- out0 -->
    <g id="node5" class="node">
    <title>out0</title>
    <polygon fill="none" stroke="black" points="411.09,-108 364.09,-108 364.09,-82 411.09,-82 423.09,-95 411.09,-108"/>
    <text text-anchor="middle" x="393.59" y="-98.8" font-family="Times,serif" font-size="14.00">y</text>
    <text text-anchor="middle" x="393.59" y="-83.8" font-family="Times,serif" font-size="14.00">(out0)</text>
    </g>
    <!-- add0&#45;&gt;out0 -->
    <g id="edge4" class="edge">
    <title>add0&#45;&gt;out0</title>
    <path fill="none" stroke="black" d="M328.23,-95C336.46,-95 345.36,-95 353.85,-95"/>
    <polygon fill="black" stroke="black" points="353.91,-98.5 363.91,-95 353.91,-91.5 353.91,-98.5"/>
    </g>
    <!-- t1 -->
    <g id="node7" class="node">
    <title>t1</title>
    <polygon fill="none" stroke="black" points="411.59,-36 375.59,-36 375.59,0 411.59,0 411.59,-36"/>
    <text text-anchor="middle" x="393.59" y="-14.3" font-family="Times,serif" font-size="14.00">t1</text>
    </g>
    <!-- t1&#45;&gt;sym2p0 -->
    <g id="edge7" class="edge">
    <title>t1&#45;&gt;sym2p0</title>
    <path fill="none" stroke="black" d="M375.35,-16.64C362.37,-15.69 344.12,-14.5 328.09,-14 298.06,-13.07 290.16,-9.27 260.49,-14 246.86,-16.18 232.42,-20.25 219.51,-24.59"/>
    <polygon fill="black" stroke="black" points="218.24,-21.33 209.96,-27.93 220.56,-27.93 218.24,-21.33"/>
    <text text-anchor="middle" x="204.96" y="-16.73" font-family="Times,serif" font-size="14.00">1</text>
    </g>
    <!-- t2&#45;&gt;t1 -->
    <g id="edge9" class="edge">
    <title>t2&#45;&gt;t1</title>
    <path fill="none" stroke="black" d="M312.68,-36.89C327.35,-33.42 348.59,-28.4 365.48,-24.41"/>
    <polygon fill="black" stroke="black" points="366.41,-27.78 375.34,-22.08 364.8,-20.97 366.41,-27.78"/>
    </g>
    </g>
    </svg>

    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 34-35

Set latencies and execution times

.. GENERATED FROM PYTHON SOURCE LINES 35-40

.. code-block:: Python

    sfg.set_latency_of_type(SymmetricTwoportAdaptor, 4)
    sfg.set_latency_of_type(Addition, 1)
    sfg.set_execution_time_of_type(SymmetricTwoportAdaptor, 1)
    sfg.set_execution_time_of_type(Addition, 1)








.. GENERATED FROM PYTHON SOURCE LINES 41-42

Simulate

.. GENERATED FROM PYTHON SOURCE LINES 42-46

.. code-block:: Python

    sim = Simulation(sfg, [Impulse()])
    sim.run_for(1000)






.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    [0.0]



.. GENERATED FROM PYTHON SOURCE LINES 47-48

Display output

.. GENERATED FROM PYTHON SOURCE LINES 48-50

.. code-block:: Python

    freqz_fir(np.array(sim.results['0']) / 2)




.. image-sg:: /examples/images/sphx_glr_thirdorderblwdf_001.png
   :alt: thirdorderblwdf
   :srcset: /examples/images/sphx_glr_thirdorderblwdf_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <Figure size 640x480 with 2 Axes>



.. GENERATED FROM PYTHON SOURCE LINES 51-52

Create and display schedule

.. GENERATED FROM PYTHON SOURCE LINES 52-54

.. code-block:: Python

    schedule = Schedule(sfg, cyclic=True)
    schedule.show()



.. image-sg:: /examples/images/sphx_glr_thirdorderblwdf_002.png
   :alt: thirdorderblwdf
   :srcset: /examples/images/sphx_glr_thirdorderblwdf_002.png
   :class: sphx-glr-single-img






.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 1.225 seconds)


.. _sphx_glr_download_examples_thirdorderblwdf.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: thirdorderblwdf.ipynb <thirdorderblwdf.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: thirdorderblwdf.py <thirdorderblwdf.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: thirdorderblwdf.zip <thirdorderblwdf.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_