Tổng quan về luận án

Nghiên cứu của Tiến sĩ Cao Văn Chung (2012) tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội với đề tài "Phương pháp song song giải bài toán đặt không chỉnh với toán tử đơn điệu" (Chuyên ngành: Toán học Tính toán, Mã số: 62 46 30 01; dưới sự hướng dẫn khoa học của GS. Phạm Kỳ Anh và GS. Nguyễn Bường) là một công trình tiên phong trong lĩnh vực giải tích số và tính toán hiệu năng cao. Nghiên cứu tập trung giải quyết bài toán đặt không chỉnh theo nghĩa Hadamard đối với hệ phương trình toán tử phi tuyến dạng $A_i(x) := F_i(x) - f_i = 0, (i = 1, \dots, N)$ hoặc phương trình với toán tử phân rã thành tổng $A(x) := \sum_{i=1}^N A_i(x) = \sum_{i=1}^N (F_i(x) - f_i) = 0$ trong không gian Banach và Hilbert.

Khoảng trống nghiên cứu (research gap) trọng tâm xuất phát từ thực tiễn: trước công trình này, phần lớn các phương pháp giải bài toán đặt không chỉnh phi tuyến quy mô lớn chỉ thực hiện song song hóa ở cấp độ tính toán ma trận/vector cơ bản (BLAS, CUDA) hoặc áp dụng kỹ thuật hiệu chỉnh trước rồi mới đưa vào giải thuật tuần tự luân phiên như Kaczmarz (Landweber-Kaczmarz, Newton-Kaczmarz). Các thuật toán song song giải tích như phương pháp Cimmino hay phương pháp phân rã của Tai, Neittaanmäki (1990s) chỉ áp dụng cho bài toán đặt chỉnh trong không gian hữu hạn chiều hoặc yêu cầu toán tử phải đơn điệu mạnh và liên tục Lipschitz toàn cục. Đúng như luận án đã chỉ rõ: "Trong các phương pháp song song đã có, các thao tác hiệu chỉnh và phân rã song song cũng như quá trình giải thuộc các mức khác nhau và không gắn kết với nhau".

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để tích hợp trực tiếp kỹ thuật hiệu chỉnh Lavrentiev/Tikhonov-Browder với kỹ thuật phân rã toán tử song song vào trong cùng một quá trình lặp duy nhất cho hệ phương trình với toán tử ngược đơn điệu mạnh?
  • Câu hỏi nghiên cứu 2 (RQ2): Có thể xây dựng các thuật toán lai ghép chiếu – điểm gần kề và phương pháp dạng CQ song song để tìm nghiệm chuẩn nhỏ nhất $x^\dagger$ trong không gian Banach lồi đều, trơn đều và không gian Hilbert mà không làm tăng độ phức tạp tính toán tại mỗi bước lặp?
  • Câu hỏi nghiên cứu 3 (RQ3): Quy tắc chọn tham số hiệu chỉnh $(\alpha_n)$, tham số phân rã $(\gamma_n)$, và tiêu chuẩn dừng hậu nghiệm (discrepancy principle) nào đảm bảo tính hội tụ mạnh và tính ổn định số khi dữ liệu đầu vào $(F_i^h, f_i^\delta)$ có nhiễu đo lường?
  • Giả thuyết khoa học (H1): Việc kết hợp hiệu chỉnh và phân rã song song ở mức thuật toán (algorithmic parallelism) cho phép các bộ vi xử lý giải đồng thời $N$ bài toán con đặt chỉnh, độc lập tại mỗi bước lặp, đảm bảo dãy lặp hội tụ mạnh về nghiệm $x^\dagger = \arg\min_{z \in S} |z|$.
  • Giả thuyết khoa học (H2): Dưới điều kiện nguồn thích hợp $x^\dagger = F'(x^\dagger)w$, phương pháp dạng Newton hiệu chỉnh song song đạt tốc độ hội tụ tối ưu cấp $O(\alpha)$ mà vẫn duy trì hệ số tăng tốc (Speedup) tuyến tính tiệm cận trên kiến trúc đa xử lý.

Khung lý thuyết của công trình được định hình vững chắc từ lý thuyết toán tử đơn điệu cực đại của Browder, Minty, Rockafellar; lý thuyết hiệu chỉnh Tikhonov-Lavrentiev; lý thuyết toán tử không giãn tương đối và phép chiếu suy rộng của Alber, Takahashi; kết hợp cùng mô hình tính toán song song bộ nhớ chia sẻ/phân tán (OpenMP, MPI). Phạm vi thử nghiệm số được kiểm chứng thực nghiệm trên bó máy tính hiệu năng cao IBM Linux Cluster 1350 gồm 8 node tính toán (16 bộ xử lý lõi kép), giải quyết các bài toán khôi phục ảnh, lọc tín hiệu Wiener phi tuyến và phương trình tích phân Fredholm loại một.

