fix wrong dimension of gemm_op in rotate_wf

Author: Cstandardlib

source/module_hsolver/diago_bpcg.cpp

@@ -224,7 +224,7 @@ void DiagoBPCG<T, Device>::rotate_wf(
     gemm_op()(this->ctx,
                 'N',
                 'N',
-                this->n_dim,        //m
+                this->n_basis,        //m
                 this->n_band,       //n
                 this->n_band,       //k
                 this->one,          //1.0