#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#    ____       ____  _      _
#   |  _ \ __ _/ ___|| |_   / \
#   | |_) / _` \___ \| __| / _ \
#   |  __/ (_| |___) | |_ / ___ \
#   |_|   \__,_|____/ \__/_/   \_\
#
#   a software for Parallel Stability Analysis
#
# Created by Jacopo Canton under the supervision of Franco Auteri and Marco
# Carini, based on Luigi Quartapelle's subroutines for FEM discretization.
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
# FEATURES:
#  - axisymmetric base flow computation
#  - continuation of solutions and bifurcation tracking (LOCA90)
#  - eigenvalues and eigenvectors computation -3d- (ARPACK library)
#  - structural sensitivity computation -3d-
#  - computation of the optimal linear perturbations -3d- (transient growth)
#  - DNS
#  - evolution of optimal linear perturbations on the nonlinear system
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
# CONTACTS:
#  jcanton@mech.kth.se
#  franco.auteri@polimi.it
#  marco.carini@polimi.it
#  luigi.quartapelle@polimi.it
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~