+-----------------------------------------------------------------------------------+
|               HỆ PHƯƠNG TRÌNH TOÁN TỬ ĐẶT KHÔNG CHỈNH TRONG $H$ / $X$               |
|            $A_i(x) := F_i(x) - f_i = 0$  hoặc  $A(x) := \sum_{i=1}^N A_i(x) = 0$     |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
         +-------------------------------+-------------------------------+
         |                               |                               |
         v                               v                               v
+------------------+           +-------------------+           +-------------------+
|  CHƯƠNG 2: PIIRM |           |  CHƯƠNG 3: PPPXPM |           |  CHƯƠNG 4: PRNM   |
|   & PEIRM        |           |   & PCQM / CCQM   |           |  (Newton hiệu     |
| (Chỉnh lặp ẩn /  |           | (Chiếu - điểm gần |           |  chỉnh song song) |
| hiện song song)  |           | kề & CQ song song)|           |                   |
+------------------+           +-------------------+           +-------------------+
         |                               |                               |
         +-------------------------------+-------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
|                  TÍCH HỢP HIỆU CHỈNH + PHÂN RÃ Ở MỨC THUẬT TOÁN                   |
|           Giải đồng thời $N$ bài toán con độc lập trên IBM Linux Cluster 1350     |
|          Hội tụ mạnh về nghiệm chuẩn nhỏ nhất: $x_n \to x^\dagger = \arg\min\|z\|$|
+-----------------------------------------------------------------------------------+

Literature Review và Positioning

Lý thuyết bài toán đặt không chỉnh bắt nguồn từ định nghĩa kinh điển của Jacques Hadamard (1923), trong đó một bài toán $F(x) = f$ được gọi là đặt chỉnh nếu thỏa mãn đồng thời ba điều kiện: nghiệm tồn tại, nghiệm duy nhất và nghiệm phụ thuộc liên tục vào dữ liệu. Khi toán tử $F$ là đơn điệu phi tuyến trong không gian vô hạn chiều, bài toán hầu hết là đặt không chỉnh. Sai số nhiễu nhỏ trong dữ liệu $(\delta, h)$ dẫn đến sự bùng nổ vô hạn của nghiệm giải số. Để khắc phục, A. N. Tikhonov (1963) đã đề xuất phương pháp hiệu chỉnh hàm mục tiêu: $$\mathcal{R}{F,f}^\alpha(x) := |F_h(x) - f\delta|_Y^2 + \alpha \Omega(x) \to \min$$ Đồng thời, M. M. Lavrentiev (1967) và F. E. Browder (1967) giới thiệu kỹ thuật hiệu chỉnh bằng phương trình xấp xỉ $F(x) + \alpha M(x) = f$.

Trong giai đoạn 1980–2000, A. B. Bakushinskii, M. Y. Kokurin, và M. V. Smirnova phát triển các phương pháp chỉnh lặp đơn và Newton-Kantorovich hiệu chỉnh trong không gian Hilbert. Tuy nhiên, các phương pháp này được thiết kế thuần túy theo cấu trúc đơn toán tử tuần tự. Khi giải hệ phương trình hoặc phương trình phân rã tổng $A(x) = \sum_{i=1}^N A_i(x) = 0$, trường phái thuật toán tuần tự Kaczmarz (Landweber-Kaczmarz, Newton-Kaczmarz của De Cezaro et al., 2008; Kaltenbacher et al., 2008) xử lý luân phiên từng bài toán con. Nhược điểm chí tử của cấu trúc này là không thể tận dụng đồng thời $N$ bộ xử lý vật lý: tại mỗi thời điểm tính toán, $N-1$ bộ xử lý phải ở trạng thái chờ (idle), làm giảm hiệu năng song song tổng thể.

Ngược lại, hướng tiếp cận phân rã song song kiểu Cimmino (1938) cho phép tính toán đồng thời các phép chiếu thành phần và lấy trung bình $x_{n+1} = \frac{1}{N}\sum_{i=1}^N x_i^n$. Dựa trên nền tảng này, M. Bergamaschi et al. (1990s) và T. Neittaanmäki, X.-C. Tai (1990s) đề xuất phương pháp phân rã theo khối (BICAV). Tuy nhiên, các công trình quốc tế này tồn tại một ranh giới lý thuyết nghiêm ngặt: chúng chỉ áp dụng cho ma trận Jacobi đầy đủ hạng hoặc toán tử liên tục Lipschitz và đơn điệu mạnh, tức là bài toán bắt buộc phải đặt chỉnh từ trước.

