From a9de4a2b5e8e2a38ea793a4724992064102696a1 Mon Sep 17 00:00:00 2001
From: Hagen Wierstorf <hagenw@posteo.de>
Date: Mon, 18 Mar 2019 12:30:10 +0100
Subject: [PATCH] Update links to match correct equation

---
 SFS_general/secondary_source_selection.m                  | 8 ++++----
 .../driving_function_mono_nfchoa_ls.m                     | 6 +-----
 .../driving_function_mono_nfchoa_ps.m                     | 4 ++--
 .../driving_function_mono_nfchoa_pw.m                     | 6 +++---
 .../driving_function_mono_sdm_kx_ps.m                     | 2 +-
 .../driving_function_mono_sdm_kx_pw.m                     | 2 +-
 .../driving_functions_mono/driving_function_mono_sdm_pw.m | 2 +-
 .../driving_functions_mono/driving_function_mono_wfs_fs.m | 6 ++----
 .../driving_functions_mono/driving_function_mono_wfs_ls.m | 7 +++----
 .../driving_functions_mono/driving_function_mono_wfs_ps.m | 8 ++++----
 .../driving_functions_mono/driving_function_mono_wfs_pw.m | 4 ++--
 SFS_monochromatic/greens_function_mono.m                  | 8 +++-----
 SFS_monochromatic/sound_field_mono.m                      | 2 +-
 .../driving_functions_imp/driving_function_imp_wfs_fs.m   | 6 +++---
 .../driving_functions_imp/driving_function_imp_wfs_ls.m   | 5 ++---
 .../driving_functions_imp/driving_function_imp_wfs_ps.m   | 4 ++--
 .../driving_functions_imp/driving_function_imp_wfs_pw.m   | 4 ++--
 SFS_time_domain/greens_function_imp.m                     | 8 +++-----
 SFS_time_domain/sound_field_imp.m                         | 2 +-
 SFS_time_domain/wfs_fir_prefilter.m                       | 4 ++--
 20 files changed, 43 insertions(+), 55 deletions(-)

diff --git a/SFS_general/secondary_source_selection.m b/SFS_general/secondary_source_selection.m
index e2390e7e..8867f1c2 100644
--- a/SFS_general/secondary_source_selection.m
+++ b/SFS_general/secondary_source_selection.m
@@ -98,7 +98,7 @@
     % a = <
     %      \ 0, else
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-wfs.pw.selection
+    % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-wfs-secondary-source-selection-plane
     %
     % Direction of plane wave (nk) is set above
     idx = nx0*nk(:) >= eps;
@@ -111,7 +111,7 @@
     % a = <
     %      \ 0, else
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-wfs.ps.selection
+    % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-wfs-secondary-source-selection-point
     %
     idx = sum(nx0.*x0,2) - nx0*xs(1:3).' >= -2*eps;
 
@@ -126,7 +126,7 @@
     % where v = x0-xs - <x0-xs,nxs > nxs,
     % and |nxs| = 1.
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-wfs.ls.selection
+    % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-wfs-secondary-source-selection-line
     %
     %NOTE: We don't check if we are in a 2D or 3D scenario and use xs(4:6)
     % whenever it is present. This can only provide problems if you use the
