Add support for off-centered waves.

Bump SONAME due to ABI break.

1 files changed, 7 insertions, 1 deletions

diff --git a/brill_data.h b/brill_data.h
index 05a3511..76401bd 100644
--- a/brill_data.h
+++ b/brill_data.h
@@ -43,7 +43,7 @@
 
 enum BDQFuncType {
     /**
-     * q(ρ, z) = A ρ^2 exp(-ρ^2 - z^2)
+     * q(ρ, z) = A ρ^2 exp(-(ρ-ρ_0)^2 - z^2)
      */
     BD_Q_FUNC_GUNDLACH,
     /**
@@ -108,6 +108,12 @@ typedef struct BDContext {
      */
     unsigned int eppley_n;
 
+    /**
+     * Center of the wave for BD_Q_FUNC_GUNDLACH.
+     * Defaults to 0.
+     */
+    double rho0;
+
     /********************************
      *        solver options        *
      ********************************/