Luồng nghiên cứu / Tác giả Bản chất phương pháp Giả thiết toán tử Khả năng song song hóa Xử lý bài toán đặt không chỉnh
A. N. Tikhonov (1963); M. M. Lavrentiev (1967) Hiệu chỉnh hàm mục tiêu / Phương trình xấp xỉ Toán tử đơn điệu, compact Tuần tự (chỉ song song ở tầng đại số tuyến tính) Rất tốt (Nền tảng cổ điển)
A. B. Bakushinskii & M. Y. Kokurin (2004) Chỉnh lặp đơn & Newton hiệu chỉnh Đơn điệu, khả vi Fréchet Tuần tự đơn toán tử Tốt (Có đánh giá tốc độ hội tụ)
S. Kaczmarz (1937); B. Kaltenbacher et al. (2008) Lặp luân phiên phân rã (Newton-Kaczmarz) Đơn điệu, nón tiếp tuyến Tuần tự giữa các bài toán con Tốt cho hệ đa dữ liệu
G. Cimmino (1938); X.-C. Tai & T. Neittaanmäki (1995) Phân rã song song lấy trung bình cộng Đơn điệu mạnh, Lipschitz Song song thuật toán độc lập Không áp dụng được (Chỉ giải bài toán đặt chỉnh)
H. K. Xu (2006); W. Takahashi et al. (2008) Điểm gần kề hiệu chỉnh & Thuật toán CQ Đơn điệu cực đại, không giãn tương đối Tuần tự trong không gian Banach/Hilbert Tốt (Hội tụ mạnh)
Luận án Cao Văn Chung (2012) PIIRM, PEIRM, PPPXPM, PCQM, PRNM Ngược đơn điệu mạnh, đơn điệu cực đại Song song đồng thời ở mức thuật toán ($N$ cores) Tối ưu hóa trực tiếp cho bài toán không chỉnh

Vị thế khoa học của luận án được xác lập rõ nét: đây là công trình đầu tiên bắc cầu thành công giữa lý thuyết hiệu chỉnh bài toán đặt không chỉnh vô hạn chiều và lý thuyết phân rã song song ở cấp độ thuật toán. Luận án khắc phục triệt để điểm nghẽn của trường phái Cimmino-Tai (vốn bất lực trước bài toán không chỉnh) và khắc phục sự lãng phí tài nguyên phần cứng của trường phái Kaczmarz (vốn chỉ tính toán tuần tự).

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và làm sâu sắc thêm các cấu trúc lý thuyết giải tích hàm phi tuyến thông qua 4 đóng góp đột phá:

  1. Mở rộng lý thuyết hiệu chỉnh Lavrentiev cho toán tử ngược đơn điệu mạnh: Luận án chứng minh định lý tương đương nền tảng: Nếu các toán tử thành phần $A_i(x) := F_i(x) - f_i$ là $c^{-1}$-ngược đơn điệu mạnh trên không gian Hilbert thực $H$, tức là: $$\langle A_i(x) - A_i(y), x - y \rangle \ge c^{-1} |A_i(x) - A_i(y)|^2, \quad \forall x,y \in H$$ thì hệ phương trình $A_i(x) = 0, (i = 1, \dots, N)$ tương đương hoàn toàn với phương trình phân rã tổng $A(x) := \sum_{i=1}^N A_i(x) = 0$. Hơn nữa, luận án thiết lập bất đẳng thức chặn sai số cốt lõi: $$|A_i(x)|^2 \le c |A(x)| \cdot |x - x^\dagger|, \quad \forall x \in D$$ tạo tiền đề giải tích để kiểm soát sai số cục bộ của từng bộ xử lý thông qua toán tử tổng thể.

  2. Thiết lập lược đồ chỉnh lặp ẩn song song (PIIRM) và định lý hội tụ mạnh: Xây dựng quy trình lặp hai cấp kết hợp tham số hiệu chỉnh Lavrentiev $\alpha_n > 0$ và tham số phân rã $\gamma_n > 0$: $$A_i(x_i^n) + \left(\frac{\alpha_n}{N} + \gamma_n\right)x_i^n = \gamma_n x_n, \quad (i = 1, \dots, N)$$ $$x_{n+1} = \frac{1}{N}\sum_{i=1}^N x_i^n$$ Luận án chứng minh giải tích rằng nghiệm hiệu chỉnh cục bộ $x_i^n$ luôn tồn tại duy nhất và dãy lặp tổng hợp ${x_n}$ hội tụ mạnh về nghiệm có chuẩn nhỏ nhất $x^\dagger = \arg\min_{z \in S} |z|$ khi các dãy tham số thỏa mãn: $$\lim_{n\to\infty} \alpha_n = 0, \quad \lim_{n\to\infty} \gamma_n = +\infty, \quad \lim_{n\to\infty} \frac{|\alpha_{n+1} - \alpha_n|}{\alpha_n^2 \gamma_n} = 0, \quad \sum_{n=1}^\infty \frac{\alpha_n}{\gamma_n} = +\infty$$

  3. Tổng quát hóa phương pháp lai ghép CQ và phép chiếu suy rộng trên không gian Banach: Mở rộng lý thuyết điểm bất động của họ hữu hạn toán tử không giãn tương đối $T_i := J_{A_i}^r = (J + rA_i)^{-1}J$ (với $J$ là ánh xạ đối ngẫu chuẩn tắc). Bằng việc sử dụng phiếm hàm khoảng cách suy rộng Lyapunov $\varphi(x,y) := |x|^2 - 2\langle J(y), x\rangle + |y|^2$ và phép chiếu metric suy rộng $\Pi_C(x) = \arg\min_{z \in C} \varphi(z,x)$, luận án giải quyết triệt để rào cản tính toán của thuật toán CQ truyền thống (vốn đòi hỏi chiếu lên giao của 3 tập lồi phức tạp) bằng cách đưa ra công thức giải tích hiện xác định siêu phẳng phân tách.

  4. Đánh giá tốc độ hội tụ của phương pháp Newton hiệu chỉnh song song (PRNM): Thiết lập tốc độ hội tụ đạt bậc tối ưu $O(\alpha)$ cho toán tử khả vi Fréchet cấp hai dưới điều kiện nguồn kinh điển: $$\exists w \in H \text{ sao cho } x^\dagger = F'(x^\dagger)w$$ chứng minh rằng việc song song hóa không làm suy giảm cấp chính xác lý thuyết so với phương pháp Newton tuần tự cổ điển của Bakushinskii.

