Changeset 626 for trunk/user/pgcd

Timestamp:

Apr 29, 2019, 7:25:09 PM (6 years ago)

Author:

alain

Message:

This version has been tested on the sort multithreaded application
for TSAR_IOB architectures ranging from 1 to 8 clusters.
It fixes three bigs bugs:
1) the dev_ioc device API has been modified: the dev_ioc_sync_read()
and dev_ioc_sync_write() function use now extended pointers on the
kernel buffer to access a mapper stored in any cluster.
2) the hal_uspace API has been modified: the hal_copy_to_uspace()
and hal_copy_from_uspace() functions use now a (cxy,ptr) couple
to identify the target buffer (equivalent to an extended pointer.
3) an implementation bug has been fixed in the assembly code contained
in the hal_copy_to_uspace() and hal_copy_from_uspace() functions.

File:

• trunk/user/pgcd/pgcd.c (modified) (3 diffs)

trunk/user/pgcd/pgcd.c

@@ -12 +12 @@
 #include <almosmkh.h>
 
+#define   INSTRUMENTATION  0
+#define   IDBG             0
+
 /////////////////
 void main( void )
@@ -17 +20 @@
     int                opx;
     int                opy;
+    int                x;
+    int                y;
     unsigned long long cycle;
     unsigned int       cxy;
@@ -24 +29 @@
     get_core( &cxy , &lid );
 
-    printf( "\n\n[pgcd] starts on core[%x,%d] / cycle %d\n",
+    printf( "\n[pgcd] starts on core[%x,%d] / cycle %d\n\n",
     cxy , lid , (unsigned int)cycle ); 
 
-    while (1) 
+    // get operand X
+    printf("operand X = ");
+    opx = get_uint32();
+    printf("\n");
+
+    // get operand Y
+    printf("operand Y = ");
+    opy = get_uint32();
+    printf("\n");
+
+    // check operands 
+    if( (opx == 0) || (opy == 0) ) 
     {
-        printf("\n*******************\n");
-        printf("operand X = ");
-        opx = get_uint32();
-        printf("\n");
-        printf("operand Y = ");
-        opy = get_uint32();
-        printf("\n");
+        printf("\n[pgcd error] operands must be strictly positive\n");
+        exit(0);
+    } 
 
-        if( (opx == 0) || (opy == 0) ) 
-        {
-           printf("operands must be positive and larger than 0 => exit\n");
-           exit( 0 );
-        } 
-        else 
-        {
-            while (opx != opy) 
-            {
-                if(opx > opy)   opx = opx - opy;
-                else            opy = opy - opx;
-            }
-            printf("pgcd      = %d", opx);
-        }
+    // compute PGCD
+    x = opx;
+    y = opy;
+    while (x != y) 
+    {
+        if(x > y)   x = x - y;
+        else        y = y - x;
     }
+
+    // display result
+    printf("pgcd      = %d\n", x);
+
+#if INSTRUMENTATION
+
+    char   name[64];
+    char   path[128];
+
+    // build a file name from X and Y values
+    snprintf( name , 64 , "pgcd_%d_%d", opx, opy ); 
+
+    // build file pathname
+    snprintf( path , 128 , "home/%s" , name );
+
+#if IDBG
+idbg();
+#endif
+
+    // open file
+    FILE * stream = fopen( path , NULL );
+
+    if( stream == NULL )
+    {
+        printf("\n[pgcd error] cannot open instrumentation file <%s>\n", name );
+        exit(0);
+    }
+
+    printf("\n[pgcd] file %s successfully open\n", path);
+
+#if IDBG
+idbg();
+#endif
+
+    // register results to file
+    int ret = fprintf( stream , "pgcd( %d , %d ) = %d\n", opx, opy, x );
+
+    if( ret < 0 )
+    {
+        printf("\n[pgcd error] cannot write to instrumentation file <%s>\n", name );
+        exit(0);
+    }
+
+    printf("\n[pgcd] file %s successfully written\n", path);
+
+    display_mapper( path , 0 , 64 );
+
+    // close instrumentation file
+
+#if IDBG
+idbg();
+#endif
+
+    if( fclose( stream ) )
+    {
+        printf("\n[pgcd error] cannot close the file <%s>\n", name );
+        exit(0);
+    }
+
+    printf("\n[pgcd] file %s successfully closed\n", path);
+
+#endif
+
+    exit(0);
+        
 } // end pgcd