|
|
|
@ -1,11 +1,3 @@
|
|
|
|
|
/*
|
|
|
|
|
License: MIT License
|
|
|
|
|
|
|
|
|
|
Changelog:
|
|
|
|
|
- 2023-03-31 Initial version by Sebastian Apel (https://github.com/SebastianApel)
|
|
|
|
|
|
|
|
|
|
*/
|
|
|
|
|
|
|
|
|
|
#include <locale.h>
|
|
|
|
|
#include "ggml.h"
|
|
|
|
|
#include <assert.h>
|
|
|
|
@ -45,7 +37,7 @@ float tensor_sum_elements(struct ggml_tensor * tensor) {
|
|
|
|
|
|
|
|
|
|
#define TENSOR_TYPE_AS_STR(TYPE) TYPE == GGML_TYPE_F32 ? "FP32" : TYPE == GGML_TYPE_F16 ? "FP16" : TYPE == GGML_TYPE_Q4_0 ? "Q4_0" : TYPE == GGML_TYPE_Q4_1 ? "Q4_1" : "UNKNOWN"
|
|
|
|
|
|
|
|
|
|
#define TENSOR_DUMP(TENSOR) printf("%15s: type = %i (%5s) ne = %5d x %5d x %5d, nb = (%5li, %5li, %5li) - ", #TENSOR, \
|
|
|
|
|
#define TENSOR_DUMP(TENSOR) printf("%15s: type = %i (%5s) ne = %5ld x %5ld x %5ld, nb = (%5li, %5li, %5li) - ", #TENSOR, \
|
|
|
|
|
TENSOR->type,TENSOR_TYPE_AS_STR(TENSOR->type),\
|
|
|
|
|
TENSOR->ne[0], TENSOR->ne[1], TENSOR->ne[2], TENSOR->nb[0], TENSOR->nb[1], TENSOR->nb[2]); \
|
|
|
|
|
{ float sum = tensor_sum_elements(TENSOR); printf("Sum of tensor %s is %6.2f\n",#TENSOR, sum); }
|
|
|
|
@ -98,12 +90,9 @@ int main(int argc, char ** argv) {
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
// create the ggml context
|
|
|
|
|
printf("Starting Test\n");
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
struct ggml_context * ctx;
|
|
|
|
|
//const int sizex = 4096;
|
|
|
|
|
//const int sizey = 11008;
|
|
|
|
@ -125,16 +114,18 @@ int main(int argc, char ** argv) {
|
|
|
|
|
#endif
|
|
|
|
|
|
|
|
|
|
//printf("Memsize required = %i\n", sizex*sizex);
|
|
|
|
|
ggml_type wtype = GGML_TYPE_F32;
|
|
|
|
|
|
|
|
|
|
size_t ctx_size = 0;
|
|
|
|
|
ctx_size += sizex*sizey*ggml_type_sizef(wtype);
|
|
|
|
|
ctx_size += sizex*sizey*ggml_type_sizef(wtype);
|
|
|
|
|
ctx_size += sizex*sizey*ggml_type_sizef(GGML_TYPE_F32);
|
|
|
|
|
ctx_size += sizex*sizeof(float);
|
|
|
|
|
ctx_size += 1024*1024*100;
|
|
|
|
|
ctx_size += sizex*sizey*ggml_type_sizef(GGML_TYPE_F32);
|
|
|
|
|
ctx_size += sizex*sizez*ggml_type_sizef(GGML_TYPE_F32);
|
|
|
|
|
ctx_size += sizex*sizey*ggml_type_sizef(GGML_TYPE_Q4_0);
|
|
|
|
|
ctx_size += sizex*sizey*ggml_type_sizef(GGML_TYPE_Q4_0);
|
|
|
|
|
ctx_size += sizex*sizey*ggml_type_sizef(GGML_TYPE_F32); // BLAS
|
|
|
|
|
ctx_size += sizex*sizey*ggml_type_sizef(GGML_TYPE_F32); // BLAS
|
|
|
|
|
ctx_size += 1024*1024*16;
|
|
|
|
|
|
|
|
|
|
printf("Allocating Memory of size %li byes, %li MB\n",ctx_size, (ctx_size/1024/1024));
|
|
|
|
|
printf("Allocating Memory of size %li bytes, %li MB\n",ctx_size, (ctx_size/1024/1024));
|
|
|
|
|
|
|
|
|
|
struct ggml_init_params params = {
|
|
|
|
|
/*.mem_size =*/ ctx_size,
|
|
|
|
@ -217,7 +208,7 @@ int main(int argc, char ** argv) {
|
|
|
|
|
const int dimz = sizez;
|
|
|
|
|
long long int flops_per_dot_product = dimy + dimy;
|
|
|
|
|
long long int flops_per_matrix = flops_per_dot_product * dimx * dimz; ;
|
|
|
|
|
printf("Matrix Multiplication of (%i,%i,%i) x (%i,%i,%i) - aboout %6.2f gFLOPS\n\n", sizex, sizey, 1, sizex, sizez, 1, 1.0f*flops_per_matrix / 1000 / 1000 / 1000);
|
|
|
|
|
printf("Matrix Multiplication of (%i,%i,%i) x (%i,%i,%i) - about %6.2f gFLOPS\n\n", sizex, sizey, 1, sizex, sizez, 1, 1.0f*flops_per_matrix / 1000 / 1000 / 1000);
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
// Let's use the F32 result from above as a reference for the q4_0 multiplication
|
|
|
|
@ -234,7 +225,6 @@ int main(int argc, char ** argv) {
|
|
|
|
|
ggml_graph_compute(ctx, &gf31);
|
|
|
|
|
long long int stop = ggml_time_us();
|
|
|
|
|
long long int usec = stop-start;
|
|
|
|
|
float sec = usec/1000000;
|
|
|
|
|
float flops_per_usec = (1.0f*flops_per_matrix)/usec;
|
|
|
|
|
printf("%9i;%8i;%6i;%6i;%6i;%15lli;%18lli;%19.2f\n",
|
|
|
|
|
i,
|