Physics 329: Introduction to Computational Physics: Utilities

• Makefile for the demonstration (driver) programs.
• p329f.f: Source code for 'lib329f.a' library routines. p329fsa.f: Sample driver for utility routines.
• fsleep.c: Code for C routine fsleep which sleeps for specified number of seconds---the routine accepts a floating point value, so fsleep(0.01), for example, sleeps for 1/100th of a second. tfsleep.c: Driver program for routine fsleep.
• sphplot.c and sphplot.h: Code for routine sphplot which plots a collection of spheres using the GL SGI library.
  • Sample Fortran usage

       real*8    positions(3,max_ncharge),   radius
       integer   ncharge
    
    c---------------------------------------------------------
    c  Set the radius for *all* displayed spheres.
    c---------------------------------------------------------
       radius = 0.03d0 
    
    c---------------------------------------------------------
    c  Call with ncharge > 0 for normal display.  A graphics 
    c  window will be opened interactively the first time 
    c  this routine is called in a program (the outline of
    c  a small window will appear on the screen; you will have 
    c  to position and size the window).
    c---------------------------------------------------------
       call sphplot(positions,ncharge,radius)
    
    c---------------------------------------------------------
    c  Call with ncharge < 0 (presumably at the end of the 
    c  simulation) for "tumbling" display.  Depress the escape 
    c  key in the graphics window to return from the routine.
    c---------------------------------------------------------
       call sphplot(positions,-ncharge,radius)
         

  • Sample linking and loading with a Fortran main program pgm.f

    % f77 -c pgm.f 
    % cc -c sphplot.c
    % f77 pgm.o sphplot.o -lsphere -lgl_s -o pgm