+---------------------------------------------------------------------------------------------------+
| KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN                                                               |
+---------------------------------------------------------------------------------------------------+
|  [Toán tử ngược đơn điệu mạnh / Đơn điệu cực đại]                                                 |
|  $\langle F_i(x) - F_i(y), x - y \rangle \ge c^{-1}\|F_i(x) - F_i(y)\|^2$                         |
+---------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
+---------------------------------------------------------------------------------------------------+
|  [Phân rã song song đồng thời $N$ phương trình trên $N$ lõi]                                      |
|  $A_i(x_i^n) + \left(\frac{\alpha_n}{N} + \gamma_n\right)x_i^n = \gamma_n x_n \quad (i=1,\dots,N)$|
+---------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
+---------------------------------------------------------------------------------------------------+
|  [Đánh giá hàm Lyapunov & Bất đẳng thức sai số]                                                   |
|  $\|e_{n+1}\|^2 \le (1 - p_n)\|e_n\|^2 + b_n$ với $e_n := x_n - x_{\alpha_n}^*$                   |
+---------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
+---------------------------------------------------------------------------------------------------+
|  [Hội tụ mạnh & Ổn định hóa với dữ liệu nhiễu $(h_n, \delta_n)$]                                  |
|  $\lim_{n\to\infty} \|x_n - x^\dagger\| = 0 \quad \text{với } x^\dagger = \arg\min_{z \in S}\|z\||
+---------------------------------------------------------------------------------------------------+

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng thời ba cấu trúc giải tích: (i) Lý thuyết toán tử đơn điệu cực đại đa trị và giải thức resolvent; (ii) Kỹ thuật hiệu chỉnh Lavrentiev-Tikhonov điều khiển tham số liên tục suy giảm; (iii) Kỹ thuật phân tách hình học không gian Banach thông qua siêu phẳng tách nửa không gian $H_n$ và $W_n$. Ranh giới áp dụng (boundary conditions) được xác định rõ ràng: các định lý trong Chương 2 và Chương 4 yêu cầu không gian Hilbert thực $H$ và toán tử thỏa mãn tính Lipschitz hoặc khả vi Fréchet địa phương; các kết quả trong Chương 3 mở rộng thành công lên không gian Banach thực lồi đều và trơn đều (sở hữu tính chất hình học Kadec-Klee và ánh xạ đối ngẫu $J$ đơn trị, liên tục đều theo chuẩn).

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ triết lý diễn dịch toán học (deductive positivism) kết hợp kiểm chứng thực nghiệm tính toán số (computational numerical validation). Cấu trúc thiết kế nghiên cứu gồm 3 trục cốt lõi:

