Minimization of Nonlinear Energies in Python Using FEM and Automatic Differentiation Tools

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Béreš, Michal, Valdman, Jan
Format: Preprint
Publié: 2024
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866911946363633664
author Béreš, Michal
Valdman, Jan
author_facet Béreš, Michal
Valdman, Jan
contents This contribution examines the capabilities of the Python ecosystem to solve nonlinear energy minimization problems, with a particular focus on transitioning from traditional MATLAB methods to Python's advanced computational tools, such as automatic differentiation. We demonstrate Python's streamlined approach to minimizing nonlinear energies by analyzing three problem benchmarks - the p-Laplacian, the Ginzburg-Landau model, and the Neo-Hookean hyperelasticity. This approach merely requires the provision of the energy functional itself, making it a simple and efficient way to solve this category of problems. The results show that the implementation is about ten times faster than the MATLAB implementation for large-scale problems. Our findings highlight Python's efficiency and ease of use in scientific computing, establishing it as a preferable choice for implementing sophisticated mathematical models and accelerating the development of numerical simulations.
format Preprint
id arxiv_https___arxiv_org_abs_2407_04706
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Minimization of Nonlinear Energies in Python Using FEM and Automatic Differentiation Tools
Béreš, Michal
Valdman, Jan
Mathematical Software
This contribution examines the capabilities of the Python ecosystem to solve nonlinear energy minimization problems, with a particular focus on transitioning from traditional MATLAB methods to Python's advanced computational tools, such as automatic differentiation. We demonstrate Python's streamlined approach to minimizing nonlinear energies by analyzing three problem benchmarks - the p-Laplacian, the Ginzburg-Landau model, and the Neo-Hookean hyperelasticity. This approach merely requires the provision of the energy functional itself, making it a simple and efficient way to solve this category of problems. The results show that the implementation is about ten times faster than the MATLAB implementation for large-scale problems. Our findings highlight Python's efficiency and ease of use in scientific computing, establishing it as a preferable choice for implementing sophisticated mathematical models and accelerating the development of numerical simulations.
title Minimization of Nonlinear Energies in Python Using FEM and Automatic Differentiation Tools
topic Mathematical Software
url https://arxiv.org/abs/2407.04706