Source Code 1994 March

home *** CD-ROM | disk | FTP | other *** search

/ Source Code 1994 March / Source_Code_CD-ROM_Walnut_Creek_March_1994.iso / compsrcs / unix / volume26 / chernikv / part01 < prev next >

Wrap

Text File | 1993-04-02 | 22.7 KB | 625 lines

Newsgroups: comp.sources.unix From: atae@spva.physics.imperial.ac.uk (Dr Ata Etemadi) Subject: v26i091: chernikov - computes Chernikov stochastic webs, Part01/01 Sender: unix-sources-moderator@vix.com Approved: paul@vix.com Submitted-By: atae@spva.physics.imperial.ac.uk (Dr Ata Etemadi) Posting-Number: Volume 26, Issue 91 Archive-Name: chernikov/part01 This program computes the stochastic webs produced by the Chernikov equations (see Nature Vol. 326, April 1987) and produces a PGM image based on occupancy of cells. The equations essentially describe the path of a non-relativistic charged particle rotating about a magnetic field line, and experiencing a periodic electric field impulse. atae@spva.physics.imperial.ac.uk (Dr Ata Etemadi) #! /bin/sh # This is a shell archive. Remove anything before this line, then unpack # it by saving it into a file and typing "sh file". To overwrite existing # files, type "sh file -c". You can also feed this as standard input via # unshar, or by typing "sh <file", e.g.. If this archive is complete, you # will see the following message at the end: # "End of archive 1 (of 1)." # Contents: COPYRIGHT Chernikov.c LICENCE Makefile README # Wrapped by atae@minie on Sat Apr 18 23:45:25 1992 PATH=/bin:/usr/bin:/usr/ucb ; export PATH if test -f 'COPYRIGHT' -a "${1}" != "-c" ; then echo shar: Will not clobber existing file \"'COPYRIGHT'\" else echo shar: Extracting \"'COPYRIGHT'\" $995 characters$ sed "s/^X//" >'COPYRIGHT' <<'END_OF_FILE' X X Chernikov X XCopyright (C) 1992 X X Ataollah Etemadi, X Space and Atmospheric Physics Group, X Blackett Laboratory, X Imperial College of Science, Technology and Medicine X Prince Consort Road, X London SW7 2BZ. X X This program is free software; you can redistribute it and/or X modify it under the terms of the GNU General Public License as X published by the Free Software Foundation; either version 1, or X any later version. X X This program is distributed in the hope that it will be useful, X but WITHOUT ANY WARRANTY; without even the implied warranty of X MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the X GNU General Public License for more details. X X You should have received a copy of the GNU General Public License X along with this program; if not, write to the Free Software X Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. X X Note: This program has been exempted from the requirement of X paragraph 2c of the General Public License. END_OF_FILE if test 995 -ne `wc -c <'COPYRIGHT'`; then echo shar: \"'COPYRIGHT'\" unpacked with wrong size! fi # end of 'COPYRIGHT' fi if test -f 'Chernikov.c' -a "${1}" != "-c" ; then echo shar: Will not clobber existing file \"'Chernikov.c'\" else echo shar: Extracting \"'Chernikov.c'\" $4245 characters$ sed "s/^X//" >'Chernikov.c' <<'END_OF_FILE' X/*++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ X XName : Chernikov XVersion: 1.0 XWritten on : 24-Apr-92 By : A. Etemadi XModified on : XDirectory : ~atae/Space X X============================================================================== X X XUsage: X Chernikov -[AKUVhsnrc] > OutputPGMImage X XOutput result : X X 0 = successful, X -1 = error, X XFunctionality: X XIterates through the Chernikov equations (Nature Vol. 326, April 1987) Xand produces an image based on occupancy of cells. X X----------------------------------------------------------------------------*/ X X#include <stdio.h> /* Standard C I/O library */ X#include <math.h> /* Standard C mathematics library */ X X#define pi 2.0*acos(0.0) X Xmain(argc,argv) X X int argc; X char **argv; X X{ X X/* X======================================================================= X===================== START OF DECLARATION ============================ X======================================================================= X*/ X X/* X Indecies X*/ X X register int i,j; X register int row,col; X X int Height,Width; X X FILE *p_OutFile; /* Pointer to start of output file */ X X/* X Create buffer for image X*/ X char OutImage[512][512]; X X/* X Chernikov parameters X*/ X float Alpha; X float K; X float U0,U1; X float V0,V1; X float Scale; X X int NIter; X X float KoverAlpha; X float cosAlpha; X float sinAlpha; X float Beta; X/* X======================================================================= X===================== START OF INITIALISATION ========================= X======================================================================= X*/ X X/* X Set defaults, process command line options X*/ X Alpha = 2.0*pi/5.0; X K = 0.9; X U0 = 0.0; X V0 = 15.0; X Scale = 1.5; X NIter = 10000; X Height = Width = 256; X X for (i=1;i<argc;i++) { X if (argv[i][0] == '-') { X switch (argv[i][1]) { X case 'A': Alpha = atof(argv[++i]); continue; X case 'K': K = atof(argv[++i]); continue; X case 'U': U0 = atof(argv[++i]); continue; X case 'V': V0 = atof(argv[++i]); continue; X case 's': Scale = atof(argv[++i]); continue; X case 'n': NIter = atoi(argv[++i]); continue; X case 'r': Height = atoi(argv[++i]); continue; X case 'c': Width = atoi(argv[++i]); continue; X default : fprintf(stderr,"unrecognized option: %s",argv[i]); return(-1); X case 'h': X fprintf(stderr,"\nUSAGE :\n"); X fprintf(stderr," Chernikov -[AKUVhsnrc] > OutImage\n\n"); X fprintf(stderr,"OPTION: DEFAULT:\n"); X fprintf(stderr,"-----------------------------------\n"); X fprintf(stderr,"-A Alpha 2pi/15\n"); X fprintf(stderr,"-K Kappa 0.9\n"); X fprintf(stderr,"-U Initial pahse 0.0\n"); X fprintf(stderr,"-V Initial Velocity 15.0\n"); X fprintf(stderr,"-s Scale 2.0\n"); X fprintf(stderr,"-n Interations 10000\n"); X fprintf(stderr,"-r Image Height 256\n"); X fprintf(stderr,"-c Image Width 256\n"); X return(-1); X } X } X } X X/* X Use standard output X*/ X p_OutFile = stdout; X/* X Initialise OutImage pixels and starting parameters X*/ X for (i = 0; i < Height; i++) { X for (j = 0; j < Width; j++) { X OutImage[i][j] = (char) 0; X } X } X X cosAlpha = cos(Alpha); X sinAlpha = sin(Alpha); X KoverAlpha = K/Alpha; X fprintf(stderr,"\n"); X X/* X Now iterate through the Chernikov equation X*/ X for (i = 0; i < NIter; i++) { X X row = (int) (U0*Scale) + Height/2; X col = (int) (V0*Scale) + Width/2; X X if ( (i+1)%10000 == 0) X fprintf(stderr,"Iteration = %d\r",i+1); X X if (row < Height && col < Width && row > 0 && col > 0 && X (OutImage[row][col] & 0377) < 255 ) X OutImage[row][col] += (char) 1; X X Beta = U0 + KoverAlpha*sin(V0); X U1 = Beta*cosAlpha + V0*sinAlpha; X V1 = -Beta*sinAlpha + V0*cosAlpha; X U0 = U1; V0 = V1; X } X/* X Write resultant image with PGM header X*/ X X fprintf(p_OutFile,"P5\n%d %d\n255\n",Height,Width); X for (i = 0; i < Height; i++) { X for (j = 0; j < Width; j++) { X putc(OutImage[i][j],p_OutFile); X } X } X X/* X Close file X*/ X fprintf(stderr,"\n"); X fclose(p_OutFile); X X return(0); X} END_OF_FILE if test 4245 -ne `wc -c <'Chernikov.c'`; then echo shar: \"'Chernikov.c'\" unpacked with wrong size! fi # end of 'Chernikov.c' fi if test -f 'LICENCE' -a "${1}" != "-c" ; then echo shar: Will not clobber existing file \"'LICENCE'\" else echo shar: Extracting \"'LICENCE'\" $12488 characters$ sed "s/^X//" >'LICENCE' <<'END_OF_FILE' X X GNU GENERAL PUBLIC LICENSE X Version 1, February 1989 X X Copyright (C) 1989 Free Software Foundation, Inc. X 675 Mass Ave, Cambridge, MA 02139, USA X Everyone is permitted to copy and distribute verbatim copies X of this license document, but changing it is not allowed. X X Preamble X X The license agreements of most software companies try to keep users Xat the mercy of those companies. By contrast, our General Public XLicense is intended to guarantee your freedom to share and change free Xsoftware--to make sure the software is free for all its users. The XGeneral Public License applies to the Free Software Foundation's Xsoftware and to any other program whose authors commit to using it. XYou can use it for your programs, too. X X When we speak of free software, we are referring to freedom, not Xprice. Specifically, the General Public License is designed to make Xsure that you have the freedom to give away or sell copies of free Xsoftware, that you receive source code or can get it if you want it, Xthat you can change the software or use pieces of it in new free Xprograms; and that you know you can do these things. X X To protect your rights, we need to make restrictions that forbid Xanyone to deny you these rights or to ask you to surrender the rights. XThese restrictions translate to certain responsibilities for you if you Xdistribute copies of the software, or if you modify it. X X For example, if you distribute copies of a such a program, whether Xgratis or for a fee, you must give the recipients all the rights that Xyou have. You must make sure that they, too, receive or can get the Xsource code. And you must tell them their rights. X X We protect your rights with two steps: (1) copyright the software, and X(2) offer you this license which gives you legal permission to copy, Xdistribute and/or modify the software. X X Also, for each author's protection and ours, we want to make certain Xthat everyone understands that there is no warranty for this free Xsoftware. If the software is modified by someone else and passed on, we Xwant its recipients to know that what they have is not the original, so Xthat any problems introduced by others will not reflect on the original Xauthors' reputations. X X The precise terms and conditions for copying, distribution and Xmodification follow. X X GNU GENERAL PUBLIC LICENSE X TERMS AND CONDITIONS FOR COPYING, DISTRIBUTION AND MODIFICATION X X 0. This License Agreement applies to any program or other work which Xcontains a notice placed by the copyright holder saying it may be Xdistributed under the terms of this General Public License. The X"Program", below, refers to any such program or work, and a "work based Xon the Program" means either the Program or any work containing the XProgram or a portion of it, either verbatim or with modifications. Each Xlicensee is addressed as "you". X X 1. You may copy and distribute verbatim copies of the Program's source Xcode as you receive it, in any medium, provided that you conspicuously and Xappropriately publish on each copy an appropriate copyright notice and Xdisclaimer of warranty; keep intact all the notices that refer to this XGeneral Public License and to the absence of any warranty; and give any Xother recipients of the Program a copy of this General Public License Xalong with the Program. You may charge a fee for the physical act of Xtransferring a copy. X X 2. You may modify your copy or copies of the Program or any portion of Xit, and copy and distribute such modifications under the terms of Paragraph X1 above, provided that you also do the following: X X a) cause the modified files to carry prominent notices stating that X you changed the files and the date of any change; and X X b) cause the whole of any work that you distribute or publish, that X in whole or in part contains the Program or any part thereof, either X with or without modifications, to be licensed at no charge to all X third parties under the terms of this General Public License (except X that you may choose to grant warranty protection to some or all X third parties, at your option). X X c) If the modified program normally reads commands interactively when X run, you must cause it, when started running for such interactive use X in the simplest and most usual way, to print or display an X announcement including an appropriate copyright notice and a notice X that there is no warranty (or else, saying that you provide a X warranty) and that users may redistribute the program under these X conditions, and telling the user how to view a copy of this General X Public License. X X d) You may charge a fee for the physical act of transferring a X copy, and you may at your option offer warranty protection in X exchange for a fee. X XMere aggregation of another independent work with the Program (or its Xderivative) on a volume of a storage or distribution medium does not bring Xthe other work under the scope of these terms. X X 3. You may copy and distribute the Program (or a portion or derivative of Xit, under Paragraph 2) in object code or executable form under the terms of XParagraphs 1 and 2 above provided that you also do one of the following: X X a) accompany it with the complete corresponding machine-readable X source code, which must be distributed under the terms of X Paragraphs 1 and 2 above; or, X X b) accompany it with a written offer, valid for at least three X years, to give any third party free (except for a nominal charge X for the cost of distribution) a complete machine-readable copy of the X corresponding source code, to be distributed under the terms of X Paragraphs 1 and 2 above; or, X X c) accompany it with the information you received as to where the X corresponding source code may be obtained. (This alternative is X allowed only for noncommercial distribution and only if you X received the program in object code or executable form alone.) X XSource code for a work means the preferred form of the work for making Xmodifications to it. For an executable file, complete source code means Xall the source code for all modules it contains; but, as a special Xexception, it need not include source code for modules which are standard Xlibraries that accompany the operating system on which the executable Xfile runs, or for standard header files or definitions files that Xaccompany that operating system. X X 4. You may not copy, modify, sublicense, distribute or transfer the XProgram except as expressly provided under this General Public License. XAny attempt otherwise to copy, modify, sublicense, distribute or transfer Xthe Program is void, and will automatically terminate your rights to use Xthe Program under this License. However, parties who have received Xcopies, or rights to use copies, from you under this General Public XLicense will not have their licenses terminated so long as such parties Xremain in full compliance. X X 5. By copying, distributing or modifying the Program (or any work based Xon the Program) you indicate your acceptance of this license to do so, Xand all its terms and conditions. X X 6. Each time you redistribute the Program (or any work based on the XProgram), the recipient automatically receives a license from the original Xlicensor to copy, distribute or modify the Program subject to these Xterms and conditions. You may not impose any further restrictions on the Xrecipients' exercise of the rights granted herein. X X 7. The Free Software Foundation may publish revised and/or new versions Xof the General Public License from time to time. Such new versions will Xbe similar in spirit to the present version, but may differ in detail to Xaddress new problems or concerns. X XEach version is given a distinguishing version number. If the Program Xspecifies a version number of the license which applies to it and "any Xlater version", you have the option of following the terms and conditions Xeither of that version or of any later version published by the Free XSoftware Foundation. If the Program does not specify a version number of Xthe license, you may choose any version ever published by the Free Software XFoundation. X X 8. If you wish to incorporate parts of the Program into other free Xprograms whose distribution conditions are different, write to the author Xto ask for permission. For software which is copyrighted by the Free XSoftware Foundation, write to the Free Software Foundation; we sometimes Xmake exceptions for this. Our decision will be guided by the two goals Xof preserving the free status of all derivatives of our free software and Xof promoting the sharing and reuse of software generally. X X NO WARRANTY X X 9. BECAUSE THE PROGRAM IS LICENSED FREE OF CHARGE, THERE IS NO WARRANTY XFOR THE PROGRAM, TO THE EXTENT PERMITTED BY APPLICABLE LAW. EXCEPT WHEN XOTHERWISE STATED IN WRITING THE COPYRIGHT HOLDERS AND/OR OTHER PARTIES XPROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESSED XOR IMPLIED, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF XMERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE ENTIRE RISK AS XTO THE QUALITY AND PERFORMANCE OF THE PROGRAM IS WITH YOU. SHOULD THE XPROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF ALL NECESSARY SERVICING, XREPAIR OR CORRECTION. X X 10. IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING XWILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MAY MODIFY AND/OR XREDISTRIBUTE THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, XINCLUDING ANY GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING XOUT OF THE USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED XTO LOSS OF DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY XYOU OR THIRD PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER XPROGRAMS), EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE XPOSSIBILITY OF SUCH DAMAGES. X X END OF TERMS AND CONDITIONS X X Appendix: How to Apply These Terms to Your New Programs X X If you develop a new program, and you want it to be of the greatest Xpossible use to humanity, the best way to achieve this is to make it Xfree software which everyone can redistribute and change under these Xterms. X X To do so, attach the following notices to the program. It is safest to Xattach them to the start of each source file to most effectively convey Xthe exclusion of warranty; and each file should have at least the X"copyright" line and a pointer to where the full notice is found. X X <one line to give the program's name and a brief idea of what it does.> X Copyright (C) 19yy <name of author> X X This program is free software; you can redistribute it and/or modify X it under the terms of the GNU General Public License as published by X the Free Software Foundation; either version 1, or (at your option) X any later version. X X This program is distributed in the hope that it will be useful, X but WITHOUT ANY WARRANTY; without even the implied warranty of X MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the X GNU General Public License for more details. X X You should have received a copy of the GNU General Public License X along with this program; if not, write to the Free Software X Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. X XAlso add information on how to contact you by electronic and paper mail. X XIf the program is interactive, make it output a short notice like this Xwhen it starts in an interactive mode: X X Gnomovision version 69, Copyright (C) 19xx name of author X Gnomovision comes with ABSOLUTELY NO WARRANTY; for details type `show w'. X This is free software, and you are welcome to redistribute it X under certain conditions; type `show c' for details. X XThe hypothetical commands `show w' and `show c' should show the Xappropriate parts of the General Public License. Of course, the Xcommands you use may be called something other than `show w' and `show Xc'; they could even be mouse-clicks or menu items--whatever suits your Xprogram. X XYou should also get your employer (if you work as a programmer) or your Xschool, if any, to sign a "copyright disclaimer" for the program, if Xnecessary. Here a sample; alter the names: X X Yoyodyne, Inc., hereby disclaims all copyright interest in the X program `Gnomovision' (a program to direct compilers to make passes X at assemblers) written by James Hacker. X X <signature of Ty Coon>, 1 April 1989 X Ty Coon, President of Vice X XThat's all there is to it! END_OF_FILE if test 12488 -ne `wc -c <'LICENCE'`; then echo shar: \"'LICENCE'\" unpacked with wrong size! fi # end of 'LICENCE' fi if test -f 'Makefile' -a "${1}" != "-c" ; then echo shar: Will not clobber existing file \"'Makefile'\" else echo shar: Extracting \"'Makefile'\" $649 characters$ sed "s/^X//" >'Makefile' <<'END_OF_FILE' X# Makefile for Chernikov X# Created by AE @ ICST 14-04-92 X# X XLIBS = -lm X XLIBDIR = /usr/lib XLIBFLG = -L${LIBDIR} X XINCDIR = /usr/include XINCFLG = -I${INCDIR} X XBINDIR = . XOBJDIR = . XSRCDIR = . X XOPT = -O X XBINARIES = ${BINDIR}/Chernikov X XOBJECTS = ${OBJDIR}/Chernikov.o X Xall: ${BINARIES} ${OBJECTS} X Xclean: X /bin/rm -fr *.*~* X Xsuperclean: X /bin/rm -fr *.*~* ${OBJECTS} X Xultraclean: X /bin/rm -fr *.*~* ${OBJECTS} ${BINARIES} X Xshar: X makekit * X X${BINDIR}/Chernikov : ${OBJDIR}/Chernikov.o X ${CC} ${OPT} ${OBJDIR}/Chernikov.o ${LIBFLG} ${LIBS} -o $@ X X${OBJDIR}/Chernikov.o: ${SRCDIR}/Chernikov.c X ${CC} -c ${CFLAGS} ${SRCDIR}/Chernikov.c ${INCFLG} -o $@ END_OF_FILE if test 649 -ne `wc -c <'Makefile'`; then echo shar: \"'Makefile'\" unpacked with wrong size! fi chmod +x 'Makefile' # end of 'Makefile' fi if test -f 'README' -a "${1}" != "-c" ; then echo shar: Will not clobber existing file \"'README'\" else echo shar: Extracting \"'README'\" $915 characters$ sed "s/^X//" >'README' <<'END_OF_FILE' XType 'make' to install in current directory. X XRunning the program with -h flag gives usage info. X XThis program computes the stochastic webs produced by the Chernikov Xequations (see Nature Vol. 326, April 1987) and produces a PGM image Xbased on occupancy of cells. The equations essentially describe the Xpath of a non-relativistic charged particle rotating about a magnetic Xfield line, and experiencing a periodic electric field impulse. You Xneed to read the paper to really understand the parameters, in brief Xhowever: X X Alpha is the period between impulses X K is the strength of the impulse it receives X U is the phase of the particle X V is the velocity of the particle X XThe other parameters are specific to the program namely: X X s is the scaling factor X n is the number of iterations X r is the number of rows in the resulting image X c is the number of columns in the resulting image X X adios amigos X Ata <(|)>. END_OF_FILE if test 915 -ne `wc -c <'README'`; then echo shar: \"'README'\" unpacked with wrong size! fi # end of 'README' fi echo shar: End of archive 1 $of 1$. cp /dev/null ark1isdone MISSING="" for I in 1 ; do if test ! -f ark${I}isdone ; then MISSING="${MISSING} ${I}" fi done if test "${MISSING}" = "" ; then echo You have the archive. rm -f ark[1-9]isdone else echo You still need to unpack the following archives: echo " " ${MISSING} fi ## End of shell archive. exit 0