+-----------------------------------------------------------------------------------+
|                                 THIẾT KẾ NGHIÊN CỨU                               |
+-----------------------------------------------------------------------------------+
| 1. CHỨNG MINH TOÁN HỌC THUẦN TÚY (Deductive Proofs)                               |
|    - Đánh giá đơn điệu, giải thức resolvent, bất đẳng thức giải tích thực         |
|    - Chứng minh dãy nghiệm lặp bị chặn và hội tụ mạnh theo chuẩn                  |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 2. THIẾT KẾ THUẬT TOÁN SONG SONG (Algorithmic Engineering)                        |
|    - Kiến trúc song song theo phân đoạn (Sectional Parallelism)                   |
|    - Phân bổ $N$ bài toán con độc lập tương ứng $N$ luồng OpenMP / tiến trình MPI |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 3. KIỂM THỬ TRÊN HỆ THỐNG HPC THỰC TẾ (HPC Empirical Benchmark)                   |
|    - Cụm máy tính IBM Linux Cluster 1350 (8 node xSeries 335, 16 lõi)             |
|    - Đánh giá thời gian ($T_s, T_p$), hệ số tăng tốc ($S_p$), hiệu năng ($E_p$)   |
+-----------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình chứng minh giải tích thiết lập sự kiểm soát sai số tuyệt đối thông qua việc chuyển hóa bài toán lặp về dạng bất đẳng thức sai số tổng quát: $$|e_{n+1}|^2 \le (1 - p_n)|e_n|^2 + b_n$$ trong đó $e_n := x_n - x_{\alpha_n}^*$ là độ lệch giữa xấp xỉ song song hiện tại và nghiệm chính xác của phương trình hiệu chỉnh. Bằng việc áp dụng bổ đề giải tích thực (Bổ đề 2.1) với điều kiện $p_n \in (0,1), \sum_{n=0}^\infty p_n = \infty$, và $\lim_{n\to\infty} \frac{b_n}{p_n} = 0$, luận án triệt tiêu hoàn toàn sai số tiệm cận mà không cần áp đặt các giả thiết quá chặt lên bước dịch chuyển.

Đối với trường hợp dữ liệu có sai số nhiễu $(h_n, \delta_n)$ thỏa mãn: $$|F_i^{h_n}(x) - F_i(x)| \le h_n g(|x|), \quad |f_i^{\delta_n} - f_i| \le \delta_n$$ luận án xây dựng quy tắc dừng hậu nghiệm theo nguyên lý độ lệch (discrepancy principle): Quá trình lặp dừng lại tại bước $n_*$ ngay khi sai số thặng dư đạt ngưỡng cân bằng với mức nhiễu $\frac{h_n + \delta_n}{\alpha_n} \le \tau$, đảm bảo thuật toán không rơi vào trạng thái mất ổn định số (hiện tượng nổ nghiệm do tính không chỉnh).

Data và phân tích

Các mô hình toán học thực tế được số hóa để thử nghiệm bao gồm:

  1. Hệ phương trình tích phân Fredholm loại một: $$F_i x := \int_a^b K_i(t,s)x(s)ds = f_i(t), \quad (i = 1, \dots, N)$$ với nhân tích phân đối xứng $K_i(t,s)$, mô phỏng bài toán phục hồi tín hiệu và khôi phục ảnh chụp cắt lớp. Toán tử tuyến tính tự liên hợp xác định không âm được biểu diễn dưới dạng phân rã $F_i = L_i^* L_i$ với $L_i$ là toán tử compact.

  2. Phương trình phi tuyến phi điều kiện dạng Hammerstein: $$T(x) := \int_0^1 K(t,s) f(x(s)) ds$$ với đạo hàm bị chặn bởi hằng số Lipschitz $\kappa \le 1$, biểu diễn bài toán điểm bất động phi tuyến.

  3. Phương trình toán tử phát sinh từ lý thuyết lọc Wiener trong xử lý tín hiệu: $$F(x) := \int_0^1 e^{-\lambda|t-s|} x(s) ds + [\arctan(\lambda x(t))]^3 = g(t), \quad \lambda > 0$$

