Avoid redundant PPCG Gram initialization

Silver-Moon-Over-Snow · Silver-Moon-Over-Snow · commit d33ddb055259 · 2026-07-04T10:11:36.000+08:00
diff --git a/source/source_hsolver/ppcg/diago_ppcg_ops.hpp b/source/source_hsolver/ppcg/diago_ppcg_ops.hpp
@@ -119,7 +119,7 @@ void DiagoPPCG<T, Device>::gram(const T* a, const T* b,
                                  std::vector<T>& out,
                                  int ld_out) const
 {
-    out.assign(ld_out * ncol_b, T(0));
+    out.resize(ld_out * ncol_b);
     const T one = T(1);
     const T zero = T(0);
     ModuleBase::gemm_op<T, Device>()('C',

Original file line number	Diff line number	Diff line change
`@@ -119,7 +119,7 @@ void DiagoPPCG<T, Device>::gram(const T* a, const T* b,`
`119`	`119`	`std::vector<T>& out,`
`120`	`120`	`int ld_out) const`
`121`	`121`	`{`
`122`		`- out.assign(ld_out * ncol_b, T(0));`
	`122`	`+ out.resize(ld_out * ncol_b);`
`123`	`123`	`const T one = T(1);`
`124`	`124`	`const T zero = T(0);`
`125`	`125`	`ModuleBase::gemm_op<T, Device>()('C',`