From 183f8a20cb2c2c7115e49c3f81274aa0bf0ecd1b Mon Sep 17 00:00:00 2001
From: dzzz2001 <dzcpku2001@163.com>
Date: Wed, 5 Mar 2025 13:16:24 +0800
Subject: [PATCH] inline function of complexarray

---
 source/module_base/complexarray.cpp | 14 --------------
 source/module_base/complexarray.h   | 21 +++++++++++++++++++--
 2 files changed, 19 insertions(+), 16 deletions(-)

diff --git a/source/module_base/complexarray.cpp b/source/module_base/complexarray.cpp
index b7d2abab01..becb8e8448 100644
--- a/source/module_base/complexarray.cpp
+++ b/source/module_base/complexarray.cpp
@@ -262,18 +262,4 @@ void point_mult(ComplexArray &in1, ComplexArray &in2, ComplexArray &out){
 		              in1.ptr[i].real() * in2.ptr[i].imag() +
 		              in1.ptr[i].imag() * in2.ptr[i].real());}
 }
-const std::complex <double> &ComplexArray::operator()(const int ind1, const int ind2, const int ind3, const int ind4) const{
-	assert(ind1>=0);	assert(ind1<bound1);
-	assert(ind2>=0);	assert(ind2<bound2);
-	assert(ind3>=0);	assert(ind3<bound3);
-	assert(ind4>=0);	assert(ind4<bound4);
-	const int ind = ((ind1 * bound2 + ind2) * bound3 + ind3) * bound4 + ind4;
-	return ptr[ind];}
-std::complex<double>& ComplexArray::operator()(const int ind1,const int ind2,const int ind3,const int ind4){
-	assert(ind1>=0);	assert(ind1<bound1);
-	assert(ind2>=0);	assert(ind2<bound2);
-	assert(ind3>=0);	assert(ind3<bound3);
-	assert(ind4>=0);	assert(ind4<bound4);
-	const int ind = ((ind1 * bound2 + ind2) * bound3 + ind3) * bound4 + ind4;
-	return ptr[ind];}
 }
diff --git a/source/module_base/complexarray.h b/source/module_base/complexarray.h
index 8fbd162d2f..9cddc9a27c 100644
--- a/source/module_base/complexarray.h
+++ b/source/module_base/complexarray.h
@@ -5,6 +5,7 @@
 #include <iostream>
 #include <fstream>
 #include <iomanip>
+#include <cassert>
 
 namespace ModuleBase
 {
@@ -58,11 +59,27 @@ class ComplexArray
 
 	/// overloaded subscript operator for non-const std::complex Array const reference return creates an lvakue
 	std::complex <double> &operator()
-		(const int ind1=0, const int ind2=0, const int ind3=0, const int ind4=0);
+		(const int ind1=0, const int ind2=0, const int ind3=0, const int ind4=0)
+		{
+			assert(ind1>=0);	assert(ind1<bound1);
+			assert(ind2>=0);	assert(ind2<bound2);
+			assert(ind3>=0);	assert(ind3<bound3);
+			assert(ind4>=0);	assert(ind4<bound4);
+			const int ind = ((ind1 * bound2 + ind2) * bound3 + ind3) * bound4 + ind4;
+			return ptr[ind];
+		};
 	//  std::complex < double> &operator()(int, int, int, int, int);
 	/// overloaded subscript operator for const std::complex Array const reference return creates an cvakue
 	const std::complex <double> &operator()
-		(const int ind1=0, const int ind2=0, const int ind3=0, const int ind4=0)const;
+		(const int ind1=0, const int ind2=0, const int ind3=0, const int ind4=0) const
+		{
+			assert(ind1>=0);	assert(ind1<bound1);
+			assert(ind2>=0);	assert(ind2<bound2);
+			assert(ind3>=0);	assert(ind3<bound3);
+			assert(ind4>=0);	assert(ind4<bound4);
+			const int ind = ((ind1 * bound2 + ind2) * bound3 + ind3) * bound4 + ind4;
+			return ptr[ind];
+		};
 	//  const std::complex < double> &operator()(int, int, int, int, int)const;
 
 	/****************************************************