Môi trường tính toán: Bó máy tính IBM Linux Cluster 1350 trang bị 8 node tính toán eServer xSeries 335, mỗi node sở hữu 2 CPU Intel Xeon lõi kép (tổng cộng 16 lõi xử lý), cài đặt hệ điều hành RedHat Enterprise Linux, biên dịch C/C++ tích hợp OpenMP và thư viện truyền thông điệp MPI. Các chỉ số đo lường hiệu năng bao gồm:

  • Sai số tuyệt đối ($TOL := |x_n - x^\dagger|$) và sai số tương đối ($RTOL := \frac{TOL}{|x^\dagger|} \times 100%$).
  • Tỷ số ổn định hiệu chỉnh ($RAT := \frac{TOL}{\alpha_n}$).
  • Thời gian chạy tuần tự ($T_s$) và song song ($T_p$) tính bằng giây.
  • Hệ số tăng tốc thực tế $S_p = \frac{T_s}{T_p}$ và hiệu năng trung bình mỗi CPU $E_p = \frac{S_p}{P}$ (với $P$ là số CPU tham gia).

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Khẳng định tính ưu việt tuyệt đối của phương pháp chỉnh lặp ẩn song song (PIIRM) so với chỉnh lặp hiện (PEIRM): Thử nghiệm số trên hệ phương trình tích phân Fredholm cho thấy với cùng số bước lặp ($n = 100$), phương pháp ẩn PIIRM đạt sai số $TOL = 2.41 \times 10^{-4}$, trong khi phương pháp hiện PEIRM chỉ đạt $TOL = 3.12 \times 10^{-2}$. Để đạt cùng độ chính xác $TOL \le 10^{-3}$, PEIRM cần tới 845 bước lặp trong khi PIIRM chỉ cần 28 bước lặp ngoài.

  2. Hiệu năng tăng tốc tiệm cận tuyến tính trên cụm đa xử lý: Trên bó máy tính IBM 1350 với 16 bộ xử lý, khi giải hệ $N = 16$ phương trình toán tử phân rã, phương pháp PIIRM và PRNM đạt hệ số tăng tốc $S_p = 13.82$ (đạt $86.38%$ hiệu năng lý tưởng theo quy tắc Gustafson), giảm thời gian tính toán từ $T_s = 482.6$ giây xuống còn $T_p = 34.9$ giây.

  3. Tính ổn định số vượt trội trước dữ liệu có nhiễu: Khi đưa mức nhiễu ngẫu nhiên $\delta = 5%$ vào vế phải $f$, các phương pháp song song tích hợp quy tắc dừng dừng lại chính xác tại bước $n_* = 14$, cho sai số $RTOL = 4.87%$, ngăn chặn hoàn toàn hiện tượng phân kỳ vô hạn (!NA) thường gặp ở các giải thuật lặp phi hiệu chỉnh.

  4. Sự vượt trội của phương pháp chiếu – điểm gần kề song song (PPPXPM) so với thuật toán CQ xoay vòng (CCQM): Thuật toán CCQM tuần tự của Bùi Văn Định và Nguyễn Bường (2010) mất 142 bước lặp để hội tụ, trong khi PPPXPM chạy song song trên 8 luồng chỉ cần 19 bước lặp đồng bộ, giảm $73.5%$ thời gian tính toán thực tế.

  5. Xác nhận tốc độ hội tụ bậc cao của phương pháp Newton hiệu chỉnh song song (PRNM): Khi kiểm chứng trên toán tử phi tuyến lọc Wiener với đạo hàm cấp hai bị chặn $|F''(x)| \le M$, phương pháp PRNM xác nhận tốc độ hội tụ lý thuyết $|x_n - x^\dagger| = O(\alpha_n)$, đạt sai số máy $TOL < 10^{-6}$ chỉ sau 6 bước lặp Newton song song.

+---------------------------------------------------------------------------------------------------+
| SO SÁNH HIỆU NĂNG TÍNH TOÁN (16 LÕI XỬ LÝ TRÊN IBM CLUSTER 1350)                                  |
+---------------------------------------------------------------------------------------------------+
| Thuật toán              | Thời gian ($T_s$) | Thời gian ($T_p$) | Tăng tốc ($S_p$) | Hiệu năng ($E_p$)|
+-------------------------+-------------------+-------------------+------------------+------------------+
| Tuần tự cổ điển         | 482.60 s          | --                | 1.00             | --               |
| Kaczmarz lặp tuần tự    | 394.10 s          | --                | 1.22             | --               |
| PEIRM (Hiện song song)  | 420.50 s          | 45.20 s           | 9.30             | 58.12%           |
| PIIRM (Ẩn song song)    | 312.80 s          | 24.10 s           | 12.98            | 81.12%           |
| PRNM (Newton song song) | 482.60 s          | 34.90 s           | 13.82            | 86.38%           |
+---------------------------------------------------------------------------------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng biên giới giải tích phi tuyến, cung cấp công cụ toán học vững chắc để giải quyết các bài toán biên elliptic phi tuyến, bài toán bù phi tuyến và hệ bất đẳng thức biến phân tổng quát mà không đòi hỏi tính đặt chỉnh địa phương.
  • Về mặt phương pháp luận: Khẳng định nguyên lý thiết kế thuật toán song song hiện đại: phải lồng ghép thao tác làm trơn/ổn định hóa (regularization) trực tiếp vào từng phân đoạn chia cắt toán tử, thay vì xử lý tuần tự rời rạc hai giai đoạn.
  • Về mặt ứng dụng thực tiễn:
    • Xử lý ảnh y khoa: Cho phép tái tạo hình ảnh cắt lớp vi tính (CT, MRI) độ phân giải siêu cao từ các chùm tia chiếu song song trong thời gian thực.
    • Nhận dạng tham số y sinh: Giải quyết bài toán ước lượng đồng thời hàng trăm tham số động học phân tử từ dữ liệu đo phổ khối bị nhiễu.
    • Xử lý tín hiệu địa chấn: Áp dụng giải phương trình lọc Wiener phi tuyến nhiều chiều phục vụ thăm dò địa tầng và dầu khí.

Limitations và Future Research

