Program: gdcm
Module: $RCSfile: TestInline.cxx,v $
Language: C++
- Date: $Date: 2005/11/18 11:46:37 $
- Version: $Revision: 1.9 $
+ Date: $Date: 2008/09/19 09:33:17 $
+ Version: $Revision: 1.21 $
Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
l'Image). All rights reserved. See Doc/License.txt or
// - with a macro : this is the quicker (any doubt ?)
// - with a function, passing the params by pointer
// - with a function, passing the params by reference (exactly same time)
-// - with an inline function described in the main()
+// - with an inline function described in the .cxx
// absolutely NO effect ?!?
// - with a function, described in the .h
// absolutely NO effect ?!?
// Which CXX_FLAGS, LINKER_FLAGS, ..., must we set to see the difference?
#include "gdcmUtil.h"
+#include <string.h> // Needed under Suse?!?
+#include <stdlib.h> // atoi
+#if defined(__BORLANDC__) || defined (_MSC_VER)
+#include <time.h>
+ #if defined(__BORLANDC__)
+ #include <stdio.h>
+ #endif
+#define GET_TIME(a) a=clock()
+#define HOW_LONG(b,a) \
+ std::cout << (double) (b-a) << std::endl
+
+#else
+#include <sys/times.h>
+#define GET_TIME(a) times(&a)
+#define HOW_LONG(b,a) \
+ std::cout \
+ << (long) ((b.tms_utime) - (a.tms_utime)) \
+ << std::endl
+#endif
#include <iostream>
inline void ifrswap(double &a, double &b);
inline void ifpswap(double *a, double *b);
-uint16_t passDirect(uint16_t a, uint16_t b);
-uint16_t passRef(uint16_t &a, uint16_t &b);
-uint16_t passPtr(uint16_t *a, uint16_t *b);
+uint8_t passDirect8(uint8_t a, uint8_t b);
+uint8_t passRef8(uint8_t &a, uint8_t &b);
+uint8_t passPtr8(uint8_t *a, uint8_t *b);
+
+uint16_t passDirect16(uint16_t a, uint16_t b);
+uint16_t passRef16(uint16_t &a, uint16_t &b);
+uint16_t passPtr16(uint16_t *a, uint16_t *b);
uint32_t passDirect32(uint32_t a, uint32_t b);
uint32_t passRef32(uint32_t &a, uint32_t &b);
double passRef(double &a, double &b);
double passPtr(double *a, double *b);
-#define \
-mswap(a, b) \
-{ \
- tmp = a; \
- a = b; \
- b = tmp; \
+#define \
+mswap(a, b) \
+{ \
+ double tmp = a;\
+ a = b; \
+ b = tmp; \
}
+// ============= no inline
+
void frswap(double &a, double &b)
{
double tmp;
tmp = a;
a = b;
b = tmp;
-
}
void fpswap(double *a, double *b)
tmp = *a;
*a = *b;
*b = tmp;
-
}
+// ============= inline
+
inline void ifpswap(double *a, double *b)
{
double tmp;
}
-uint16_t passRef(uint16_t &a, uint16_t &b)
+uint16_t passRef16(uint16_t &a, uint16_t &b)
+{
+ return a + b;
+}
+uint16_t passPtr16(uint16_t *a, uint16_t *b)
+{
+ return *a + *b;
+}
+uint16_t passDirect16(uint16_t a, uint16_t b)
+{
+ return a + b;
+}
+
+uint8_t passRef8(uint8_t &a, uint8_t &b)
{
return a + b;
}
-uint16_t passPtr(uint16_t *a, uint16_t *b)
+uint8_t passPtr8(uint8_t *a, uint8_t *b)
{
return *a + *b;
}
-uint16_t passDirect(uint16_t a, uint16_t b)
+uint8_t passDirect8(uint8_t a, uint8_t b)
+{
+ return a + b;
+}
+
+float passRefFloat(float &a, float &b)
+{
+ return a + b;
+}
+float passPtrFloat(float *a, float *b)
+{
+ return *a + *b;
+}
+float passDirectFloat(float a, float b)
{
return a + b;
}
int TestInline(int argc, char *argv[])
{
+#if defined(__BORLANDC__) || defined (_MSC_VER)
+ clock_t tms1, tms2;
+#else
+ struct tms tms1, tms2;
+#endif
+
+// ====================================================================
+
+ std::cout << std::endl << std::endl
+ << "Just to be sure : sizes of native types" << std::endl
+ << "======================================="
+ << std::endl << std::endl;
+ // just to know, on every proc
+ std::cout << "Size of char " << sizeof(char) << std::endl;
+ std::cout << "Size of short int " << sizeof(short int) << std::endl;
+ std::cout << "Size of int " << sizeof(int) << std::endl;
+ std::cout << "Size of long " << sizeof(long) << std::endl;
+ std::cout << "Size of float " << sizeof(float) << std::endl;
+ std::cout << "Size of double " << sizeof(double) << std::endl;
+ std::cout << "Size of bool " << sizeof(bool) << std::endl;
+ std::cout << std::endl;
+ std::cout << "Size of char* " << sizeof(char*) << std::endl;
+ std::cout << "Size of short int*" << sizeof(short int*)<< std::endl;
+ std::cout << "Size of int* " << sizeof(int*) << std::endl;
+ std::cout << "Size of double* " << sizeof(double*) << std::endl;
+ std::cout << "-----------------" << std::endl;
+
+ // ====================================================================
+
unsigned int nbLoop;
unsigned int i;
if (argc > 1)
nbLoop = atoi(argv[1]);
else
- nbLoop = 10000000;
+ nbLoop = 100000000;
- //clock_t r1, r2;
- struct tms tms1, tms2;
-
+ uint8_t x8 =1, y8 =2;
+ uint16_t x16=1, y16=2;
+ uint32_t x32=1, y32=2;
+ float fx =1.0f, fy=1.0f;
+ double dx =1.0 , dy=1.0;
double a = 1, b = 2;
- double tmp;
- uint16_t x=1, y=2;
+ // ====================================================================
+
+ std::cout << std::endl << std::endl
+ << "Check different ways of passing scalars to a function "<< nbLoop << " times" << std::endl
+ << "====================================================="
+ << std::endl << std::endl;
+
+ std::cout << "Pass uint_8 param directly"
+ << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ passDirect8 (x8, y8);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
// ----------------------------------------
- std::cout << "Use a macro "<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
+ std::cout << "Pass uint_8 param as ref"
+ << std::endl;
+
+ GET_TIME(tms1);
for(i = 0 ; i< nbLoop ; i++)
{
- mswap (a,b);
+ passRef8 (x8, y8);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
-
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
// ----------------------------------------
- std::cout << "Use reference function" << std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+ std::cout << "Pass uint_8 param as ptr"
+ << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- frswap (a,b);
+ passPtr8 (&x8, &y8);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
-
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
// ----------------------------------------
-
- std::cout << "Use pointer function" << std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+ std::cout << std::endl;
+ std::cout << "Pass uint_16 param directly"
+ << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- fpswap (&a, &b);
+ passDirect16 (x16, y16);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
-
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
// ----------------------------------------
- std::cout << "Use inline, main-defined reference function" << std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+ std::cout << "Pass uint_16 param as ref"
+ << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- ifrswap (a, b);
+ passRef16 (x16, y16);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
-
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
// ----------------------------------------
- std::cout << "Use inline, main-defined pointer function" << std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+ std::cout << "Pass uint_16 param as ptr"
+ << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- ifpswap (&a, &b);
+ passPtr16 (&x16, &y16);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
// ----------------------------------------
-
-//To check the 2 following cases, we just put the two 'static' functions
-//hifpswap and hNoifpswap in gdcmUtil.h
-
- std::cout << "Use inline, .h defined, WITH inline keyword pointer function"
+ std::cout << std::endl;
+ std::cout << "Pass uint_32 param directly"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- //gdcm::Util::hifpswap (&a, &b);
- hifpswap (&a, &b);
+ passDirect32 (x32, y32);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
-
// ----------------------------------------
-
- std::cout << "Use inline, .h defined, NO inline keyword pointer function"
+
+ std::cout << "Pass uint32_t param as ref"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- //gdcm::Util::hNoifpswap (&a, &b);
- hNoifpswap (&a, &b);
+ passRef32 (x32, y32 );
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
-
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
// ----------------------------------------
- std::cout << "Pass uint_16 param directly"
+ std::cout << "Pass uint_32 param as ptr"
+ << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ passPtr32 (&x32, &y32);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+ // ----------------------------------------
+ std::cout << std::endl;
+ std::cout << "Pass float param directly"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passDirect (x, y);
+ passDirectFloat (fx, fy);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
// ----------------------------------------
- std::cout << "Pass uint_16 param as ref"
+ std::cout << "Pass float param as ref"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passRef (x, y);
+ passRefFloat (fx, fy);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
// ----------------------------------------
- std::cout << "Pass uint_16 param as ptr"
+ std::cout << "Pass float param as ptr"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passPtr (&x, &y);
+ passPtrFloat (&fx, &fy);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+ // ----------------------------------------
+ std::cout << std::endl;
+ std::cout << "Pass double param directly"
+ << std::endl;
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ passDirectDouble (dx, dy);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
// ----------------------------------------
+
+ std::cout << "Pass double param as ref"
+ << std::endl;
- uint32_t m =1, n=2;
- std::cout << "Pass uint_32 param directly"
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ passRefDouble (dx, dy);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+ // ----------------------------------------
+
+ std::cout << "Pass double param as ptr"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passDirect32 (m, n);
+ passPtrDouble (&dx, &dy);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+// ====================================================================
+
+ std::cout << std::endl;
+ std::cout << "Exchange 2 scalars " << nbLoop << " times" << std::endl
+ << "==================="
+ << std::endl << std::endl;
+
// ----------------------------------------
- std::cout << "Pass uint32_t param as ref"
- << std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+ std::cout << "Direct "<< std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passRef32 (m, n);
+ double tmp;
+ tmp=a;
+ a=b;
+ b=tmp;
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
// ----------------------------------------
- std::cout << "Pass uint_32 param as ptr"
- << std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+ std::cout << "Use a macro "<< std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passPtr32 (&m, &n);
+ mswap (a,b);
}
- //r2 = times(&tms2);
- times(&tms2);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+ // ----------------------------------------
+
+ std::cout << std::endl;
+ std::cout << "Use a function, param passed by reference" << std::endl;
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ frswap (a,b);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+ // ----------------------------------------
+
+ std::cout << "Use a function, param passed by pointer" << std::endl;
+
+ GET_TIME(tms1); ;
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ fpswap (&a, &b);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
// ----------------------------------------
+ std::cout << std::endl;
- double dx=1.0, dy=2.0;
+ std::cout << "Use inline, .cxx-defined function, param passed by reference" << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ ifrswap (a, b);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
- std::cout << "Pass double param directly"
+ // ----------------------------------------
+
+ std::cout << "Use inline, .cxx-defined function, param passed by pointer" << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ ifpswap (&a, &b);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+ // ----------------------------------------
+
+ std::cout << std::endl;
+
+//To check the 2 following cases, we just put the two 'static' functions
+//hifpswap and hNoifpswap in gdcmUtil.h
+
+ std::cout << "Use inline, .h defined, WITH inline keyword, param passed by pointer"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GDCM_NAME_SPACE::Util util;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passDirectDouble (dx, dy);
+ util.hifpswap (&a, &b);
}
- //r2 = times(&tms2);
- times(&tms2);
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+
+ // ----------------------------------------
+
+ std::cout << "Use inline, .h defined, NO inline keyword, param passed by pointer"
+ << std::endl;
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ util.hNoifpswap (&a, &b);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
// ----------------------------------------
-
- std::cout << "Pass double param as ref"
+ std::cout << std::endl;
+
+ std::cout << "Use inline, .h defined, WITH inline keyword, param passed by pointer STATIC function"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passRefDouble (dx, dy);
+ GDCM_NAME_SPACE::Util::sthifpswap (&a, &b);
}
- //r2 = times(&tms2);
- times(&tms2);
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+ // ----------------------------------------
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+ std::cout << "Use inline, .h defined, NO inline keyword, param passed by pointer STATIC function"
+ << std::endl;
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ GDCM_NAME_SPACE::Util::sthNoifpswap (&a, &b);
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
// ----------------------------------------
- std::cout << "Pass double param as ptr"
+ // Just to point out that playing with pointers doesn't save so much time ...
+
+ std::cout << "Play with arrays\n================" << std::endl;
+
+ nbLoop=1000;
+
+ std::cout << "Copy 2 arrays [i][j]"
+ << std::endl;
+
+
+ unsigned short int z1[128][3118], z2[128][3118];
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ unsigned short int *pv1=&z1[0][0], *pv2=&z2[0][0];
+ for (int j=0;j<3118;j++)
+ for(int i=0; i<128;i++)
+ z2[i][j] = z1[i][j];
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+ std::cout << "Copy 2 arrays ([i][j], pointer)"
<< std::endl;
- //r1 = times(&tms1);
- times(&tms1);
- for(i = 0 ; i< nbLoop ; i++)
+
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
{
- passPtrDouble (&dx, &dy);
- }
- //r2 = times(&tms2);
- times(&tms2);
+ unsigned short int *pv1=&z1[0][0], *pv2=&z2[0][0];
+ for (int j=0;j<3118;j++)
+ for(int i=0; i<128;i++)
+ z2[i][j] = *pv1++;
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+ std::cout << "Copy 2 arrays (2 pointers)" << std::endl;
+ GET_TIME(tms1);
+ // unsigned short int w1[3118*128], w2[3118][128];
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ unsigned short int *pw1=&z1[0][0], *pw2=&z2[0][0];
+ for (int j=0;j<3118;j++)
+ for(int i=0; i<128;i++)
+ *pw2++ = *pw1++;
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+
+ std::cout << "Copy 2 arrays (memcpy)" << std::endl;
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ unsigned short int *pw1=&z1[0][0], *pw2=&z2[0][0];
+ memcpy(pw2,pw1,3118*128*sizeof(unsigned short int));
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+
+
+ std::cout << "Transpose 2 arrays [i][j]"
+ << std::endl;
+
+ unsigned short int t1[3118][128], t2[128][3118];
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ unsigned short int *pv1=&t1[0][0], *pv2=&t2[0][0];
+ for (int j=0;j<3118;j++)
+ for(int i=0; i<128;i++)
+ t2[i][j] = t1[j][i];
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+
+ std::cout << "Transpose 2 arrays ([i][j], pointer)"
+ << std::endl;
+
+ unsigned short int w1[3118*128], w2[3118][128];
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ unsigned short int *pw1=w1, *pw2=&w2[0][0];
+ for (int j=0;j<3118;j++)
+ for(int i=0; i<128;i++)
+ w2[i][j] = *pw1++;
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
+
+
+ std::cout << "Transpose 2 arrays (2 pointers)"
+ << std::endl;
+
+ unsigned short int v1[3118*128], v2[128*3118];
+ GET_TIME(tms1);
+ for(i = 0 ; i< nbLoop ; i++)
+ {
+ unsigned short int *pv1=v1, *pv2=v2;
+ for (int j=0;j<3118;j++)
+ for(int i=0; i<128;i++)
+ *(pv2+i*128+j) = *pv1++;
+ }
+ GET_TIME(tms2);
+ HOW_LONG(tms2,tms1);
- std::cout
- << (long) ((tms2.tms_utime) - (tms1.tms_utime))
- << std::endl;
+
+ //return 1; // will generate an error,
+ // just to allow us to see the full log in the dashboard
return 0;
}
git://git.creatis.insa-lyon.fr / gdcm.git / blobdiff