From ab5d463760cd76ad1f24a149cf7f181cc54b82dc Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Thu, 11 Sep 2014 13:08:16 +0200
Subject: Cosmetics.

---
 src/brill.c | 13 ++++++++++---
 1 file changed, 10 insertions(+), 3 deletions(-)

(limited to 'src/brill.c')

diff --git a/src/brill.c b/src/brill.c
index b48980d..cc45587 100644
--- a/src/brill.c
+++ b/src/brill.c
@@ -45,7 +45,7 @@ typedef struct BasisSet {
 /*
  * The basis of even (n = 2 * idx) SB functions (Boyd 2000, Ch 17.9)
  * SB(x, n) = sin((n + 1) arccot(|x| / L))
- * They are symmetric wrt beginning and decay as 1/x in infinity.
+ * They are symmetric wrt origin and decay as 1/x in infinity.
  */
 
 static CCTK_REAL scale_factor;
@@ -232,10 +232,16 @@ static int solve_svd(gsl_matrix *mat, gsl_vector **px, gsl_vector **prhs)
     return 0;
 }
 
-/* solve the equation
+/*
+ * Solve the equation
  * Δψ + ¼ ψ (∂²q/∂r² + ∂²q/∂z²) = 0
  * for the coefficients of spectral approximation of ψ:
- * ψ(r, z) = ΣaᵢⱼTᵢ(r)Tⱼ(z)
+ * ψ(r, z) = 1 + ΣaᵢⱼTᵢ(r)Tⱼ(z)
+ * where i =  { 0, ... , bd->nb_coeffs_x };
+ *       j =  { 0, ... , bd->nb_coeffs_z };
+ * q(r, z) and Tᵢ(x) are defined by bd->q_func, bd->laplace_q_func and
+ * bd->basis.
+ *
  * The cofficients are exported in bd->psi_coeffs
  */
 static int brill_solve(BrillDataContext *bd)
@@ -283,6 +289,7 @@ static int brill_solve(BrillDataContext *bd)
     else
         ret = solve_linear(mat, &coeffs, &rhs);
 
+    /* export the result to the caller */
     bd->psi_coeffs = coeffs;
 
 #if ACC_TEST
-- 
cgit v1.2.3