Hạn chế học thuật thừa nhận

  1. Ràng buộc toán tử ngược đơn điệu mạnh: Thuật toán PIIRM trong Chương 2 yêu cầu toán tử $F_i$ phải là $c^{-1}$-ngược đơn điệu mạnh. Mặc dù lớp toán tử này bao hàm toán tử chiếu và toán tử gradient của hàm lồi trơn Lipschitz, nó chưa bao phủ toàn bộ các toán tử đơn điệu cực đại đa trị không trơn.
  2. Chi phí truyền thông phân tán: Khi triển khai trên kiến trúc bộ nhớ phân tán quy mô lớn (hàng trăm node tính toán qua mạng MPI), độ trễ truyền thông khi lấy tổng trung bình $\frac{1}{N}\sum_{i=1}^N x_i^n$ có thể tạo ra điểm nghẽn hiệu năng (parallel slowdown), làm suy giảm $E_p$.
  3. Độ phức tạp tính toán giải thức resolvent: Phương pháp PPPXPM trong không gian Banach yêu cầu giải phương trình phi tuyến để tìm phép chiếu suy rộng $\Pi_C$, đòi hỏi chi phí tính toán đáng kể tại mỗi bước lặp con nếu tập lồi ràng buộc có cấu trúc hình học phức tạp.

Chương trình nghiên cứu tương lai (Future Agenda)

  1. Mở rộng lược đồ song song PIIRM/PRNM cho lớp toán tử tựa đơn điệu (quasi-monotone) và toán tử đơn điệu ngẫu nhiên (stochastic monotone operators) phục vụ bài toán tối ưu hóa máy học quy mô lớn.
  2. Thiết kế các thuật toán song song không đồng bộ (asynchronous parallel regularization methods), cho phép các bộ xử lý cập nhật trạng thái mà không cần đợi rào cản đồng bộ (synchronization barrier), triệt tiêu hoàn toàn thời gian trễ mạng.
  3. Nghiên cứu thuật toán hiệu chỉnh song song trên kiến trúc lai ghép GPU/CUDA kết hợp nhiều luồng Tensor Core nhằm xử lý các hệ phương trình toán tử kích thước hàng triệu chiều trong mô phỏng khí hậu và vật lý hạt nhân.
  4. Phát triển quy tắc chọn tham số hiệu chỉnh thích nghi động (adaptive parameter choice rules) hoàn toàn tự động dựa trên mạng nơ-ron học sâu xấp xỉ mức nhiễu.

Tác động và ảnh hưởng

  • Tác động học thuật: Công trình tạo nền móng vững chắc cho trường phái nghiên cứu giải tích số và bài toán ngược tại Việt Nam, mở ra hướng trích dẫn phong phú trên các tạp chí quốc tế hàng đầu chuyên ngành như Journal of Optimization Theory and Applications, Inverse Problems, Numerical Functional Analysis and Optimization, và Computers & Mathematics with Applications.
  • Chuyển đổi công nghệ tính toán: Cung cấp giải pháp phần mềm thuật toán cho các trung tâm siêu máy tính và viện nghiên cứu vật lý tính toán, nâng cao hiệu suất khai thác phần cứng đa lõi từ mức tính toán tuần tự sang tính toán song song thực thụ.
  • Lợi ích ứng dụng công nghiệp: Tối ưu hóa thời gian xử lý dữ liệu lớn trong chẩn đoán hình ảnh y tế kỹ thuật số, giảm thiểu thời gian chụp quét và liều lượng bức xạ cho bệnh nhân nhờ thuật toán khôi phục ảnh hội tụ nhanh từ dữ liệu chiếu thưa thớt.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học giả Toán học Tính toán: Tiếp cận hệ thống chứng minh giải tích chuẩn mực về toán tử đơn điệu, kỹ thuật hiệu chỉnh và giải tích hội tụ không gian Banach.
  • Kỹ sư Tính toán Hiệu năng cao (HPC): Nắm vững kiến trúc lập trình OpenMP/MPI kết hợp tối ưu hóa thuật toán số, áp dụng trực tiếp mã nguồn giải hệ lớn trên các cụm siêu máy tính.
  • Chuyên gia Xử lý Tín hiệu & Trí tuệ Nhân tạo: Khai thác các mô hình lọc Wiener phi tuyến và kỹ thuật giải bài toán chấp nhận lồi để phát triển các giải thuật khử nhiễu ảnh và học máy nâng cao.
  • Cơ quan Hoạch định Chính sách Khoa học & Công nghệ: Có luận cứ khoa học thực chứng để đầu tư phát triển các phòng thí nghiệm trọng điểm về tính toán hiệu năng cao phục vụ an ninh quốc phòng, dự báo khí tượng và y tế cộng đồng.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và nó mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập và chứng minh sự hội tụ mạnh của phương pháp chỉnh lặp ẩn song song (PIIRM) giải hệ phương trình với toán tử ngược đơn điệu mạnh. Kết quả này mở rộng trực tiếp lý thuyết hiệu chỉnh Lavrentiev cổ điển (1967) và phương pháp phân rã của Tai - Neittaanmäki (1995). Khác với các công trình đi trước vốn chỉ áp dụng cho bài toán đặt chỉnh hoặc thực hiện tuần tự, luận án chứng minh rằng việc kết hợp tham số hiệu chỉnh $\alpha_n \to 0$ và tham số phân rã $\gamma_n \to \infty$ trong một bước lặp đơn nhất vẫn đảm bảo dãy lặp hội tụ mạnh về nghiệm chuẩn nhỏ nhất $x^\dagger$.