@@ -152,7 +152,7 @@
     % a = <
     %      \ 0, else
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-wfs.fs.selection
+    % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-wfs-secondary-source-selection-focused
     %
     nxs = xs(4:6);  % vector for orientation of focused source
     xs = xs(1:3);  % vector for position of focused source
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_ls.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_ls.m
index 70169c9d..07204c46 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_ls.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_ls.m
@@ -97,7 +97,7 @@
         %             2pi r0  m=-N   (2)
         %                           Hm(w/c r0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.nfchoa.ls.2D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-nfchoa-line-2d
         %
         for m=-N:N
             D = D + 1 ./ (2.*pi.*r0) ...
@@ -126,8 +126,6 @@
         %             2r0  m=-N        (2)
         %                         w/c h|m|(w/c r0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.nfchoa.ls.2.5D
-        %
         for m=-N:N
             D = D + 1 ./ (2.*r0) ...
                 .* win(abs(m)+1) ...
@@ -166,8 +164,6 @@
         %                m
         %             x Yn(beta0,alpha0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.nfchoa.ls.3D
-        %
         for n=-N:N
             for m=-n:n
                 D = D + 1 ./ (2.*r0.^2) ...
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_ps.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_ps.m
index b643e534..86f4de4a 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_ps.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_ps.m
@@ -113,7 +113,7 @@
         %             2pi r0  m=-N   (2)
         %                           h|m|(w/c r0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.nfchoa.ps.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-nfchoa-point-25d
         for m=-N:N
             D = D + 1 ./ (2.*pi.*r0) ...
                 .* win(abs(m)+1) ...
@@ -144,7 +144,7 @@
         %                       m
         %                    x Yn(theta0,phi0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.nfchoa.ps.3D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-nfchoa-point-3d
         %
         for n=0:N
             for m=-n:n
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_pw.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_pw.m
index 37a2e028..412d7788 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_pw.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_nfchoa_pw.m
@@ -93,7 +93,7 @@
         %               pi r0  m=-N  (2)
         %                           Hm(w/c r0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.nfchoa.pw.2D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-nfchoa-plane-2d
         %
         for m=-N:N
             D = D - 2.*1i ./ (pi.*r0) ...
@@ -122,7 +122,7 @@
         %                   r0  m=-N       (2)
         %                             i w/c h|m|(w/c r0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.nfchoa.pw.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-nfchoa-plane-25d
         %
         for m=-N:N
             D = D - 2./r0 ...
@@ -154,7 +154,7 @@
         %                       m
         %                    x Yn(theta0,phi0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.nfchoa.pw.3D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-nfchoa-plane-3d
         %
         for n=0:N
             for m=-n:n
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_kx_ps.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_kx_ps.m
index 02b04cae..0bbe0832 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_kx_ps.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_kx_ps.m
@@ -108,7 +108,7 @@
         %                     \ ----------_-_-_-_-_-_---------,       |kx|>|w/c|
         %                          K1( \|kx^2-(w/c)^2 yref )
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.sdm.ps.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-sdm-point-25d
         %
         D(idxpr) =  exp(1i*kx(idxpr)*xs(1)) .* ...
             besselh(1,2,sqrt( (omega/c)^2 - kx(idxpr).^2 )*abs(xref(2)-xs(2))) ./ ...
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_kx_pw.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_kx_pw.m
index 9eca647f..a84d6b86 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_kx_pw.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_kx_pw.m
@@ -100,7 +100,7 @@
         %                    H0  | - nky*xrefy |
         %                         \c          /
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.sdm.pw.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-sdm-plane-25d
         %
         idx = find(kx>=omega/c*nk(:,1),1,'first');
         D(idx) = 4*1i*exp(-1i*omega/c*nk(2).*xref(2)) / ...
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_pw.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_pw.m
index cc12b502..179d9481 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_pw.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_sdm_pw.m
@@ -103,7 +103,7 @@
         %                    H0  | - nky*xrefy |
         %                         \c          /
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.sdm.pw.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_nfchoa/#equation-fd-sdm-plane-25d
         %
         D = 4*1i.*exp(-1i*omega/c.*nk(:,2).*xref(:,2)) ./ ...
             besselh(0,2,omega/c.*nk(:,2).*xref(:,2)) .* ...
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_fs.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_fs.m
index 70c7ef0d..8021ced2 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_fs.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_fs.m
@@ -107,7 +107,7 @@
         % D(x0,w) = --- --- ----------- e^(i w/c |x0-xs|)
         %           2pi  c   |x0-xs|^2
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.fs
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-focused-3d
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -141,8 +141,6 @@
         % D(x0,w) =  - -- --------- H1  | - |x0-xs| |
         %              2c  |x0-xs|      \ c         /
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.fs.ls
-        %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
         % Driving signal
@@ -199,7 +197,7 @@
         % D_2.5D(x0,w) = g0  _ |---  _ |--- ------------- e^(i w/c |x0-xs|)
         %                     \|2pi   \| c  |x0-xs|^(3/2)
         %
-        % See https://sfs.rtfd.io/en/3.2/en/update_wfs_ps/#equation-D.wfs.fs.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-focused-25d
         %
         % Driving signal
         D = 1./sqrt(2.*pi) .* sqrt(-1i.*omega./c) .* g0 ...
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ls.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ls.m
index 8da73774..540e0b0f 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ls.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ls.m
@@ -85,7 +85,7 @@
         % D(x0,w) =  - -- -----------  H1  | - |x0-xs| |
         %              2c   |x0-xs|        \ c         /
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.ls
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-line
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -122,7 +122,7 @@
         % D_2.5D(x0,w) =  - -- _ |---  -----------  H1  | - |x0-xs| |
         %                   2   \| c    |x0-xs|         \ c         /
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.ls.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-line-25d
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -152,8 +152,7 @@
         % where v = x0-xs - <x0-xs,nxs > nxs,
         % and |nxs| = 1.
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.ls
-        % and https://sfs.rtfd.io/en/3.2/#equation-v.ls
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-line
         %
         % v = (I - nxs'nxs)(x0-xs)
         % r = |v|
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ps.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ps.m
index a0382543..b7a7ce9f 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ps.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ps.m
@@ -94,7 +94,7 @@
         % D(x0,w) = --- --- ----------- e^(-i w/c |x0-xs|)
         %           2pi  c   |x0-xs|^2
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.ps
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-point
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -110,7 +110,7 @@
         % D(x0,w) = --- | --- - ------- |  ----------- e^(-i w/c |x0-xs|)
         %           2pi \  c    |x0-xs| /   |x0-xs|^2
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.ps.woapprox
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-point-woapprox
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -186,7 +186,7 @@
         % g0 _ |--- _ |--- ------------- e^(-i w/c |x0-xs|)
         %     \|2pi  \| c  |x0-xs|^(3/2)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.ps.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-point-25d
         %
         % Driving signal
         D = 1./sqrt(2*pi) .* sqrt(1i.*omega./c) .* g0 ...
@@ -214,7 +214,7 @@
         % g0 _ |--- _ |--- ------------- e^(-i w/c |x0-xs|)
         %     \|2pi  \| c  |x0-xs|^(3/2)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.ps.2.5D.refline
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-point-25d-refline
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
diff --git a/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_pw.m b/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_pw.m
index 6d2a1e51..896ae1a5 100644
--- a/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_pw.m
+++ b/SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_pw.m
@@ -86,7 +86,7 @@
         % D(x0,w) =  2 --- nk nx0  e^(-i w/c nk x0)
         %               c
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.pw
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-plane
         %
         D = 2.*1i.*omega./c .* vector_product(nk,nx0,2) ...
             .* exp(-1i.*omega./c.*vector_product(nk,x0,2));
@@ -119,7 +119,7 @@
         % D_2.5D(x0,w) = 2g0 nk nx0 _ |---  e^(-i w/c nk x0)
         %                            \| c
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-D.wfs.pw.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-fd-wfs-plane-25d
         %
         D = 2.*g0 .* vector_product(nk,nx0,2) .* sqrt(1i.*omega./c) ...
             .* exp(-1i.*omega./c.*vector_product(nk,x0,2));
diff --git a/SFS_monochromatic/greens_function_mono.m b/SFS_monochromatic/greens_function_mono.m
index 9f180bf9..9c341c48 100644
--- a/SFS_monochromatic/greens_function_mono.m
+++ b/SFS_monochromatic/greens_function_mono.m
@@ -74,7 +74,7 @@
     % G(x-xs,w) = --- -----------------
     %             4pi      |x-xs|
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-S.ps
+    % https://sfs.rtfd.io/en/3.2/sources/#equation-fd-point
     %
     G = 1/(4*pi) * exp(-1i*omega/c .* sqrt((x-xs(1)).^2+(y-xs(2)).^2+(z-xs(3)).^2)) ./ ...
             sqrt((x-xs(1)).^2+(y-xs(2)).^2+(z-xs(3)).^2);
@@ -86,8 +86,6 @@
     % ---- G(x-xs,w) = --- | ----- + ------- | ----------- e^(-i w/c |x-xs|)
     % d ns             4pi  \  c     |x-xs| /   |x-xs|^2
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-S.dps
-    %
     % r = |x-xs|
     r = sqrt((x-xs(1)).^2+(y-xs(2)).^2+(z-xs(3)).^2);
     % scalar = (x-xs) nxs
@@ -102,7 +100,7 @@
     % G(x-xs,w) =  - -  H0  |  - |x-xs|  |
     %                4       \ c        /
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-S.ls
+    % https://sfs.rtfd.io/en/3.2/sources/#equation-fd-line
     %
     G = -1i/4 * besselh(0,2,omega/c* ...
         sqrt( (x-xs(1)).^2 + (y-xs(2)).^2 + (z-xs(3)).^2 ));
@@ -112,7 +110,7 @@
     %
     % G(x,w) = e^(-i w/c n x)
     %
-    % See: https://sfs.rtfd.io/en/3.2/#equation-S.pw
+    % https://sfs.rtfd.io/en/3.2/sources/#equation-fd-plane
     %
     % Direction of plane wave
     nxs = xs(:,1:3) / norm(xs(:,1:3));
diff --git a/SFS_monochromatic/sound_field_mono.m b/SFS_monochromatic/sound_field_mono.m
index c70b1a8f..f1c1fe9b 100644
--- a/SFS_monochromatic/sound_field_mono.m
+++ b/SFS_monochromatic/sound_field_mono.m
@@ -130,7 +130,7 @@
     % P(x,omega) = | D(x0,omega) G(x-x0,omega) dx0
     %              /
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-single-layer
+    % https://sfs.rtfd.io/en/3.2/problem/#equation-single-layer
     %
     % x0(ii,7) is a weight for the single secondary sources which includes for
     % example a tapering window for WFS or a weighting of the sources for
diff --git a/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_fs.m b/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_fs.m
index e93dfcdd..bcd52340 100644
--- a/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_fs.m
+++ b/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_fs.m
@@ -95,7 +95,7 @@
         % d(x0,t) = h(t) * --- ----------- delta(t+|x0-xs|/c)
         %                  2pi  |x0-xs|^2
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.fs
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-focused-3d
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -110,7 +110,7 @@
         % d(x0,t) = h(t) * _ |--- ------------- delta(t+|x0-xs|/c)
         %                   \|2pi |x0-xs|^(3/2)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.fs.ls
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-focused-2d
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -167,7 +167,7 @@
         % d_2.5D(x0,t) = h_pre(-t) * g0 _ |---  ------------- delta(t+|x0-xs|/c)
         %                                \|2pi  |x0-xs|^(3/2)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.fs.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-focused-25d
         %
         % Delay and amplitude weight
         delay = -1./c .* r;
diff --git a/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_ls.m b/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_ls.m
index 57f5ebf0..3b062a9a 100644
--- a/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_ls.m
+++ b/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_ls.m
@@ -88,7 +88,7 @@
         % d(x0,t) = h(t) * _ |--- ------------- delta(t-|x0-xs|/c)
         %                   \|2pi |x0-xs|^(3/2)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.ls
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-line
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -134,8 +134,7 @@
         % where v = x0-xs - <x0-xs,nxs > nxs,
         % and |nxs| = 1.
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.ls
-        % and https://sfs.rtfd.io/en/3.2/#equation-v.ls
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-line
         %
         % v = (I - nxs'nxs)(x0-xs)
         % r = |v|
diff --git a/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_ps.m b/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_ps.m
index 5a3fe3bd..49d1b014 100644
--- a/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_ps.m
+++ b/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_ps.m
@@ -85,7 +85,7 @@
         % d(x0,t) = h(t) * --- ----------- delta(t-|x0-xs|/c)
         %                  2pi  |x0-xs|^2
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.ps
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-point
         %
         % r = |x0-xs|
         r = vector_norm(x0-xs,2);
@@ -142,7 +142,7 @@
          % d_2.5D(x0,t) = h_pre(t) * g0 _ |---  ------------- delta(t-|x0-xs|/c)
          %                               \|2pi  |x0-xs|^(3/2)
          %
-         % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.ps.2.5D
+         % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-point-25d
          %
          % Delay and amplitude weight
          delay = 1./c .* r;
diff --git a/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_pw.m b/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_pw.m
index 9f4c5535..cda3f770 100644
--- a/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_pw.m
+++ b/SFS_time_domain/driving_functions_imp/driving_function_imp_wfs_pw.m
@@ -80,7 +80,7 @@
         %
         % d_2D(x0,t) = h(t) * 2 nk nx0 delta(t - 1/c nk x0)
         %
-        % See https://sfs.rtfd.io/en/3.2/#equation-d.wfs.pw
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-plane
         %
         % Delay and amplitude weight
         delay = 1./c .* vector_product(nk,x0,2);
@@ -111,7 +111,7 @@
         %
         % d_2.5D(x0,t) = h(t) * 2 g0 nk nx0 delta(t - 1/c nk x0)
         %
-        % See https://sfs.rtfd.io/en/3.2/en/update_wfs_ps/#equation-d.wfs.pw.2.5D
+        % https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-plane-25d
         %
         % Delay and amplitude weight
         delay = 1./c .* vector_product(nk,x0,2);
diff --git a/SFS_time_domain/greens_function_imp.m b/SFS_time_domain/greens_function_imp.m
index de2e5253..33ca4fd7 100644
--- a/SFS_time_domain/greens_function_imp.m
+++ b/SFS_time_domain/greens_function_imp.m
@@ -76,7 +76,7 @@
     % g(x-xs,t) = ---------- delta(t - |x-xs|/c)
     %             4pi |x-xs|
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-s.ps
+    % https://sfs.rtfd.io/en/3.2/sources/#equation-td-point
     %
     r = sqrt((x-xs(1)).^2+(y-xs(2)).^2+(z-xs(3)).^2);
     g = 1./(4*pi.*r);
@@ -89,8 +89,6 @@
     % g(x-xs,ns,t) = --- | F  | ---- | + ------  | --------- delta(t - |x-xs|/c)
     %                4pi  \    \ c  /    |x-xs| /  |x-xs|^2
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-s.dps
-    %
     to_be_implemented(mfilename);
 
 elseif strcmp('ls',src)
@@ -100,7 +98,7 @@
     % g(x-xs,t) = F |--  |  - |---  --_-_-_- delta(t - |x-xs|/c)
     %                \iw/    \|8pi  \||x-xs|
     %
-    % See https://sfs.rtfd.io/en/3.2/en/latest/#equation-s.ls
+    % https://sfs.rtfd.io/en/3.2/sources/#equation-td-line
     % Note, that the filter F^-1 is not implemented!!!!
     %
     r = sqrt((x-xs(1)).^2+(y-xs(2)).^2+(z-xs(3)).^2);
@@ -112,7 +110,7 @@
     %
     % g(x,t) = delta(t - nx/c)
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-s.pw
+    % https://sfs.rtfd.io/en/3.2/sources/#equation-td-plane
     %
     % direction of plane wave
     nxs = xs / norm(xs);
diff --git a/SFS_time_domain/sound_field_imp.m b/SFS_time_domain/sound_field_imp.m
index 4ad8ea5c..571fb556 100644
--- a/SFS_time_domain/sound_field_imp.m
+++ b/SFS_time_domain/sound_field_imp.m
@@ -171,7 +171,7 @@
     % p(x,t) = | d(x0,t) g(x-x0,t) dx0
     %          /
     %
-    % See https://sfs.rtfd.io/en/3.2/#equation-single-layer
+    % Compare https://sfs.rtfd.io/en/3.2/problem/#equation-fd-single-layer
     %
     % x0(ii,7) is a weight for the single secondary sources which includes for
     % example a tapering window for WFS or a weighting of the sources for
diff --git a/SFS_time_domain/wfs_fir_prefilter.m b/SFS_time_domain/wfs_fir_prefilter.m
index 70b78f4d..e7c71f35 100644
--- a/SFS_time_domain/wfs_fir_prefilter.m
+++ b/SFS_time_domain/wfs_fir_prefilter.m
@@ -93,7 +93,7 @@
     %  H(f) = \ |------ for flow<=f<=fhigh
     %          \|  c
     %
-    %  See https://sfs.rtfd.io/en/3.2/#equation-h.wfs.2.5D
+    %  https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-prefilter-25d
     %
     H(idxflow:idxfhigh) = sqrt(2*pi*f(idxflow:idxfhigh)/c);
     H(idxfhigh:end) = H(idxfhigh);
@@ -103,7 +103,7 @@
     %  H(f) = ------ for flow<=f<=fhigh
     %           c
     %
-    %  See https://sfs.rtfd.io/en/3.2/#equation-h.wfs
+    %  https://sfs.rtfd.io/en/3.2/d_wfs/#equation-td-wfs-prefilter
     %
     H(idxflow:idxfhigh) = 2*pi*f(idxflow:idxfhigh)/c;
     H(idxfhigh:end) = H(idxfhigh);