// clang-format off
inline void dsqmm1(double*a,double*b,double*c,int n){for(int i=0;i<n;++i){for(int j=0;j<n;++j){for(int k=0;k<n;++k){c[i+j*n]+=a[i+k*n]*b[k+j*n];}}}}