Initial multithreaded version.

Author: jteresco

pthreads/Makefile

@@ -0,0 +1,18 @@
+#
+# Makefile for Jacobi iteration on an adaptive quadtree, pthreads version
+#
+# Jim Teresco, CSIS-335, Siena College, Fall 2021
+#
+CFILES = quadtree.c quadtree_solver.c com_abort.c
+OFILES = $(CFILES:.c=.o)
+CC=gcc
+CFLAGS=-g
+
+quadtree_solver:	$(OFILES)
+	$(CC) $(CFLAGS) -o quadtree_solver $(OFILES) -lm -lpthread
+
+quadtree.o:	quadtree.h macros.h
+quadtree_solver.o:	quadtree.h
+
+clean::
+	/bin/rm -f quadtree_solver $(OFILES)

pthreads/com_abort.c

@@ -0,0 +1,45 @@
+/* 
+   com_abort.c, based on version from
+
+   Parallel Communications Library
+
+   Jim Teresco
+
+   Rensselaer Polytechnic Institute
+   Department of Computer Science
+   Scientific Computation Research Center
+
+   Adapted from PMDB, JDT, Wed Nov 22 15:36:53 EST 1995
+
+   Last change: 
+   Mon Nov 27 15:08:24 EST 1995
+
+   $Id$
+   $Log: com_abort.c,v $
+   Revision 1.1  2006/03/06 05:28:23  terescoj
+   First prelim version
+
+   Revision 1.1.1.1  2002/06/15 01:04:29  terescoj
+   CVSing scorec util library
+
+   Revision 1.2  2000/08/14 19:02:23  oklaas
+   update to latest version of com
+
+   Revision 1.1.1.1  1996/12/18 14:53:28  jteresco
+   initial sources
+
+*/
+
+#include <stdio.h>
+#include <stdlib.h>
+
+void com_abort(char *function, char *msg) {
+
+  if (function || msg) {
+    fflush(stdout);
+    fprintf(stderr,"ABORT in %s!\n",function);
+    if (msg)
+      fprintf(stderr,"  ErrMsg: %s\n",msg);
+  }
+  abort();
+}

pthreads/macros.h

@@ -0,0 +1,129 @@
+/*
+   macros.h
+
+   Based on util_macros originally written by
+
+   Louis Ziantz
+   Jim Teresco
+
+   Rensselaer Polytechnic Institute
+   Department of Computer Science
+   Scientific Computation Research Center
+
+   Macros which are generally useful.
+
+   Created:
+   Wed Nov 29 17:05:11 EST 1995
+
+   Last change:
+   Mon Jan 22 13:02:18 EST 1996
+
+   $Id$
+
+   Modification History
+   07/16/1998  Added in FIELD_SET and OUTPUT_IN_ORDER from pmdb
+*/
+
+#ifndef __H_MACROS
+#define __H_MACROS
+
+#include <stdlib.h>
+#include <stdio.h>
+#include <errno.h>
+
+#ifndef FALSE
+#define FALSE 0
+#endif
+#ifndef TRUE
+#define TRUE  1
+#endif
+
+extern void com_abort(char *, char *);
+
+#define SAFE_MALLOC(v,type,size) \
+ {  v = (type) malloc(size) ; \
+    if ( v == NULL) { \
+       fflush(stdout); \
+       fprintf(stderr,"in file %s, line %d, failed to allocate %ld bytes",\
+               __FILE__,__LINE__,size); \
+       com_abort(NULL,NULL); \
+    } \
+ }
+
+#define SAFE_REALLOC(v,type,size) \
+ {  v = (type) realloc(v,size) ; \
+    if ( v == NULL) { \
+       fflush(stdout); \
+       fprintf(stderr,"in file %s, line %d, failed to reallocate %ld bytes",\
+               __FILE__,__LINE__,size); \
+       com_abort(NULL,NULL); \
+    } \
+ }
+
+#define GRACEFULLY_OPEN(fp, filename, mode) {                      \
+  fp = fopen(filename,mode) ;                                      \
+  if ( fp == NULL ) {                                              \
+    char msg[FILENAME_MAX+25];                                     \
+    sprintf(msg, "Could not open file %s.\n", filename);           \
+    perror(" fopen");                                              \
+    com_abort("GRACEFULLY_OPEN", msg);                             \
+  }                                                                \
+}
+
+#define FILE_EXISTS(fp,filename)                                   \
+        ( ((fp = fopen(filename,"r")) != NULL)                     \
+        ? (! fclose(fp)) /* TRUE */ : FALSE) 
+
+/*
+ * Macro for checking return from fprintf
+ *
+ * Should only be used if routine doesn't need to clean up before
+ * returning.  Calling routine should decare "int ret=0;" and do a
+ * "return(ret)" at the end.  This should not be used in main, as it
+ * will result in the program terminating and returning 0 on failure.
+ *
+ */
+
+#define SAFE_FPRINTF  if (ret==EOF) return(0); else ret=fprintf
+
+/*
+ * Macro for checking return from fclose
+ *
+ * Should only be used if routine doesn't need to clean up before
+ * returning.  Calling routine must return an int, 0 for failure.
+ * It is the caller's repsonsibility to abort if necessary.
+ *
+ */
+
+
+#define SAFE_FCLOSE(fp, fname) \
+  if (fclose(fp)) { \
+      fprintf(stderr,"error closing file %s\n",fname); \
+      perror(" fclose"); \
+      return(0); \
+    }
+
+#define GRACEFULLY_CLOSE(fp, fname) \
+  if (fclose(fp)) { \
+      char msg[FILENAME_MAX+25]; \
+      sprintf(msg,"Error closing file %s\n",fname); \
+      perror(" fclose"); \
+      com_abort("GRACEFULLY_CLOSE", msg); \
+    }
+
+#define ASSERT(cond) \
+  if (!(cond)) { \
+      fflush(stdout); \
+      fprintf(stderr,"ASSERT failed in file %s at line %d\n", \
+         __FILE__,__LINE__); \
+      com_abort(NULL,NULL); \
+  }
+
+#define FAIL \
+      fflush(stdout); \
+      fprintf(stderr,"FAIL in file %s at line %d\n", \
+         __FILE__,__LINE__); \
+      com_abort(NULL,NULL)
+
+/* do not add anything below this line! */
+#endif /* __H_MACROS */