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

.. only:: html

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

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

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

.. _sphx_glr_examples_gallery_advanced_pdes_heterogeneous_pde.py:


Heterogeneous PDE
=================

This example loads an example image and uses it as the source field for a simple
reaction-diffusion equation.

.. GENERATED FROM PYTHON SOURCE LINES 8-25



.. image-sg:: /examples_gallery/advanced_pdes/images/sphx_glr_heterogeneous_pde_001.png
   :alt: heterogeneous pde
   :srcset: /examples_gallery/advanced_pdes/images/sphx_glr_heterogeneous_pde_001.png
   :class: sphx-glr-single-img


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

 .. code-block:: none


      0%|          | 0/100.0 [00:00<?, ?it/s]
    Initializing:   0%|          | 0/100.0 [00:00<?, ?it/s]
      0%|          | 0/100.0 [00:00<?, ?it/s]              
      0%|          | 0.27743/100.0 [00:08<50:20, 30.29s/it]
      0%|          | 0.37315/100.0 [00:08<37:23, 22.52s/it]
      2%|▏         | 2.45746/100.0 [00:08<05:34,  3.42s/it]
     23%|██▎       | 23.29858/100.0 [00:08<00:27,  2.76it/s]
     23%|██▎       | 23.29858/100.0 [00:08<00:28,  2.73it/s]
    100%|██████████| 100.0/100.0 [00:08<00:00, 11.72it/s]   
    100%|██████████| 100.0/100.0 [00:08<00:00, 11.72it/s]






|

.. code-block:: Python


    import inspect
    from pathlib import Path

    from pde import PDE, ScalarField

    # load a field relative to the current file
    package_path = Path(inspect.getfile(lambda: None)).parents[2]
    img_path = package_path / "docs" / "source" / "_images" / "logo_small.png"

    background = ScalarField.from_image(img_path)  # create source field from image
    state = ScalarField(background.grid)  # generate initial condition

    # define the pde
    eq = PDE({"c": "laplace(c) + 0.2 * source - 0.1 * c"}, consts={"source": background})
    result = eq.solve(state, t_range=100, adaptive=True)
    result.plot()


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

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


.. _sphx_glr_download_examples_gallery_advanced_pdes_heterogeneous_pde.py:

.. only:: html

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

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

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

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

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

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

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