#!/bin/bash

if [ ! -e damp.done ] ; then
    mpiexec -np 18 Pelegant damp.ele
    sddsprocess damp.twi -process=s,max,%sMax -define=param,T0,"sMax c_mks /",units=s
fi

# Analytical solution:
# Sdelta^2 = (SdeltaEq)^2 + (SdeltaIni^2 - SdeltaEq^2)*exp(-2*t/tau)

sddssequence -pipe=out -define=Pass,type=long -seq=begin=0,end=30000,delta=1 \
    | sddsxref damp.twi -leave=* -transfer=par,* -pipe \
    | sddsprocess -pipe=in analytical.sdds \
    "-define=param,SdeltaIni,4e-3" \
    "-define=col,Sdelta,Sdelta0 sqr SdeltaIni sqr Sdelta0 sqr - 2 T0 * taudelta / Pass * exp / + sqrt"

# Compare

sddsplot  -col=Pass,Sdelta damp001.w -graph=dot -col=Pass,Sdelta analytical.sdds -graph=line,type=1


