add components

Author: James Michael DuPont

Makefile

@@ -1,10 +1,14 @@
 #
-all : strgcomp
+all : strgcomp componentsreadbin
 	echo done
 #bgl 
 #-save-temps
 bgl :
 	mpiCC -o bgl  bgl.cpp -lboost_system -lboost_mpi 
 
 strgcomp : strongcomponents.cpp
-	mpiCC -g -o strgcomp  strongcomponents.cpp -lboost_system  -lboost_graph
+	mpiCC -g -o strgcomp  strongcomponents.cpp -lboost_system  -lboost_graph
+
+componentsreadbin : componentsreadbin.c
+	gcc -o componentsreadbin componentsreadbin.c
+# use this to create a pack file perl process_components.pl

componentsreadbin.c

@@ -0,0 +1,31 @@
+#include <stdio.h>
+int main ()
+{
+  FILE * pFile;
+  pFile = fopen ( "components.bin" , "r" );
+  if (!pFile)   {
+    printf("file could not be opened\n");
+    return 2;
+  }
+  int nodecount;
+  int i=0;
+  int br=fread ((void*)&nodecount , sizeof(int), 1 , pFile ); // count of ways
+  printf("Reading bytes %d nodes %d",br,nodecount);
+  int nodecomponents[nodecount];
+  br= fread ((void*)&nodecomponents , sizeof(int), nodecount , pFile ); // index
+  fclose (pFile);
+
+
+  if (br != nodecount)
+    {
+      printf("Reading way wanted %d got %d\n",nodecount,br);
+      return 1;
+    }
+
+  for (i =0; i < nodecount; i++)    {
+    printf("%d -> %d\n",i, (int)nodecomponents[i]);
+  }
+
+
+}
+

process_components.pl

@@ -0,0 +1,22 @@
+#my %components;
+my @components;
+
+open IN,"components.txt";
+while(<IN>)
+{
+    if (/(\d+)\t(\d+)/)
+    {
+#        my $nodeindex=$1;
+        my $component=$2;
+        push @components,$2;
+#        $components{$nodeindex}=$component;
+    }
+}
+close IN;
+
+open(my $out, '>:raw', 'components.bin') or die "Unable to open: $!";
+my $count=scalar(@components);
+warn "Count is $count";
+print $out pack ("L",$count);
+print $out pack ("L*",@components);
+close $out;