双11狂欢大促

盲盒赢 定制羽绒服

开启惊喜盲盒，抽取你的专属好礼

— 惊喜奖品抢先看 —

超级年卡50元券

周边70元券

手机支架

一拖三数据线

雨伞

运动毛巾

电动牙刷

VIP月卡

机械键盘

定制羽绒服

超级年卡50元券

周边70元券

手机支架

一拖三数据线

雨伞

运动毛巾

电动牙刷

VIP月卡

机械键盘

定制羽绒服

优化矩阵乘法程序求bug解决方法

2014-01-09 01:44:17 226 阅读 AI 写同款 · GEO 优化AI 选题数据监测一键分发 ›

chenlei1700 墨衍会员 · AI 创作全网分发关注 + 加入社区

C++实现Strassen矩阵乘法：从理论到工程优化的完整指南矩阵乘法是计算机科学和工程计算中的基础操作，其标准实现的时间复杂度为O(n³)，在处理大规模数据时成为性能瓶颈。Strassen算法通过巧妙的数学变换，将8次子矩阵乘法减少为7次，将理论复杂度降至O(n^2.807)，实现了从概念到原理的突破。该算法的技术价值在于为大规模稠密矩阵运算提供了加速可能，尤其在机器学习、计算机图形学和高性能计算等应用场景中潜力巨大。然而，其工程实现面临递归阈值选择、内存布局优化和缓存友好性等核心挑战。本文通过C++实践，深入探讨了如何设计高效的Matrix类、管理递归过程中的临时阅读详情

我照着资料伪代码写了个矩阵优化算法思路就是充分利用1级缓存资源把大矩阵分成若干个小矩阵进行乘算
问题是我算10*10的矩阵没有问题算1000*1000就出bug了
我估计是内存分配或者指针越界的问题但是我找不出来了
求高手帮我找找到底哪里出错了。
代码如下（ubuntu下编译的）





#include<sys/types.h>

#include <stdio.h>

#include <math.h>

#include <stdlib.h>

#include <string.h>

#include <unistd.h>

#include <sched.h>

#include <time.h>

#include <sys/time.h>

#include <pthread.h>

#include <sys/wait.h>

#define N 100

#define THREAD_NUM 1



void thread_func(void *arg){

	printf("thread start \n"); 

	int i,j,k,i2,j2,k2;



	double **a;   

	double **b;

	double **c;

	double *aa,*bb,*cc;        

        

 

    a = (double**)malloc(sizeof(int*)*N);   

    for (i=0; i<N; i++)   

        a[i] = (double*)malloc(sizeof(double)*N);



    b = (double**)malloc(sizeof(double*)*N);   

    for (i=0; i<N; i++)   

        b[i] = (double*)malloc(sizeof(double)*N);

    

    c = (double**)malloc(sizeof(double*)*N);   

    for (i=0; i<N; i++)   

        c[i] = (double*)malloc(sizeof(double)*N);

	





	srand((unsigned)time(NULL));

	for(i = 0; i < N; i++)

		for(j = 0; j <N; j++)

			*(*(a+i)+j) = 555.123;

	for(i = 0; i < N; i++)

			for(j = 0; j <N; j++){

				*(*(b+i)+j) = 666.456;

			}

	for(i = 0; i < N; i++)

		for(j = 0; j <N; j++){

			*(*(c+i)+j) = 0.0;

			if(j==N-1&&i==N-1)

					printf("*(*(b+i)+j) =%f \n",*(*(b+i)+j));

		}



	for(i = 0; i < N; i+=8)

	for(j = 0; j < N; j+=8)	

	for(k = 0; k < N; k+=8)

		for(i2=0,cc=&(*(*(c+i)+j)),aa=&(*(*(a+i)+k));i2<8;++i2,cc+=N,aa+=N)

		for(k2=0,bb=&(*(*(b+k)+j));k2<8;++k2,bb+=N)

		for(j2=0;j2<8;++j2){

			*(cc+j2)+=*(aa+k2)*(*(bb+j2));	

			//printf("cc ,*(cc+j2)%d %d\n",cc,*(cc+j2));

			

		}





 printf("ready to free \n"); 



for (i=0; i<N; i++)   

        free(a[i]); free(b[i]);free(c[i]);  

   free(a);free(b);free(c);

 printf("thread done \n"); 	



}









int main(void) 

{ 

 printf("program start \n"); 

int i;



pthread_t handle[THREAD_NUM];



for( i=0;i<THREAD_NUM;i++)

pthread_create(&handle[i],NULL,(void*)thread_func,NULL);

 

for(i=0;i<THREAD_NUM;i++)

	pthread_join(handle[i],NULL);



perror("a*b leaveing");

	return 0; 

}