Taras' Numerical Modelling Course (the Python Way)

Conversion of Taras Matlab models to Python 3

Here you'll find Python 3 translations of some of the models originally developed in Matlab as part of Taras Gerya's introduction to Numerical Modelling (2017).

The aim is to offer an open source alternative to interested students, and compare computation times and code readability between Python and Matlab.

This repository is a work in progress.

Solvers

Global indexing in Python

The basic global indexing relations introduced in Taras course to relate the numerical (grid) and algebraic indices assume that indexing starts at 1. Therefore, they have to adapted because in Python, indexing starts at 0.

So in Matlab,

$gk = (j-1)\times Ny+i$

becomes in Python:

$gk = (j)\times Ny+i$

And similarly,

$gk = ((j-1)\times Ny+i-1)\times U+1$

becomes:

$gk = ((j)\times Ny+i)\times U$

brixel-1984/TarasNumericalModellingCourse

Taras' Numerical Modelling Course (the Python Way)

Solvers

Global indexing in Python

Boundary Conditions in Python