cc-1098 cc: ERROR File = test_quad.c, Line = 7
An array cannot have elements of the indicated type.
void matrix_mult(long double A[][], long double B[][], long double C[][], int n);
^
cc-1098 cc: ERROR File = test_quad.c, Line = 7
An array cannot have elements of the indicated type.
void matrix_mult(long double A[][], long double B[][], long double C[][], int n);
^
cc-1098 cc: ERROR File = test_quad.c, Line = 7
An array cannot have elements of the indicated type.
void matrix_mult(long double A[][], long double B[][], long double C[][], int n);
^
cc-1098 cc: ERROR File = test_quad.c, Line = 8
An array cannot have elements of the indicated type.
void matrix_print(long double M[][], FILE * nout, int n);
^
cc-1098 cc: ERROR File = test_quad.c, Line = 9
An array cannot have elements of the indicated type.
void matrix_rand(long double M[][], int n);
^
cc-1241 cc: ERROR File = test_quad.c, Line = 27
A declaration cannot appear after an executable statement in a block.
long double A[MAX_DIM][MAX_DIM];
^
cc-1241 cc: ERROR File = test_quad.c, Line = 28
A declaration cannot appear after an executable statement in a block.
long double B[MAX_DIM][MAX_DIM];
^
cc-1241 cc: ERROR File = test_quad.c, Line = 29
A declaration cannot appear after an executable statement in a block.
long double C[MAX_DIM][MAX_DIM];
^
cc-1551 cc: WARNING File = test_quad.c, Line = 34
The variable "t1" is used before its value is set.
t1 += (double)tp.tv_sec+(1.e-6)*tp.tv_usec;
^
cc-1551 cc: WARNING File = test_quad.c, Line = 37
The variable "t2" is used before its value is set.
t2 += (double)tp.tv_sec+(1.e-6)*tp.tv_usec;
^
cc-1552 cc: WARNING File = test_quad.c, Line = 14
The variable "rtn" is set but never used.
int rtn, i;
^
cc-1174 cc: WARNING File = test_quad.c, Line = 14
The variable "i" was declared but never referenced.
int rtn, i;
~:cc -64 -IPA -Ofast=IP30 test_quad.c -o test_quad
cc-1178 cc: WARNING File = test_quad.c, Line = 34
Argument is incompatible with the corresponding format string conversion.
printf("size of long double: %d\n", sizeof(long double));
^
~:test_quad
Using default value, n = 128
size of long double: 16
Elapsed time to multiply two matrices of order 128: 0.242758
Code: Select all
#include <math.h>
#include <stdio.h>
#include <sys/time.h>
#include <stdlib.h>
#define MAX_DIM 256
void matrix_mult(long double A[][MAX_DIM], long double B[][MAX_DIM], long double C[][MAX_DIM], int n);
void matrix_print(long double M[][MAX_DIM], FILE * nout, int n);
void matrix_rand(long double M[][MAX_DIM], int n);
int main(int argc, char ** argv){
int n;
double t1=0.0, t2=0.0, elapsed;
timespec_t tp[2];
FILE *nout;
long double A[MAX_DIM][MAX_DIM];
long double B[MAX_DIM][MAX_DIM];
long double C[MAX_DIM][MAX_DIM];
memset( C, 0, MAX_DIM*MAX_DIM*sizeof(long double) );
nout = fopen("/dev/null", "r+");
if ( !nout ) {
fprintf(stderr,"open failed\n");
exit(-1);
}
n = 128; /* default order of matrix */
if(argc < 2)
printf("Using default value, n = 128\n");
else{
n = atoi(argv[1]);
printf("Using n = %d\n", n);
}
printf("size of long double: %d\n", sizeof(long double));
matrix_rand(A, n);
matrix_rand(B, n);
clock_gettime(CLOCK_SGI_CYCLE, &tp[0]);
matrix_mult(A, B, C, n);
clock_gettime(CLOCK_SGI_CYCLE, &tp[1]);
t1 += (double)tp[0].tv_sec+(1.e-9)*tp[0].tv_nsec;
t2 += (double)tp[1].tv_sec+(1.e-9)*tp[1].tv_nsec;
matrix_print(C, nout, n);
elapsed = t2 -t1;
printf("\nElapsed time to multiply two matrices of order %d: %f\n", n, elapsed);
return(0);
}
/*****************************************************************************/
void matrix_mult(long double A[][MAX_DIM], long double B[][MAX_DIM], long double C[][MAX_DIM], int n){
int i, j, k;
for(i = 0; i < n; i++)
for(j = 0; j < n; j++)
for(k = 0; k < n; k++)
C[i][j] += A[i][k]*B[k][j];
}
/*****************************************************************************/
void matrix_print(long double M[][MAX_DIM], FILE * nout, int n){
int i, j;
for(i = 0; i < n; i++){
for(j = 0; j < n; j++){
fprintf(nout, "%5.5Lf ", M[i][j]);
}
}
}
/*****************************************************************************/
void matrix_rand(long double M[][MAX_DIM], int n){
int i, j;
for(i = 0; i < n; i++)
for(j = 0; j < n; j++)
M[i][j] = (double)rand()/12.3;
}
jan-jaap wrote: Algorithm optimization rules -- you can't beat that with compiler switches!
man INTRO_BLAS3 wrote: NOTE: SCSL supports two different C interfaces to the BLAS:
* The C interface described in this man page and in individual BLAS man
pages follows the same conventions used for the C interface to the
SCSL signal processing library.
* SCSL also supports the C interface to the legacy BLAS set forth by
the BLAS Technical Forum. This interface supports row-major storage
of multidimensional arrays; see INTRO_CBLAS(3S) for details.
Code: Select all
~:diff test_double.c test_double.orig.c
5,6c5,6
< /* #include <cblas.h> */
< #include <scsl_blas.h>
---
> #include <cblas.h>
>
42c42
< dgemm("N", "N", n, n, n, alpha, &a[0], n, &b[0], n, beta, &c[0], n);
---
> dgemm(NoTranspose, NoTranspose, n, n, n, alpha, &a[0], n, &b[0], n, beta, &c[0], n);
cc -fullwarn -O3 -Ofast=IP35 -IPA -o test_double test_double.c -lscs -lftn
test_double
Using default value, n = 128
Elapsed time to multiply two matrices of order 128: 0.003320
~:cc -fullwarn -O3 -Ofast=IP35 -IPA -o test_quad-double test_quad-double.c
~:test_quad-double Using default value, n = 128
size of double: 8
Elapsed time to multiply two matrices of order 128: 0.002489
cc -64 -fullwarn -O3 -Ofast=IP35 -IPA -o test_double test_double.c -lscs -lftn
test_double
Using default value, n = 128
Elapsed time to multiply two matrices of order 128: 0.003326
cc -n32 -fullwarn -O3 -Ofast=IP35 -IPA -o test_double test_double.c -lscs_mp -lftn
test_double
Using default value, n = 128
Elapsed time to multiply two matrices of order 128: 0.045474
cc -n32 -fullwarn -O3 -Ofast=IP35 -IPA -o test_double test_double.c -lscs -lftn
Using default value, n = 1024
Elapsed time to multiply two matrices of order 1024: 1.320950
cc -n32 -fullwarn -O3 -Ofast=IP35 -IPA -o test_double test_double.c -lscs_mp -lftn
test_double
Using default value, n = 1024
Elapsed time to multiply two matrices of order 1024: 0.941582
Code: Select all
xset +fp /usr/freeware/share/emacs/fonts & /usr/nekoware/bin/iconbar & /usr/nekoware/bin/xscreensaver-command -exit; /usr/nekoware/bin/xscreensaver -nosplash
).
joerg wrote: I hope you got something wrong.
Btw:
If i would like so see something in /current i would place it there.
)
nekonoko wrote:squeen wrote: I'll try and repackage neko_subversion with the lastest release using both a new snvserve bin location, neko_svnserve chkconfig flag, and rc2.d startup script if no one minds.
Did you ever get that pulled together? Would be an interesting addition!
GeneratriX wrote:
neko_gimp-2.2.12.tardist (which bombs out as soon I open an image)
neko_mozilla-1.8a5.tardist (which bombs out randomly)
BTW: ...I also have neko_dia-0.95.tardist installed, which seems to work just fine... so, is there any major bug of which I was not awared of, I mean, it was removed from nekoware?
schleusel wrote: uploaded to beta:
neko_dvdauthor-0.6.11.tardist - command line DVD authoring tools
neko_libdvdread-0.9.6.tardist - library for reading DVD-video images
neko_mjpegtools-1.8.0.tardist - tools for MJPEG/MPEG capture/editing/compression
neko_mpgtx-1.3.tardist - command line MPEG audio/video toolbox
neko_dvdplusrw_tools-6.1.tardist DVD mastering tools (does not interfere with neko_cdrtools. It mereley adds some additional tools)
neko_dvdstyler-1.4.tardist wxwindows based GUI frontend for the dirty dozen
I couldn't test it much yet. Creating a very simple project and letting it create the DVD structures and image of it at least looked promising though..
DVDstyler is built against wxgtk. I was hoping to be able to build it against wxmotif but that didn't work out due to apparently missing drag and drop support in wxmotif