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

.. only:: html

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

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

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

.. _sphx_glr__gallery_mesh_surface_gaussian.py:


Gaussian kernel as a surface
============================

Example showing a gaussian kernel as a surface mesh

.. GENERATED FROM PYTHON SOURCE LINES 7-45



.. image-sg:: /_gallery/mesh/images/sphx_glr_surface_gaussian_001.webp
   :alt: surface gaussian
   :srcset: /_gallery/mesh/images/sphx_glr_surface_gaussian_001.webp
   :class: sphx-glr-single-img


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

 .. code-block:: none

    /opt/hostedtoolcache/Python/3.12.12/x64/lib/python3.12/site-packages/pygfx/objects/_ruler.py:400: RuntimeWarning: divide by zero encountered in divide
      screen_full = (ndc_full[:, :2] / ndc_full[:, 3:4]) * half_canvas_size
    /opt/hostedtoolcache/Python/3.12.12/x64/lib/python3.12/site-packages/pygfx/objects/_ruler.py:400: RuntimeWarning: invalid value encountered in divide
      screen_full = (ndc_full[:, :2] / ndc_full[:, 3:4]) * half_canvas_size
    /opt/hostedtoolcache/Python/3.12.12/x64/lib/python3.12/site-packages/pygfx/objects/_ruler.py:412: RuntimeWarning: invalid value encountered in divide
      screen_sel = (ndc_sel[:, :2] / ndc_sel[:, 3:4]) * half_canvas_size
    /home/runner/work/fastplotlib/fastplotlib/fastplotlib/graphics/features/_base.py:18: UserWarning: casting float64 array to float32
      warn(f"casting {array.dtype} array to float32")







|

.. code-block:: Python


    # test_example = true

    import fastplotlib as fpl
    import numpy as np


    figure = fpl.Figure(size=(700, 560), cameras="3d", controller_types="orbit")


    def gaus2d(x=0, y=0, mx=0, my=0, sx=1, sy=1):
        return (
            1.0
            / (2.0 * np.pi * sx * sy)
            * np.exp(
                -((x - mx) ** 2.0 / (2.0 * sx**2.0) + (y - my) ** 2.0 / (2.0 * sy**2.0))
            )
        )


    r = np.linspace(0, 10, num=200)
    x, y = np.meshgrid(r, r)
    z = gaus2d(x, y, mx=5, my=5, sx=1, sy=1) * 50

    mesh = figure[0, 0].add_surface(
        np.dstack([x, y, z]), mode="phong", cmap="jet"
    )

    # figure[0, 0].axes.grids.xy.visible = True
    figure[0, 0].camera.show_object(mesh.world_object, (-2, 2, -2), up=(0, 0, 1))
    figure.show()


    # NOTE: fpl.loop.run() should not be used for interactive sessions
    # See the "JupyterLab and IPython" section in the user guide
    if __name__ == "__main__":
        print(__doc__)
        fpl.loop.run()


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

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


.. _sphx_glr_download__gallery_mesh_surface_gaussian.py:

.. only:: html

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

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

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

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

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


.. only:: html

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

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