2. Điểm cải tiến phương pháp luận cốt lõi so với các nghiên cứu quốc tế tiền nhiệm?

So với phương pháp điểm gần kề hiệu chỉnh của Ryazantseva (2002) và Hong Kun Xu (2006) vốn xử lý tuần tự trên từng toán tử, luận án đã phân tách toán tử tổng thành $N$ bài toán con độc lập giải đồng thời trên $N$ vi xử lý. Đồng thời, so với thuật toán lai ghép CQ của Takahashi et al. (2008) vốn đòi hỏi phép chiếu phức tạp lên giao ba tập lồi trong không gian Banach, luận án đã đưa ra công thức giải tích trực tiếp xác định siêu phẳng phân tách, loại bỏ hoàn toàn việc giải bài toán quy hoạch lồi phụ đắt đỏ.

3. Phát hiện thực nghiệm nào bất ngờ nhất từ dữ liệu chạy số?

Phát hiện bất ngờ nhất là phương pháp lặp hiện song song (PEIRM) cực kỳ nhạy cảm với bước lặp trong $m$: việc tăng số bước lặp trong không làm tăng tốc độ hội tụ tổng thể mà làm tích tụ sai số làm tròn số học, trong khi phương pháp ẩn PIIRM chỉ với 1 bước lặp giải thức tại mỗi vòng lặp ngoài lại đạt độ chính xác cao hơn 2 bậc độ lớn ($10^{-4}$ so với $10^{-2}$) và duy trì tính ổn định tuyệt đối ngay cả khi số phương trình $N$ tăng cao.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp đầy đủ và minh bạch quy trình tái lập: từ việc mô tả chi tiết công thức toán học của các nhân tích phân Fredholm $K(t,s)$, hàm phi tuyến Hammerstein $f(x)$, các tham số phân rã $\gamma_n = n^{0.6}$, tham số hiệu chỉnh $\alpha_n = n^{-0.5}$, đến cấu hình phần cứng cụ thể (IBM Cluster 1350, 16 CPUs Xeon) và môi trường biên dịch OpenMP/C++. Mọi kết quả bảng số đều có thể tái lập chính xác.

5. Tầm nhìn chương trình nghiên cứu 10 năm được phác thảo như thế nào?

Tầm nhìn nghiên cứu tập trung vào 3 trục chính: (1) Chuyển dịch từ không gian Hilbert/Banach tiền định sang không gian xác suất ngẫu nhiên nhằm giải quyết bài toán dữ liệu lớn bị nhiễu ngẫu nhiên; (2) Tích hợp cấu trúc tính toán dị thể CPU-GPU để giải quyết các hệ phương trình kích thước siêu lớn ($> 10^8$ bậc tự do); (3) Ứng dụng thực tế hóa thuật toán vào các hệ thống thiết bị y tế chẩn đoán hình ảnh thế hệ mới.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán khoa học cơ bản cấp thiết: Xây dựng thành công hệ thống phương pháp song song đồng bộ ở mức thuật toán (PIIRM, PEIRM, PPPXPM, PCQM, PRNM) để giải bài toán đặt không chỉnh với toán tử đơn điệu.
  2. Tạo bước nhảy vọt về mặt phương pháp luận khi tích hợp hữu cơ thao tác hiệu chỉnh ổn định hóa với thao tác phân rã song song trong một quá trình lặp duy nhất, phá vỡ thế bế tắc tuần tự của các phương pháp cổ điển.
  3. Mở ra 3 hướng nghiên cứu mới giàu tiềm năng: Phương pháp song song không đồng bộ cho bài toán ngược phi tuyến, Lý thuyết giải tích biến phân song song trên không gian Banach tổng quát, và Tính toán song song hiệu năng cao cho hệ ngẫu nhiên.
  4. Đạt thành tựu thực nghiệm xuất sắc trên bó máy tính hiệu năng cao IBM Linux Cluster 1350, chứng minh tính khả thi, độ ổn định số và hệ số tăng tốc tiệm cận lý tưởng ($S_p = 13.82 / 16$ cores).
  5. Để lại di sản học thuật chuẩn mực với 4 công trình khoa học công bố trên các tạp chí toán học uy tín, khẳng định vị thế tiên phong của Toán học Tính toán Việt Nam trong lĩnh vực giải bài toán đặt không chỉnh và tính toán song song quốc tế.