GNU DDD 3.3.9 (i686-pld-linux-gnu)
Copyright (C) 1995-1999 Technische Universität Braunschweig, Germany.
Copyright (C) 1999-2001 Universität Passau, Germany.
Copyright (C) 2001 Universität des Saarlandes, Germany.
Copyright (C) 2001-2004 Free Software Foundation, Inc.

Compiled with GCC 3.3.4 (PLD Linux), GNU libc 2.3
Requires X11R6, Xt11R6, Motif 2.2.3 (Motif Version 2.2.3)
Includes XPM 3.4.11, Athena Panner (7000002L), DDD core
Built 2004-07-08 by PLD  builder <address@hidden>.
$  ddd
+  /bin/sh -c 'exec gdb -q -fullname'
#  Hello, world!
#  Starting GDB...
#  Enabling core dumps...
#  Enabling core dumps...done.
#  Running GDB (pid 3402, tty /dev/pts/35)...
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "set prompt (gdb) \n"
#  Current language: C/C++
#  Searching "vsllib/ddd.vsl"...
#  Trying "/home/zbyszek/.ddd/vsllib/ddd.vsl"
#  Trying "/usr/share/ddd-3.3.9/vsllib/ddd.vsl"
#  Searching "vsllib/ddd.vsl"..."/usr/share/ddd-3.3.9/vsllib/ddd.vsl".
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "set height 0\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "set width 0\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "set annotate 1\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> " set verbose off\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "info line\n"
2004.07.14 03:27:05
<- "No line number information available.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "list\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line\n"
2004.07.14 03:27:05
<- "No line number information available.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "output 4711\n"
2004.07.14 03:27:05
<- "4711(gdb) "
2004.07.14 03:27:05
-> "show language\n"
2004.07.14 03:27:05
<- "The current source language is \"auto; currently c\".\n"
   "(gdb) "
2004.07.14 03:27:05
-> "show version\n"
2004.07.14 03:27:05
<- "GNU gdb 6.1.1\n"
   "Copyright 2004 Free Software Foundation, Inc.\n"
   "GDB is free software, covered by the GNU General Public License, and you are\n"
   "welcome to change it and/or distribute copies of it under certain conditions.\n"
   "Type \"show copying\" to see the conditions.\n"
   "There is absolutely no warranty for GDB.  Type \"show warranty\" for details.\n"
   "This GDB was configured as \"--host= --target=athlon-pld-linux\".\n"
   "(gdb) "
2004.07.14 03:27:05
-> "pwd\n"
2004.07.14 03:27:05
<- "Working directory /home/zbyszek/programs/shocks.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info breakpoints\n"
2004.07.14 03:27:05
<- "No breakpoints or watchpoints.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "show history filename\n"
2004.07.14 03:27:05
<- "The filename in which to record the command history is \"/home/zbyszek/programs/shocks/.gdb_history\".\n"
   "(gdb) "
2004.07.14 03:27:05
-> "show history size\n"
2004.07.14 03:27:05
<- "The size of the command history is 1000.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info sources\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line main\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line MAIN\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info files\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "info line main_\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info program\n"
2004.07.14 03:27:05
<- "The program being debugged is not being run.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line MAIN_\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line main__\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line MAIN__\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line _main\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line _MAIN\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line __main\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line __MAIN\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line my_main\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line the_main\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line vxworks_main\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info line main_vxworks\n"
2004.07.14 03:27:05
<- "No symbol table is loaded.  Use the \"file\" command.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "set confirm off\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "source /tmp/dddYWXFo2\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "info breakpoints\n"
2004.07.14 03:27:05
<- "No breakpoints or watchpoints.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "# reset\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "info breakpoints\n"
2004.07.14 03:27:05
<- "No breakpoints or watchpoints.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "display\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "info display\n"
2004.07.14 03:27:05
<- "There are no auto-display expressions now.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "set environment TERM dumb\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "info source\n"
2004.07.14 03:27:05
<- "No current source file.\n"
   "(gdb) "
2004.07.14 03:27:05
-> "info files\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "info program\n"
2004.07.14 03:27:05
<- "The program being debugged is not being run.\n"
   "(gdb) "
#  Starting GDB...done.
#  Welcome to DDD 3.3.9 "Perky" (i686-pld-linux-gnu)
2004.07.14 03:27:05
-> "help detach\n"
2004.07.14 03:27:05
<- "Detach a process or file previously attached.\n"
   "If a process, it is no longer traced, and it continues its execution.  If\n"
   "you were debugging a file, the file is closed and gdb no longer accesses it.\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "help run\n"
2004.07.14 03:27:05
<- "Start debugged program.  You may specify arguments to give it.\n"
   "Args may include \"*\", or \"[...]\"; they are expanded using \"sh\".\n"
   "Input and output redirection with \">\", \"<\", or \">>\" are also allowed.\n"
   "\n"
   "With no arguments, uses arguments last specified (with \"run\" or \"set args\").\n"
   "To cancel previous arguments and run with no arguments,\n"
   "use \"set args\" without arguments.\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "help step\n"
2004.07.14 03:27:05
<- "Step program until it reaches a different source line.\n"
   "Argument N means do this N times (or till program stops for another reason).\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "help stepi\n"
2004.07.14 03:27:05
<- "Step one instruction exactly.\n"
   "Argument N means do this N times (or till program stops for another reason).\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "help next\n"
2004.07.14 03:27:05
<- "Step program, proceeding through subroutine calls.\n"
   "Like the \"step\" command as long as subroutine calls do not happen;\n"
   "when they do, the call is treated as one instruction.\n"
   "Argument N means do this N times (or till program stops for another reason).\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "help nexti\n"
2004.07.14 03:27:05
<- "Step one instruction, but proceed through subroutine calls.\n"
   "Argument N means do this N times (or till program stops for another reason).\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "help until\n"
2004.07.14 03:27:05
<- "Execute until the program reaches a source line greater than the current\n"
   "or a specified location (same args as break command) within the current frame.\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:05
-> "help finish\n"
2004.07.14 03:27:05
<- "Execute until selected stack frame returns.\n"
   "Upon return, the value returned is printed and put in the value history.\n"
2004.07.14 03:27:05
<- "(gdb) "
2004.07.14 03:27:06
-> "help cont\n"
2004.07.14 03:27:06
<- "Continue program being debugged, after signal or breakpoint.\n"
   "If proceeding from breakpoint, a number N may be used as an argument,\n"
   "which means to set the ignore count of that breakpoint to N - 1 (so that\n"
   "the breakpoint won\'t break until the Nth time it is reached).\n"
2004.07.14 03:27:06
<- "(gdb) "
2004.07.14 03:27:06
-> "help signal\n"
2004.07.14 03:27:06
<- "Continue program giving it signal specified by the argument.\n"
   "An argument of \"0\" means continue program without giving it a signal.\n"
2004.07.14 03:27:06
<- "(gdb) "
2004.07.14 03:27:06
-> "help kill\n"
2004.07.14 03:27:06
<- "Kill execution of program being debugged.\n"
2004.07.14 03:27:06
<- "(gdb) "
2004.07.14 03:27:06
-> "help up\n"
2004.07.14 03:27:06
<- "Select and print stack frame that called this one.\n"
   "An argument says how many frames up to go.\n"
2004.07.14 03:27:06
<- "(gdb) "
2004.07.14 03:27:06
-> "help down\n"
2004.07.14 03:27:06
<- "Select and print stack frame called by this one.\n"
   "An argument says how many frames down to go.\n"
2004.07.14 03:27:06
<- "(gdb) "
#  Filtering files...
#  Filtering files...done.
2004.07.14 03:27:16
-> "info files\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "info program\n"
2004.07.14 03:27:16
<- "The program being debugged is not being run.\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "file /home/zbyszek/programs/shocks/bogobox\n"
2004.07.14 03:27:16
<- "Using host libthread_db library \"/lib/libthread_db.so.1\".\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "info line\n"
2004.07.14 03:27:16
<- "No line number information available.\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "list\n"
2004.07.14 03:27:16
<- "67\tFILE *out1, *in1;\n"
2004.07.14 03:27:16
<- "68\t\n"
2004.07.14 03:27:16
<- "69\tint main(int argc, char **argv){\n"
2004.07.14 03:27:16
<- "70\t\tint i, j, k, l;\n"
2004.07.14 03:27:16
<- "71\t\tchar buff[20];\n"
2004.07.14 03:27:16
<- "72\t\tdouble x;\n"
2004.07.14 03:27:16
<- "73\t\tint SoB;\n"
2004.07.14 03:27:16
<- "74\t\tdouble tmp;\n"
2004.07.14 03:27:16
<- "75\t\n"
2004.07.14 03:27:16
<- "76\t\tprintf(\"Using %d box modes\\n"
   "\", NrOfBoxModes);\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "info line\n"
2004.07.14 03:27:16
<- "Line 76 of \"BogolubovBox.c\" starts at address 0x8048948 <main+16> and ends at 0x804895c <main+36>.\n"
2004.07.14 03:27:16
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:76:1489:beg:0x8048948\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "info breakpoints\n"
2004.07.14 03:27:16
<- "No breakpoints or watchpoints.\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "display\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "info display\n"
2004.07.14 03:27:16
<- "There are no auto-display expressions now.\n"
2004.07.14 03:27:16
<- "(gdb) "
#  Reading file "/home/zbyszek/programs/shocks/src/BogolubovBox.c"...
#  Reading file "/home/zbyszek/programs/shocks/src/BogolubovBox.c"...done.
#  File "/home/zbyszek/programs/shocks/src/BogolubovBox.c" 1381 lines, 44648 characters
#  Setting buttons...
#  Setting buttons...done.
2004.07.14 03:27:16
-> "info source\n"
2004.07.14 03:27:16
<- "Current source file is BogolubovBox.c\n"
2004.07.14 03:27:16
<- "Compilation directory is /home/zbyszek/programs/shocks/src\n"
2004.07.14 03:27:16
<- "Located in /home/zbyszek/programs/shocks/src/BogolubovBox.c\n"
2004.07.14 03:27:16
<- "Contains 1381 lines.\n"
2004.07.14 03:27:16
<- "Source language is c.\n"
2004.07.14 03:27:16
<- "Compiled with DWARF 2 debugging format.\n"
2004.07.14 03:27:16
<- "Does not include preprocessor macro info.\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "info files\n"
2004.07.14 03:27:16
<- "Symbols from \"/home/zbyszek/programs/shocks/bogobox\".\n"
2004.07.14 03:27:16
<- "Local exec file:\n"
2004.07.14 03:27:16
<- "\t`/home/zbyszek/programs/shocks/bogobox\', file type elf32-i386.\n"
2004.07.14 03:27:16
<- "\tEntry point: 0x8048880\n"
2004.07.14 03:27:16
<- "\t0x08048134 - 0x08048147 is .interp\n"
2004.07.14 03:27:16
<- "\t0x08048148 - 0x08048168 is .note.ABI-tag\n"
2004.07.14 03:27:16
<- "\t0x08048168 - 0x08048234 is .hash\n"
2004.07.14 03:27:16
<- "\t0x08048234 - 0x08048434 is .dynsym\n"
2004.07.14 03:27:16
<- "\t0x08048434 - 0x080485a5 is .dynstr\n"
2004.07.14 03:27:16
<- "\t0x080485a6 - 0x080485e6 is .gnu.version\n"
2004.07.14 03:27:16
<- "\t0x080485e8 - 0x08048638 is .gnu.version_r\n"
2004.07.14 03:27:16
<- "\t0x08048638 - 0x08048648 is .rel.dyn\n"
2004.07.14 03:27:16
<- "\t0x08048648 - 0x080486e8 is .rel.plt\n"
2004.07.14 03:27:16
<- "\t0x080486e8 - 0x080486ff is .init\n"
2004.07.14 03:27:16
<- "\t0x08048700 - 0x08048850 is .plt\n"
2004.07.14 03:27:16
<- "\t0x08048880 - 0x0804ef58 is .text\n"
2004.07.14 03:27:16
<- "\t0x0804ef58 - 0x0804ef73 is .fini\n"
2004.07.14 03:27:16
<- "\t0x0804ef80 - 0x0804f7b3 is .rodata\n"
2004.07.14 03:27:16
<- "\t0x0804f7b4 - 0x0804f7b8 is .eh_frame\n"
2004.07.14 03:27:16
<- "\t0x080507b8 - 0x080507c0 is .ctors\n"
2004.07.14 03:27:16
<- "\t0x080507c0 - 0x080507c8 is .dtors\n"
2004.07.14 03:27:16
<- "\t0x080507c8 - 0x080507cc is .jcr\n"
2004.07.14 03:27:16
<- "\t0x080507cc - 0x080508c4 is .dynamic\n"
2004.07.14 03:27:16
<- "\t0x080508c4 - 0x080508c8 is .got\n"
2004.07.14 03:27:16
<- "\t0x080508c8 - 0x08050924 is .got.plt\n"
2004.07.14 03:27:16
<- "\t0x08050924 - 0x08050930 is .data\n"
2004.07.14 03:27:16
<- "\t0x08050930 - 0x08050ab4 is .bss\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "info program\n"
2004.07.14 03:27:16
<- "The program being debugged is not being run.\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:16
-> "output fopen\n"
2004.07.14 03:27:16
<- "No symbol \"fopen\" in current context.\n"
2004.07.14 03:27:16
<- "(gdb) "
2004.07.14 03:27:18
-> "output exceeds\n"
2004.07.14 03:27:18
<- "No symbol \"exceeds\" in current context."
2004.07.14 03:27:18
<- "\n"
2004.07.14 03:27:18
<- "(gdb) "
2004.07.14 03:27:37
-> "break main\n"
2004.07.14 03:27:37
<- "Breakpoint 1 at 0x8048948: file BogolubovBox.c, line 76.\n"
2004.07.14 03:27:37
<- "(gdb) "
2004.07.14 03:27:37
-> "info breakpoints\n"
2004.07.14 03:27:37
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:27:37
<- "1   breakpoint     keep y   0x08048948 in main at BogolubovBox.c:76\n"
2004.07.14 03:27:37
<- "(gdb) "
2004.07.14 03:27:37
-> "output readParameters\n"
2004.07.14 03:27:37
<- "{void (char *)} 0x8048bab <readParameters>"
2004.07.14 03:27:37
<- "(gdb) "
#  : waiting until GDB gets ready
2004.07.14 03:27:42
-> "break BogolubovBox.c:91\n"
#  Delete: waiting until GDB gets ready
2004.07.14 03:27:42
<- "Breakpoint 2 at 0x8048a05: file BogolubovBox.c, line 91.\n"
   "(gdb) "
2004.07.14 03:27:42
-> "info breakpoints\n"
2004.07.14 03:27:42
<- "Num Type           Disp Enb Address    What\n"
   "1   breakpoint     keep y   0x08048948 in main at BogolubovBox.c:76\n"
   "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
   "(gdb) "
2004.07.14 03:27:42
-> "delete 1\n"
2004.07.14 03:27:42
<- "(gdb) "
2004.07.14 03:27:42
-> "info breakpoints\n"
2004.07.14 03:27:42
<- "Num Type           Disp Enb Address    What\n"
   "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:27:42
<- "(gdb) "
2004.07.14 03:27:44
-> "set environment TERM dumb\n"
2004.07.14 03:27:44
<- "(gdb) "
2004.07.14 03:27:44
-> "run\n"
2004.07.14 03:27:45
<- "Using 101 box modes\n"
2004.07.14 03:27:45
<- "\n"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "E"
2004.07.14 03:27:45
<- "l"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- "c"
2004.07.14 03:27:45
<- "t"
2004.07.14 03:27:45
<- "r"
2004.07.14 03:27:45
<- "i"
2004.07.14 03:27:45
<- "c"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "F"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- "n"
2004.07.14 03:27:45
<- "c"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "2"
2004.07.14 03:27:45
<- "."
2004.07.14 03:27:45
<- "2"
2004.07.14 03:27:45
<- "."
2004.07.14 03:27:45
<- "0"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "C"
2004.07.14 03:27:45
<- "o"
2004.07.14 03:27:45
<- "p"
2004.07.14 03:27:45
<- "y"
2004.07.14 03:27:45
<- "r"
2004.07.14 03:27:45
<- "i"
2004.07.14 03:27:45
<- "g"
2004.07.14 03:27:45
<- "h"
2004.07.14 03:27:45
<- "t"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "("
2004.07.14 03:27:45
<- "C"
2004.07.14 03:27:45
<- ")"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "1"
2004.07.14 03:27:45
<- "9"
2004.07.14 03:27:45
<- "8"
2004.07.14 03:27:45
<- "7"
2004.07.14 03:27:45
<- "-"
2004.07.14 03:27:45
<- "1"
2004.07.14 03:27:45
<- "9"
2004.07.14 03:27:45
<- "9"
2004.07.14 03:27:45
<- "9"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "B"
2004.07.14 03:27:45
<- "r"
2004.07.14 03:27:45
<- "u"
2004.07.14 03:27:45
<- "c"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "P"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- "r"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- "n"
2004.07.14 03:27:45
<- "s"
2004.07.14 03:27:45
<- " "
2004.07.14 03:27:45
<- "<"
2004.07.14 03:27:45
<- "b"
2004.07.14 03:27:45
<- "r"
2004.07.14 03:27:45
<- "u"
2004.07.14 03:27:45
<- "c"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- "@"
2004.07.14 03:27:45
<- "p"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- "r"
2004.07.14 03:27:45
<- "e"
2004.07.14 03:27:45
<- "n"
2004.07.14 03:27:45
<- "s"
2004.07.14 03:27:45
<- "."
2004.07.14 03:27:45
<- "c"
2004.07.14 03:27:45
<- "o"
2004.07.14 03:27:45
<- "m"
2004.07.14 03:27:45
<- ">"
2004.07.14 03:27:45
<- "\n"
2004.07.14 03:27:45
<- "\n"
   "Breakpoint 2, main (argc=1, argv=0xbffffbe4) at BogolubovBox.c:91\n"
   "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:91:1829:beg:0x8048a05\n"
2004.07.14 03:27:45
<- "(gdb) "
2004.07.14 03:27:45
-> "info breakpoints\n"
2004.07.14 03:27:45
<- "Num Type           Disp Enb Address    What\n"
   "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
   "\tbreakpoint already hit 1 time\n"
2004.07.14 03:27:45
<- "(gdb) "
2004.07.14 03:28:00
-> "stepi\n"
2004.07.14 03:28:00
<- "allocateAndFillUpArrays () at BogolubovBox.c:271\n"
2004.07.14 03:28:00
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:271:5278:beg:0x8048c3d\n"
2004.07.14 03:28:00
<- "(gdb) "
2004.07.14 03:28:00
-> "info breakpoints\n"
2004.07.14 03:28:00
<- "Num Type           Disp Enb Address    What\n"
   "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
   "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:00
<- "(gdb) "
#  Reading file "/home/zbyszek/programs/shocks/BogolubovBox.c"...
#  Reading file "/home/zbyszek/programs/shocks/BogolubovBox.c"...No such file or directory.
2004.07.14 03:28:00
-> "info source\n"
2004.07.14 03:28:00
<- "Current source file is BogolubovBox.c\n"
2004.07.14 03:28:00
<- "Compilation directory is /home/zbyszek/programs/shocks/src\n"
2004.07.14 03:28:00
<- "Located in /home/zbyszek/programs/shocks/src/BogolubovBox.c\n"
2004.07.14 03:28:00
<- "Contains 1381 lines.\n"
2004.07.14 03:28:00
<- "Source language is c.\n"
2004.07.14 03:28:00
<- "Compiled with DWARF 2 debugging format.\n"
2004.07.14 03:28:00
<- "Does not include preprocessor macro info.\n"
2004.07.14 03:28:00
<- "(gdb) "
#  Reading file "BogolubovBox.c" from GDB...
2004.07.14 03:28:00
-> "list BogolubovBox.c:1,1000000\n"
2004.07.14 03:28:00
<- "1\t#include<stdio.h>\n"
2004.07.14 03:28:00
<- "2\t#include<stdlib.h>\n"
2004.07.14 03:28:00
<- "3\t#include<math.h>\n"
2004.07.14 03:28:00
<- "4\t\n"
2004.07.14 03:28:00
<- "5\t#define NrOfBoxModes 101\n"
2004.07.14 03:28:00
<- "6\t\n"
2004.07.14 03:28:00
<- "7\tvoid readParameters(char *fileName);\n"
2004.07.14 03:28:00
<- "8\tvoid allocateAndFillUpArrays(void);\n"
2004.07.14 03:28:00
<- "9\tvoid writeState(char *fileName, double *phi);\n"
2004.07.14 03:28:00
<- "10\tvoid fillHamiltonian(double *H);\n"
2004.07.14 03:28:00
<- "11\tint diagonalize(void);\n"
2004.07.14 03:28:00
<- "12\tint diagonalize2(void);\n"
2004.07.14 03:28:00
<- "13\tvoid outVector(int nr);\n"
2004.07.14 03:28:00
<- "14\tvoid outSortVector(int nr);\n"
2004.07.14 03:28:00
<- "15\tvoid overlap(char *fileName);\n"
2004.07.14 03:28:00
<- "16\tvoid projectUV(void);\n"
2004.07.14 03:28:00
<- "17\tvoid sortUV(void);\n"
2004.07.14 03:28:00
<- "18\tvoid checkoutOrt(void);\n"
2004.07.14 03:28:00
<- "19\textern void indexx(unsigned long n, double arr[], unsigned long indx[]);\n"
2004.07.14 03:28:00
<- "20\tdouble normSign(int nr);\n"
2004.07.14 03:28:00
<- "21\tvoid buildMatricesUV(void);\n"
2004.07.14 03:28:00
<- "22\tvoid buildMatricesUVOrthoBasis(void);\n"
2004.07.14 03:28:00
<- "23\tvoid invertU(void);\n"
2004.07.14 03:28:00
<- "24\tvoid invertUOrthoBasis(void);\n"
2004.07.14 03:28:00
<- "25\textern void gaussj(double *a, int n, double *b, int m);\n"
2004.07.14 03:28:00
<- "26\tvoid createZM(void);\n"
2004.07.14 03:28:00
<- "27\tint findBogoModes(void);\n"
2004.07.14 03:28:00
<- "28\tvoid outBogoModes(void);\n"
2004.07.14 03:28:00
<- "29\tvoid testH(void);\n"
2004.07.14 03:28:00
<- "30\tvoid fillBogoModesX(void);\n"
2004.07.14 03:28:00
<- "31\tvoid outputs(void);\n"
2004.07.14 03:28:00
<- "32\tvoid fillBoxModes(void);\n"
2004.07.14 03:28:00
<- "33\tvoid outUVX(void);\n"
2004.07.14 03:28:00
<- "34\tint family(int nr);\n"
2004.07.14 03:28:00
<- "35\tvoid checkForSymmetry(void);\n"
2004.07.14 03:28:00
<- "36\tvoid reorganizeEigenvectors(void); //becouse of lapack convetions\n"
2004.07.14 03:28:00
<- "37\tvoid checkDum60(void);\n"
2004.07.14 03:28:00
<- "38\t\n"
2004.07.14 03:28:00
<- "39\tint N;\n"
2004.07.14 03:28:00
<- "40\tdouble xL;\n"
2004.07.14 03:28:00
<- "41\tdouble xR;\n"
2004.07.14 03:28:00
<- "42\tint M;\n"
2004.07.14 03:28:00
<- "43\tdouble Planck;\n"
2004.07.14 03:28:00
<- "44\tdouble U0;\n"
2004.07.14 03:28:00
<- "45\tdouble s;\n"
2004.07.14 03:28:00
<- "46\tint NoP;\n"
2004.07.14 03:28:00
<- "47\tdouble g;\n"
2004.07.14 03:28:00
<- "48\tdouble acc;\n"
2004.07.14 03:28:00
<- "49\t\n"
2004.07.14 03:28:00
<- "50\t\n"
2004.07.14 03:28:00
<- "51\tint SizeOfBase;\n"
2004.07.14 03:28:00
<- "52\tdouble *phi, *H;\n"
2004.07.14 03:28:00
<- "53\tdouble *V;\n"
2004.07.14 03:28:00
<- "54\tdouble mu;\n"
2004.07.14 03:28:00
<- "55\tdouble dx, dxdx;\n"
2004.07.14 03:28:00
<- "56\t\n"
2004.07.14 03:28:00
<- "57\tdouble  *wr, *wi, *z;\n"
2004.07.14 03:28:00
<- "58\tdouble *uv;\n"
2004.07.14 03:28:00
<- "59\tdouble *UM, *VM, *ZM;\n"
2004.07.14 03:28:00
<- "60\tunsigned long *indx; \n"
2004.07.14 03:28:00
<- "61\tdouble *BoxModes;\n"
2004.07.14 03:28:00
<- "62\tdouble *Lambda, *BogoModes, *BogoModesX;\n"
2004.07.14 03:28:00
<- "63\tint UVSize;\n"
2004.07.14 03:28:00
<- "64\tint NoM;\n"
2004.07.14 03:28:00
<- "65\tdouble *boxPhi;\n"
2004.07.14 03:28:00
<- "66\t\n"
   "67\tFILE *out1, *in1;\n"
   "68\t\n"
   "69\tint main(int argc, char **argv){\n"
   "70\t\tint i, j, k, l;\n"
   "71\t\tchar buff[20];\n"
   "72\t\tdouble x;\n"
   "73\t\tint SoB;\n"
   "74\t\tdouble tmp;\n"
   "75\t\n"
   "76\t\tprintf(\"Using %d box modes\\n"
   "\", NrOfBoxModes);\n"
   "77\t\treadParameters(\"param.txt\");\n"
   "78\t\tdx = (xR-xL)/N; //bo periodyczne warunki brzegowe.\n"
   "79\t\tdxdx = dx*dx;\n"
   "80\t\tNoM = NrOfBoxModes;\n"
   "81\t\tUVSize = floor((double)NoM*0.9);\n"
   "82\t//\tUVSize = 4;\n"
   "83\t\tif (UVSize>NoM){\n"
   "84\t\t\tfprintf(stderr,\"Size of UV matrix exceeds the size of grid base!\\n"
   "\");\n"
   "85\t\t\texit(1);\n"
   "86\t\t}\n"
   "87\t\n"
   "88\t\n"
   "89\t\tif (argc == 1){\t\n"
   "90\t\t\n"
   "91\t\t\tallocateAndFillUpArrays();\n"
   "92\t\t\tSoB = 2*SizeOfBase;\n"
   "93\t\n"
   "94\t\t\tdiagonalize();\n"
   "95\t\t/*\tout1 = fopen(\"Freq.dat\", \"wt\");\n"
   "96\t\t\tfor(i=0; i<2*NoM; ++i){\n"
   "97\t\t\t\tfprintf(out1, \"%d\\t%.16lf\\t%.16lf\\n"
   "\", i,wr[i], wi[i]);\n"
   "98\t\t\t}\n"
   "99\t\t\tfclose(out1);\n"
   "100\t\t\tcheckDum60();\n"
   "101\t\t\texit(1);\n"
   "102\t\t*/\tsortUV();\n"
   "103\t//\t\toutUVX();\n"
   "104\t\t\tout1 = fopen(\"Freq.dat\", \"wt\");\n"
   "105\t\t\tfor(i=0; i<SoB; ++i){\n"
   "106\t\t\t\tj = family(indx[i]-1);\n"
   "107\t\t\t\tif(1){\n"
   "108\t\t\t\t\tfprintf(out1, \"%d\\t%.16lf\\t%.16lf\\t%d\\n"
   "\", i,wr[indx[i]-1], wi[indx[i]-1], j);\n"
   "109\t\t\t\t}\n"
   "110\t\t\t}\n"
   "111\t\t\tfclose(out1);\n"
   "112\t\t\texit(0);\n"
   "113\t\n"
   "114\t\t\tout1 = fopen(\"freq.dat\", \"wt\");\n"
   "115\t\t\ti=0;\n"
   "116\t\t\tk=0;\n"
   "117\t\t\twhile (i<UVSize && k<(2*SizeOfBase*9)/10-4) {\n"
   "118\t\t\t\tl=k;\n"
   "119\t\t\t\ttmp = normSign(indx[l]-1); // also does projection on ortho. to phi.\n"
   "120\t\t\t\tif(!tmp){\n"
   "121\t\t\t\t\tk+=4;\n"
   "122\t\t\t\t\tcontinue;\n"
   "123\t\t\t\t}\n"
   "124\t\t\t\tif (tmp > 0){\n"
   "125\t\t\t\t\tfor (j=0; j<SoB; ++j) {\n"
   "126\t\t\t\t\t\tuv[i*SoB+j] = z[(indx[l]-1)*SoB+j];\n"
   "127\t\t\t\t\t}\n"
   "128\t\t\t\t}\n"
   "129\t\t\t\telse{\n"
   "130\t\t\t\t\tl=k+1;\n"
   "131\t\t\t\t\ttmp = normSign(indx[l]-1); \n"
   "132\t\t\t\t\tif (tmp > 0){\n"
   "133\t\t\t\t\t\tfor (j=0; j<SoB; ++j) {\n"
   "134\t\t\t\t\t\t\tuv[i*SoB+j] = z[(indx[l]-1)*SoB+j];\n"
   "135\t\t\t\t\t\t}\n"
   "136\t\t\t\t\t}\n"
   "137\t\t\t\t/*\telse{\n"
   "138\t\t\t\t\t\tfprintf(stderr, \"Error in uv selection!\\n"
   "\");\n"
   "139\t\t\t\t\t\texit(1);\n"
   "140\t\t\t\t\t}\t*/\n"
   "141\t\t\t\t\telse{\n"
   "142\t\t\t\t\t\tl=k+2;\n"
   "143\t\t\t\t\t\ttmp = normSign(indx[l]-1); \n"
   "144\t\t\t\t\t\tif (tmp > 0){\n"
   "145\t\t\t\t\t\t\tfor (j=0; j<SoB; ++j) {\n"
   "146\t\t\t\t\t\t\t\tuv[i*SoB+j] = z[(indx[l]-1)*SoB+j];\n"
   "147\t\t\t\t\t\t\t}\n"
   "148\t\t\t\t\t\t}\n"
   "149\t\t\t\t\t\telse{\n"
   "150\t\t\t\t\t\t\tl=k+3;\n"
   "151\t\t\t\t\t\t\ttmp = normSign(indx[l]-1); \n"
   "152\t\t\t\t\t\t\tif (tmp > 0){\n"
   "153\t\t\t\t\t\t\t\tfor (j=0; j<SoB; ++j) {\n"
   "154\t\t\t\t\t\t\t\t\tuv[i*SoB+j] = z[(indx[l]-1)*SoB+j];\n"
   "155\t\t\t\t\t\t\t\t}\n"
   "156\t\t\t\t\t\t\t}\n"
   "157\t\t\t\t\t\t\telse{\n"
   "158\t\t\t\t\t\t\t\tfprintf(stderr, \"Error in uv selection!\\n"
   "\");\n"
   "159\t\t\t\t\t\t\t\texit(1);\n"
   "160\t\t\t\t\t\t\t}\t\n"
   "161\t\t\t\t\t\t}\n"
   "162\t\t\t\t\t}\n"
   "163\t\t\t\t}\n"
   "164\t\t\t\tfprintf(out1, \"%d\\t%.16lf\\t%.16lf\\n"
   "\", l, wr[(indx[l]-1)], wi[(indx[l]-1)]);\n"
   "165\t\t\t\t++i;\n"
   "166\t\t\t\tk+=4;\n"
   "167\t\t\t}\n"
   "168\t\t\tfclose(out1);\n"
   "169\t\n"
   "170\t\t\tif(i<UVSize){\n"
   "171\t\t\t\tUVSize=i;\n"
   "172\t\t\t\tprintf(\"Actual UVSize=%d\\n"
   "\", i);\n"
   "173\t\t\t}\n"
   "174\t\t\tif(UVSize<M)  M = UVSize;\n"
   "175\t\t//\toutUVX();\n"
   "176\t\t/*\t\n"
   "177\t\t\tfor(i=0; i<20; i++) buff[i]=0;\n"
   "178\t\t\tfor (k=0; k<4; ++k) {\n"
   "179\t\t\t\tsprintf(buff,\"zero%d.dat\", k);\n"
   "180\t\t\t\tout1=fopen(buff, \"wt\");\n"
   "181\t\t\t\tfor(i=0; i<NoM; i++){\n"
   "182\t\t\t\t\tl=indx[k]-1;\n"
   "183\t\t\t\t\tfprintf(out1, \"%d\\t%lf\\t%lf\\t%lf\\t%lf\\n"
   "\", i, z[l*SoB+i], z[l*SoB+i+2*NoM], z[l*SoB+i+NoM], z[l*SoB+i+3*NoM]);\n"
   "184\t\t\t\t}\n"
   "185\t\t\t\tfclose(out1);\n"
   "186\t\t\t}\n"
   "187\t\t*/\t\n"
   "188\t\t\tfree(z);\n"
   "189\t\t}\n"
   "190\t\telse if (argc == 2) {\n"
   "191\t\t\tallocateAndFillUpArrays();\n"
   "192\t\t\tSoB = SizeOfBase;\n"
   "193\t\n"
   "194\t\t\tin1 = fopen(argv[1], \"rb\");\n"
   "195\t\t\tif (fread(uv, 8, SoB*UVSize, in1) != SoB*UVSize) {\n"
   "196\t\t\t\tfprintf(stderr, \"bogo: Error reading uv! \\n"
   "\");\n"
   "197\t\t\t\texit(1);\n"
   "198\t\t\t}\n"
   "199\t\t\tfclose(in1);\n"
   "200\t\n"
   "201\t\t}\n"
   "202\t\n"
   "203\t\tcheckDum60();\n"
   "204\t\tcheckoutOrt();\n"
   "205\t\texit(0);\n"
   "206\t\tprojectUV();\n"
   "207\t//\ttestH();\n"
   "208\t//\tbuildMatricesUV();\n"
   "209\t\tbuildMatricesUVOrthoBasis();\n"
   "210\t//\tinvertU();\n"
   "211\t\tinvertUOrthoBasis();\n"
   "212\t\texit(0);\n"
   "213\t\tcreateZM();\n"
   "214\t\tcheckForSymmetry();\n"
   "215\t//\tfindBogoModes();\n"
   "216\t\tdiagonalize2();\n"
   "217\t\tfillBogoModesX();\n"
   "218\t\n"
   "219\t\n"
   "220\t\toutputs();\n"
   "221\t\n"
   "222\t\tfree(H);\n"
   "223\t\tfree(phi);\n"
   "224\t\tfree(V);\n"
   "225\t\tfree(indx);\n"
   "226\t\tfree(uv);\n"
   "227\t\tfree(UM);\n"
   "228\t\tfree(VM);\n"
   "229\t\tfree(ZM);\n"
   "230\t\tfree(Lambda);\n"
   "231\t\tfree(BogoModes);\n"
   "232\t\n"
   "233\t\treturn 0;\n"
   "234\t}\n"
   "235\t\n"
   "236\t\n"
   "237\t\n"
   "238\t\n"
   "239\t\n"
   "240\tvoid readParameters(char *fileName){\n"
   "241\t\t\n"
   "242\t\tin1 = fopen(fileName, \"rt\");\n"
   "243\t\t\n"
   "244\t\tfscanf(in1, \"N=%d\\t%*s\\n"
   "xL=%lf\\t%*s\\n"
   "xR=%lf\\t%*s\\n"
   "M=%d\\t%*s\\n"
   "\"\n"
   "245\t\t\t\"Planck=%lf\\t%*s\\n"
   "aL=%*f\\t%*s\\n"
   "aR=%*f\\t%*s\\n"
   "aN=%*d\\t%*s\\n"
   "\"\n"
   "246\t\t\t\t\"U0=%lf\\t%*s\\n"
   "s=%lf\\t%*s\\n"
   "NoP=%d\\t%*s\\n"
   "g=%lf\\t%*s\\n"
   "\"\n"
   "247\t\t\t\t\"speed=%*f"
2004.07.14 03:28:00
<- "\\t%*s\\n"
   "acc=%lf\\t%*s\\n"
   "atA=%*f\\t%*s\\n"
   "PROG=%*f\\t%*s\\n"
   "\",\n"
2004.07.14 03:28:00
<- "248\t\t\t\t&N, &xL, &xR, &M, &Planck, &U0, &s, &NoP, &g, &acc);\n"
2004.07.14 03:28:00
<- "249\t\t\n"
2004.07.14 03:28:00
<- "250\t\tfclose(in1);\n"
2004.07.14 03:28:00
<- "251\t}\n"
2004.07.14 03:28:00
<- "252\t\n"
2004.07.14 03:28:00
<- "253\t\n"
2004.07.14 03:28:00
<- "254\t/*\n"
2004.07.14 03:28:00
<- "255\tvoid writeState(char *fileName, double *phi){\n"
2004.07.14 03:28:00
<- "256\t    int i,j;\n"
2004.07.14 03:28:00
<- "257\t\n"
2004.07.14 03:28:00
<- "258\t    out1=fopen(fileName, \"wt\");\n"
2004.07.14 03:28:00
<- "259\t    for(i=0; i<SizeOfBase; i++){ \n"
2004.07.14 03:28:00
<- "260\t\tfprintf(out1, \"%d\\t%lf\\t%lf\\t%lf\\t|\", i, phi[i], phi[i+SizeOfBase],\n"
2004.07.14 03:28:00
<- "261\t\t\tphi[i]*phi[i]+phi[i+SizeOfBase]*phi[i+SizeOfBase]);\n"
2004.07.14 03:28:00
<- "262\t    \tfor(j=0; j<M-1; j++) fprintf(out1,\"%d,\", Base[i][j]);\n"
2004.07.14 03:28:00
<- "263\t\tfprintf(out1, \"%d>\\n"
   "\", Base[i][j]);\n"
2004.07.14 03:28:00
<- "264\t    }\n"
2004.07.14 03:28:00
<- "265\t    fclose(out1);\n"
2004.07.14 03:28:00
<- "266\t\n"
2004.07.14 03:28:00
<- "267\t}\n"
2004.07.14 03:28:00
<- "268\t*/\n"
2004.07.14 03:28:00
<- "269\t\n"
2004.07.14 03:28:00
<- "270\t\n"
2004.07.14 03:28:00
<- "271\tvoid allocateAndFillUpArrays(void){\n"
2004.07.14 03:28:00
<- "272\t\tint i, j, k, l;\n"
2004.07.14 03:28:00
<- "273\t\tdouble x;\n"
2004.07.14 03:28:00
<- "274\t\tdouble mem;\n"
2004.07.14 03:28:00
<- "275\t\tdouble sdx;\n"
2004.07.14 03:28:00
<- "276\t\tdouble tmp;\n"
2004.07.14 03:28:00
<- "277\t\tdouble tmpr, tmpi;\n"
2004.07.14 03:28:00
<- "278\t\n"
2004.07.14 03:28:00
<- "279\t\tSizeOfBase = 2*NoM;\n"
2004.07.14 03:28:00
<- "280\t\tphi = (double *)malloc((size_t)2*N*sizeof(double));\n"
2004.07.14 03:28:00
<- "281\t\t\n"
2004.07.14 03:28:00
<- "282\t\tin1=fopen(\"gp.bin\", \"rb\");\n"
2004.07.14 03:28:00
<- "283\t\tif(fread(phi, 8, 2*N, in1)!=2*N) {\n"
2004.07.14 03:28:00
<- "284\t\t\t    fprintf(stderr, \"Reading gp.bin error!\\n"
   "\");\n"
2004.07.14 03:28:00
<- "285\t\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "286\t\t}\n"
2004.07.14 03:28:00
<- "287\t\tif(fread(&mu, 8, 1, in1)!=1) {\n"
2004.07.14 03:28:00
<- "288\t\t\t    fprintf(stderr, \"Reading gp.bin error!\\n"
   "\");\n"
2004.07.14 03:28:00
<- "289\t\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "290\t\t}\n"
2004.07.14 03:28:00
<- "291\t\tif(fread(&tmp, 8, 1, in1)!=1) {\n"
2004.07.14 03:28:00
<- "292\t\t\t    fprintf(stderr, \"Reading gp.bin error!\\n"
   "\");\n"
2004.07.14 03:28:00
<- "293\t\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "294\t\t}\n"
2004.07.14 03:28:00
<- "295\t\tfclose(in1);\n"
2004.07.14 03:28:00
<- "296\t\n"
2004.07.14 03:28:00
<- "297\t\tprintf(\"GP energy = %lf\\n"
   "\", mu);\n"
2004.07.14 03:28:00
<- "298\t\n"
2004.07.14 03:28:00
<- "299\t\ttmp=0.0;\n"
2004.07.14 03:28:00
<- "300\t\tfor(i=0; i<2*N; ++i){\n"
2004.07.14 03:28:00
<- "301\t\t\ttmp += phi[i]*phi[i];\n"
2004.07.14 03:28:00
<- "302\t\t}\n"
2004.07.14 03:28:00
<- "303\t\ttmp = 1.0/sqrt(tmp*dx);\n"
2004.07.14 03:28:00
<- "304\t\tfor(i=0; i<2*N; ++i){\n"
2004.07.14 03:28:00
<- "305\t\t\tphi[i]*=tmp;\n"
2004.07.14 03:28:00
<- "306\t\t}\n"
2004.07.14 03:28:00
<- "307\t\n"
2004.07.14 03:28:00
<- "308\t\n"
2004.07.14 03:28:00
<- "309\t\tBoxModes = (double *)malloc((size_t)N*NoM*sizeof(double));\n"
2004.07.14 03:28:00
<- "310\t\tfillBoxModes();\n"
2004.07.14 03:28:00
<- "311\t\n"
2004.07.14 03:28:00
<- "312\t\tout1 = fopen(\"boxMod.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "313\t\tfor(i=0; i<N; ++i){\n"
2004.07.14 03:28:00
<- "314\t\t\tx=xL+i*dx;\n"
2004.07.14 03:28:00
<- "315\t\t\tfprintf(out1, \"%lf\\t%lf\\n"
   "\", x, BoxModes[1*N+i]);\n"
2004.07.14 03:28:00
<- "316\t\t}\n"
2004.07.14 03:28:00
<- "317\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "318\t\t\n"
2004.07.14 03:28:00
<- "319\t/*\tin1=fopen(\"modeVectors.bin\", \"rb\");\n"
2004.07.14 03:28:00
<- "320\t\tif(fread(BoxModes, 8, N*NrOfBoxModes, in1)!=N*NrOfBoxModes) {\n"
2004.07.14 03:28:00
<- "321\t\t\t    fprintf(stderr, \"Reading modeVectors.bin error!\\n"
   "\");\n"
2004.07.14 03:28:00
<- "322\t\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "323\t\t}\n"
2004.07.14 03:28:00
<- "324\t\tfclose(in1);\n"
2004.07.14 03:28:00
<- "325\t*/\n"
2004.07.14 03:28:00
<- "326\t\t\n"
2004.07.14 03:28:00
<- "327\t\tsdx = sqrt(dx);\n"
2004.07.14 03:28:00
<- "328\t\tout1=fopen(\"gb.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "329\t\tboxPhi = (double *)malloc((size_t)NoM*2*sizeof(double));\n"
2004.07.14 03:28:00
<- "330\t\tfor(i=0; i<NoM; i++){\n"
2004.07.14 03:28:00
<- "331\t\t\ttmpr = 0.0;\n"
2004.07.14 03:28:00
<- "332\t\t\ttmpi = 0.0;\n"
2004.07.14 03:28:00
<- "333\t\t\tfor(j=0; j<N; j++){\n"
2004.07.14 03:28:00
<- "334\t\t\t\ttmpr += BoxModes[i*N+j]*phi[j];\n"
2004.07.14 03:28:00
<- "335\t\t\t\ttmpi += BoxModes[i*N+j]*phi[j+N];\n"
2004.07.14 03:28:00
<- "336\t\t\t}\n"
2004.07.14 03:28:00
<- "337\t\t\tboxPhi[i] = tmpr*sdx;\n"
2004.07.14 03:28:00
<- "338\t\t\tboxPhi[i+NoM] = tmpi*sdx;\n"
2004.07.14 03:28:00
<- "339\t\t\tfprintf(out1, \"%d\\t%lf\\t%lf\\t%lf\\n"
   "\", i, tmpr*sdx, tmpi*sdx, (tmpr*tmpr+tmpi*tmpi)*dx);\n"
2004.07.14 03:28:00
<- "340\t\t}\n"
2004.07.14 03:28:00
<- "341\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "342\t\t\n"
2004.07.14 03:28:00
<- "343\t\ttmpr = 0.0;\n"
2004.07.14 03:28:00
<- "344\t\tfor(i=0; i<NoM; ++i){\n"
2004.07.14 03:28:00
<- "345\t\t\ttmpr += boxPhi[i]*boxPhi[i]+boxPhi[i+NoM]*boxPhi[i+NoM];\n"
2004.07.14 03:28:00
<- "346\t\t}\n"
2004.07.14 03:28:00
<- "347\t\tprintf(\"Norma = %.10lf\\n"
   "\", tmpr);\n"
2004.07.14 03:28:00
<- "348\t\ttmpr = 1.0/sqrt(tmpr);\n"
2004.07.14 03:28:00
<- "349\t\tfor(i=0; i<NoM; ++i){\n"
2004.07.14 03:28:00
<- "350\t\t\tboxPhi[i]*=tmpr;\n"
2004.07.14 03:28:00
<- "351\t\t\tboxPhi[i+NoM]*=tmpr;\n"
2004.07.14 03:28:00
<- "352\t\t}\n"
2004.07.14 03:28:00
<- "353\t\n"
2004.07.14 03:28:00
<- "354\t\t\n"
2004.07.14 03:28:00
<- "355\t\tV = (double *)malloc((size_t)N*sizeof(double));\n"
2004.07.14 03:28:00
<- "356\t\tfor (i=0; i<N; ++i){\n"
2004.07.14 03:28:00
<- "357\t\t\tx = xL+i*dx;\n"
2004.07.14 03:28:00
<- "358\t\t\tV[i] = 0.0*x*x;\n"
2004.07.14 03:28:00
<- "359\t\t}\n"
2004.07.14 03:28:00
<- "360\t\t\n"
2004.07.14 03:28:00
<- "361\t\n"
2004.07.14 03:28:00
<- "362\t\tH = (double *)malloc((size_t)4*SizeOfBase*SizeOfBase*sizeof(double));\n"
2004.07.14 03:28:00
<- "363\t\tfor(i=0; i<4*SizeOfBase*SizeOfBase; i++) H[i]=0.0;\n"
2004.07.14 03:28:00
<- "364\t\n"
2004.07.14 03:28:00
<- "365\t\tfillHamiltonian(H);\n"
2004.07.14 03:28:00
<- "366\t\n"
2004.07.14 03:28:00
<- "367\t\tout1 = fopen(\"h11r.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "368\t\tfor(i=0; i<NoM; ++i){\n"
2004.07.14 03:28:00
<- "369\t\t\tfor(j=0; j<NoM; ++j){\n"
2004.07.14 03:28:00
<- "370\t\t\t\tfprintf(out1, \"%.4lf  \", H[i*2*SizeOfBase+j]);\n"
2004.07.14 03:28:00
<- "371\t\t\t}\n"
2004.07.14 03:28:00
<- "372\t\t\tfprintf(out1,\"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "373\t\t}\n"
2004.07.14 03:28:00
<- "374\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "375\t\n"
2004.07.14 03:28:00
<- "376\t\tout1 = fopen(\"h11i.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "377\t\tfor(i=0; i<NoM; ++i){\n"
2004.07.14 03:28:00
<- "378\t\t\tfor(j=0; j<NoM; ++j){\n"
2004.07.14 03:28:00
<- "379\t\t\t\tfprintf(out1, \"%.4lf  \", H[i*2*SizeOfBase+j+2*NoM]);\n"
2004.07.14 03:28:00
<- "380\t\t\t}\n"
2004.07.14 03:28:00
<- "381\t\t\tfprintf(out1,\"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "382\t\t}\n"
2004.07.14 03:28:00
<- "383\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "384\t\n"
2004.07.14 03:28:00
<- "385\t\tout1 = fopen(\"h12r.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "386\t\tfor(i=0; i<NoM; ++i){\n"
2004.07.14 03:28:00
<- "387\t\t\tfor(j=0; j<NoM; ++j){\n"
2004.07.14 03:28:00
<- "388\t\t\t\tfprintf(out1, \"%.4lf  \", H[i*2*SizeOfBase+j+NoM]);\n"
2004.07.14 03:28:00
<- "389\t\t\t}\n"
2004.07.14 03:28:00
<- "390\t\t\tfprintf(out1,\"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "391\t\t}\n"
2004.07.14 03:28:00
<- "392\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "393\t\n"
2004.07.14 03:28:00
<- "394\t\tout1 = fopen(\"h12i.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "395\t\tfor(i=0; i<NoM; ++i){\n"
2004.07.14 03:28:00
<- "396\t\t\tfor(j=0; j<NoM; ++j){\n"
2004.07.14 03:28:00
<- "397\t\t\t\tfprintf(out1, \"%.4lf  \", H[i*2*SizeOfBase+j+3*NoM]);\n"
2004.07.14 03:28:00
<- "398\t\t\t}\n"
2004.07.14 03:28:00
<- "399\t\t\tfprintf(out1,\"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "400\t\t}\n"
2004.07.14 03:28:00
<- "401\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "402\t\n"
2004.07.14 03:28:00
<- "403\t\n"
2004.07.14 03:28:00
<- "404\t\tz = (double *)malloc((size_t)4*SizeOfBase*SizeOfBase*sizeof(double));\n"
2004.07.14 03:28:00
<- "405\t\twr = (double *)malloc((size_t)2*SizeOfBase*sizeof(double));\n"
2004.07.14 03:28:00
<- "406\t\twi = (double *)malloc((size_t)2*SizeOfBase*sizeof(double));\n"
2004.07.14 03:28:00
<- "407\t\t\n"
2004.07.14 03:28:00
<- "408\t\tindx = (unsigned long *)malloc((size_t)2*SizeOfBase*sizeof(unsigned long));\n"
2004.07.14 03:28:00
<- "409\t\tuv = (double *)malloc((size_t)2*SizeOfBase*UVSize*sizeof(double));\n"
2004.07.14 03:28:00
<- "410\t\n"
2004.07.14 03:28:00
<- "411\t\tmem = ((2.0*SizeOfBase +UVSize + 2)*SizeOfBase+N+N*NoM+N+NoM)*8;\n"
2004.07.14 03:28:00
<- "412\t\tmem = mem/1024/1024;\n"
2004.07.14 03:28:00
<- "413\t\tprintf(\"Memory size demand: %.1lfMB\\n"
   "\", mem);\n"
2004.07.14 03:28:00
<- "414\t}\n"
2004.07.14 03:28:00
<- "415\t\n"
2004.07.14 03:28:00
<- "416\tvoid fillHamiltonian(double *H){\n"
2004.07.14 03:28:00
<- "417\t\tint i,j,k, l;\n"
2004.07.14 03:28:00
<- "418\t\tdouble *st1, *st2;\n"
2004.07.14 03:28:00
<- "419\t\tdouble tmpr, tmpi;\n"
2004.07.14 03:28:00
<- "420\t\tdouble c=-Planck*Planck*0.5/dxdx;\n"
2004.07.14 03:28:00
<- "421\t\t\n"
2004.07.14 03:28:00
<- "422\t\tfor(i=0; i<NoM; ++i){\n"
2004.07.14 03:28:00
<- "423\t\t\tst2 = BoxModes+i*N;\n"
2004.07.14 03:28:00
<- "424\t\t\tfor(j=0; j<NoM; ++j){\n"
2004.07.14 03:28:00
<- "425\t\t\t\tst1 = BoxModes+j*N;\n"
2004.07.14 03:28:00
<- "426\t\t\t\ttmpr = 0.0;\n"
2004.07.14 03:28:00
<- "427\t\n"
2004.07.14 03:28:00
<- "428\t\t\t\tk=0;\n"
2004.07.14 03:28:00
<- "429\t\t   \t\ttmpr += st2[k]*(c*(st1[k+1]+st1[N-1])+(-2.0*c+V[k] - mu)*st1[k]);\n"
2004.07.14 03:28:00
<- "430\t\t\t\tfor(k=1; k<N-1; ++k){\n"
2004.07.14 03:28:00
<- "431\t\t   \t\t\ttmpr += st2[k]*(c*(st1[k+1]+st1[k-1])+(-2.0*c+V[k] - mu)*st1[k]);\n"
2004.07.14 03:28:00
<- "432\t\t\t\t}\n"
2004.07.14 03:28:00
<- "433\t\t\t\tk=N-1;\n"
2004.07.14 03:28:00
<- "434\t\t   \t\ttmpr += st2[k]*(c*(st1[0]+st1[k-1])+(-2.0*c+V[k] - mu)*st1[k]);\n"
2004.07.14 03:28:00
<- "435\t\n"
2004.07.14 03:28:00
<- "436\t\t\t\tfor(k=0; k<N; ++k){\n"
2004.07.14 03:28:00
<- "437\t\t   \t\t\ttmpr += st2[k]*2.0*g*(phi[k]*phi[k]+phi[k+N]*phi[k+N])*st1[k];\n"
2004.07.14 03:28:00
<- "438\t\t\t\t}\n"
2004.07.14 03:28:00
<- "439\t\t   \t\tH[i + 2*SizeOfBase*j] = tmpr;\n"
2004.07.14 03:28:00
<- "440\t\t   \t\tH[i+NoM + 2*SizeOfBase*(j+NoM)] = -tmpr;\n"
2004.07.14 03:28:00
<- "441\t\t   \t\tH[i+2*NoM + 2*SizeOfBase*(j+2*NoM)] = tmpr;\n"
2004.07.14 03:28:00
<- "442\t\t   \t\tH[i+3*NoM + 2*SizeOfBase*(j+3*NoM)] = -tmpr;\n"
2004.07.14 03:28:00
<- "443\t\t\t\t\n"
2004.07.14 03:28:00
<- "444\t\t\t\ttmpr = 0.0;\n"
2004.07.14 03:28:00
<- "445\t\t\t\ttmpi = 0.0;\n"
2004.07.14 03:28:00
<- "446\t\t\t\tfor(k=0; k<N; ++k){\n"
2004.07.14 03:28:00
<- "447\t\t   \t\t\ttmpr += st2[k]*g*(phi[k]*phi[k]-phi[k+N]*phi[k+N])*st1[k];\n"
2004.07.14 03:28:00
<- "448\t\t   \t\t\ttmpi += st2[k]*g*2.0*phi[k]*phi[k+N]*st1[k];\n"
2004.07.14 03:28:00
<- "449\t\t\t\t}\n"
2004.07.14 03:28:00
<- "450\t\t   \t\tH[i + 2*SizeOfBase*(NoM + j)] = tmpr;\n"
2004.07.14 03:28:00
<- "451\t\t   \t\tH[i+NoM + 2*SizeOfBase*j] = -tmpr;\n"
2004.07.14 03:28:00
<- "452\t\t   \t\tH[i+2*NoM+ 2*SizeOfBase*(3*NoM + j)] = tmpr;\n"
2004.07.14 03:28:00
<- "453\t\t   \t\tH[i+3*NoM + 2*SizeOfBase*(j+2*NoM)] = -tmpr;\n"
2004.07.14 03:28:00
<- "454\t\t\t\t\n"
2004.07.14 03:28:00
<- "455\t\t   \t\tH[i + 2*SizeOfBase*(3*NoM + j)] = tmpi;\n"
2004.07.14 03:28:00
<- "456\t\t   \t\tH[i+NoM + 2*SizeOfBase*(2*NoM+j)] = tmpi;\n"
2004.07.14 03:28:00
<- "457\t\t   \t\tH[i+2*NoM+ 2*SizeOfBase*(NoM + j)] = -tmpi;\n"
2004.07.14 03:28:00
<- "458\t\t   \t\tH[i+3*NoM + 2*SizeOfBase*j] = -tmpi;\n"
2004.07.14 03:28:00
<- "459\t\t\t}\n"
2004.07.14 03:28:00
<- "460\t\t}\n"
2004.07.14 03:28:00
<- "461\t}\n"
2004.07.14 03:28:00
<- "462\t\n"
2004.07.14 03:28:00
<- "463\t\n"
2004.07.14 03:28:00
<- "464\tint diagonalize(void){\n"
2004.07.14 03:28:00
<- "465\t\tchar balanc=\'N\', jobvl=\'N\';\n"
2004.07.14 03:28:00
<- "466\t\tchar jobvr=\'V\', sense=\'N\';\n"
2004.07.14 03:28:00
<- "467\t\tint il, iu, info, n=2*SizeOfBase; //il>=1 iu=il+SizeOfBase-1\n"
2004.07.14 03:28:00
<- "468\t\tint lda=n, lwork=8*n;\n"
2004.07.14 03:28:00
<- "469\t\tdouble abnrm;\n"
2004.07.14 03:28:00
<- "470\t\tint *iwork;\n"
2004.07.14 03:28:00
<- "471\t\tdouble *scale,*rconde, *rcondv, *work;\n"
2004.07.14 03:28:00
<- "472\t\n"
2004.07.14 03:28:00
<- "473\t\tiwork=(int *)malloc((size_t)2*n*sizeof(int));\n"
2004.07.14 03:28:00
<- "474\t\tscale=(double *)malloc((size_t)n*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "475\t\trconde=(double *)malloc((size_t)n*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "476\t\trcondv=(double *)malloc((size_t)n*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "477\t\twork=(double *)malloc((size_t)lwork*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "478\t\n"
2004.07.14 03:28:00
<- "479\t\tdgeevx_( &balanc, &jobvl, &jobvr, &sense, &n, H, &lda, wr, wi, z, &n, z, &n, &il, &iu, scale, &abnrm, rconde, rcondv, work, &lwork, iwork, &info );\n"
2004.07.14 03:28:00
<- "480\t\n"
2004.07.14 03:28:00
<- "481\t\tif(info){\n"
2004.07.14 03:28:00
<- "482\t\t\tfprintf(stderr,\"diagonalize(): Problem with diagonalization! \\n"
   "\");\n"
2004.07.14 03:28:00
<- "483\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "484\t\t}\n"
2004.07.14 03:28:00
<- "485\t\n"
2004.07.14 03:28:00
<- "486\t\treorganizeEigenvectors();\n"
2004.07.14 03:28:00
<- "487\t\n"
2004.07.14 03:28:00
<- "488\t\tfree(scale);\n"
2004.07.14 03:28:00
<- "489\t\tfree(rconde);\n"
2004.07.14 03:28:00
<- "490\t\tfree(rcondv);\n"
2004.07.14 03:28:00
<- "491\t\tfree(iwork);\n"
2004.07.14 03:28:00
<- "492\t\tfree(work);\n"
2004.07.14 03:28:00
<- "493\t\treturn info;\n"
2004.07.14 03:28:00
<- "494\t}\n"
2004.07.14 03:28:00
<- "495\t\n"
2004.07.14 03:28:00
<- "496\tvoid projectUV(void){\n"
2004.07.14 03:28:00
<- "497\t\tint i,j,k;\t\n"
2004.07.14 03:28:00
<- "498\t\tdouble prour, proui, provr, provi;\n"
2004.07.14 03:28:00
<- "499\t\tdouble *ur, *ui, *vr, *vi;\n"
2004.07.14 03:28:00
<- "500\t\tdouble *br, *bi;\n"
2004.07.14 03:28:00
<- "501\t\tdouble normu, normv, norm;\n"
2004.07.14 03:28:00
<- "502\t\tdouble tmp, sdx;\n"
2004.07.14 03:28:00
<- "503\t\n"
2004.07.14 03:28:00
<- "504\t\tbr = boxPhi;\n"
2004.07.14 03:28:00
<- "505\t\tbi = boxPhi+NoM;\n"
2004.07.14 03:28:00
<- "506\t\n"
2004.07.14 03:28:00
<- "507\t\tout1 = fopen(\"proj.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "508\t\tfor(i=0; i<UVSize; i++){\n"
2004.07.14 03:28:00
<- "509\t\t\tprour = 0.0;\n"
2004.07.14 03:28:00
<- "510\t\t\tproui = 0.0;\n"
2004.07.14 03:28:00
<- "511\t\t\tprovr = 0.0;\n"
2004.07.14 03:28:00
<- "512\t\t\tprovi = 0.0;\n"
2004.07.14 03:28:00
<- "513\t\t\tur = uv+i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "514\t\t\tvr = ur+NoM;\n"
2004.07.14 03:28:00
<- "515\t\t\tui = vr+NoM;\n"
2004.07.14 03:28:00
<- "516\t\t\tvi = ui+NoM;\n"
2004.07.14 03:28:00
<- "517\t\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "518\t\t\t\tprour += ur[j]*br[j]+ui[j]*bi[j];\n"
2004.07.14 03:28:00
<- "519\t\t\t\tproui += ur[j]*bi[j]-ui[j]*br[j];\n"
2004.07.14 03:28:00
<- "520\t\t\t\tprovr += vr[j]*br[j]+vi[j]*bi[j];\n"
2004.07.14 03:28:00
<- "521\t\t\t\tprovi += vr[j]*bi[j]-vi[j]*br[j];\n"
2004.07.14 03:28:00
<- "522\t\t\t}\n"
2004.07.14 03:28:00
<- "523\t\t\tfprintf(out1, \"%lf\\t%lf\\t%lf\\t%lf\\n"
   "\", prour, proui, provr, provi);\n"
2004.07.14 03:28:00
<- "524\t\t}\n"
2004.07.14 03:28:00
<- "525\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "526\t\t\n"
2004.07.14 03:28:00
<- "527\t\t/*\n"
2004.07.14 03:28:00
<- "528\t//\t\tprou *= dx;\n"
2004.07.14 03:28:00
<- "529\t//\t\tprov *= dx;\n"
2004.07.14 03:28:00
<- "530\t\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "531\t\t\t\tu[j] -= prou*boxPhi[j];\n"
2004.07.14 03:28:00
<- "532\t\t\t\tv[j] -= prov*boxPhi[j];\n"
2004.07.14 03:28:00
<- "533\t\t\t}\n"
2004.07.14 03:28:00
<- "534\t\t\tnormu = 0.0;\n"
2004.07.14 03:28:00
<- "535\t\t\tnormv = 0.0;\n"
2004.07.14 03:28:00
<- "536\t\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "537\t\t\t\tnormu += u[j]*u[j];\n"
2004.07.14 03:28:00
<- "538\t\t\t\tnormv += v[j]*v[j];\n"
2004.07.14 03:28:00
<- "539\t\t\t}\n"
2004.07.14 03:28:00
<- "540\t//\t\tnorm = fabs(normu-normv)*dx;\n"
2004.07.14 03:28:00
<- "541\t\t\tnorm = fabs(normu-normv);\n"
2004.07.14 03:28:00
<- "542\t\t\tif(norm < 1e-14) {\n"
2004.07.14 03:28:00
<- "543\t\t\t\tprintf(\"n: %d\\n"
   "\", i);\n"
2004.07.14 03:28:00
<- "544\t\t\t\tcontinue;\n"
2004.07.14 03:28:00
<- "545\t\t\t}\n"
2004.07.14 03:28:00
<- "546\t\t\tnorm = 1.0/sqrt(norm);\n"
2004.07.14 03:28:00
<- "547\t\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "548\t\t\t\tu[j] *= norm;\n"
2004.07.14 03:28:00
<- "549\t\t\t\tv[j] *= norm;\n"
2004.07.14 03:28:00
<- "550\t\t\t}\n"
2004.07.14 03:28:00
<- "551\t\t}*/\n"
2004.07.14 03:28:00
<- "552\t}\n"
2004.07.14 03:28:00
<- "553\t\n"
2004.07.14 03:28:00
<- "554\tvoid sortUV(void){\n"
2004.07.14 03:28:00
<- "555\t\tint i;\n"
2004.07.14 03:28:00
<- "556\t\tdouble *tab;\n"
2004.07.14 03:28:00
<- "557\t\n"
2004.07.14 03:28:00
<- "558\t\ttab = (double *)malloc((size_t)2*SizeOfBase*sizeof(double));\n"
2004.07.14 03:28:00
<- "559\t\n"
2004.07.14 03:28:00
<- "560\t\tfor(i=0; i<2*SizeOfBase; ++i) {\n"
2004.07.14 03:28:00
<- "561\t\t\ttab[i] = wr[i]*wr[i]+wi[i]*wi[i];\n"
2004.07.14 03:28:00
<- "562\t\t}\n"
2004.07.14 03:28:00
<- "563\t\tindexx(2*SizeOfBase, tab-1, indx-1);\n"
2004.07.14 03:28:00
<- "564\t\t\n"
2004.07.14 03:28:00
<- "565\t/*\tfor(i=0; i<2*SizeOfBase; ++i) {\n"
2004.07.14 03:28:00
<- "566\t\t\tprintf(\"%lf\\n"
   "\", wr[indx[i]-1]);\n"
2004.07.14 03:28:00
<- "567\t\t}\n"
2004.07.14 03:28:00
<- "568\t*/\n"
2004.07.14 03:28:00
<- "569\t\tfree(tab);\n"
2004.07.14 03:28:00
<- "570\t}\n"
2004.07.14 03:28:00
<- "571\t\n"
2004.07.14 03:28:00
<- "572\t\n"
2004.07.14 03:28:00
<- "573\tvoid checkoutOrt(void){\n"
2004.07.14 03:28:00
<- "574\t\tint i,j;\t\n"
2004.07.14 03:28:00
<- "575\t\tdouble tmpr, tmpi;\n"
2004.07.14 03:28:00
<- "576\t\tdouble *ur, *vr, *ui, *vi;\n"
2004.07.14 03:28:00
<- "577\t\tdouble *bPr, *bPi;\n"
2004.07.14 03:28:00
<- "578\t\n"
2004.07.14 03:28:00
<- "579\t\n"
2004.07.14 03:28:00
<- "580\t\tbPr = boxPhi;\n"
2004.07.14 03:28:00
<- "581\t\tbPi = bPr + NoM;\n"
2004.07.14 03:28:00
<- "582\t\t\n"
2004.07.14 03:28:00
<- "583\t\tout1 = fopen(\"ortU.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "584\t\tfor(i=0; i<UVSize; ++i) {\n"
2004.07.14 03:28:00
<- "585\t\t\tur = uv+i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "586\t\t\tui = ur + 2*NoM;\n"
2004.07.14 03:28:00
<- "587\t\t\ttmpr = tmpi = 0.0;\n"
2004.07.14 03:28:00
<- "588\t\t\tfor(j=0; j<NoM; ++j) {\n"
2004.07.14 03:28:00
<- "589\t\t\t\ttmpr += bPr[j]*ur[j]+bPi[j]*ui[j];\n"
2004.07.14 03:28:00
<- "590\t\t\t\ttmpi += bPr[j]*ui[j]-bPi[j]*ur[j];\n"
2004.07.14 03:28:00
<- "591\t\t\t}\n"
2004.07.14 03:28:00
<- "592\t\t\tfprintf(out1, \"%d\\t%.16lf\\t%.16lf\\n"
   "\", i, tmpr, tmpi);\n"
2004.07.14 03:28:00
<- "593\t\t}\n"
2004.07.14 03:28:00
<- "594\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "595\t\n"
2004.07.14 03:28:00
<- "596\t\tout1 = fopen(\"ortV.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "597\t\tfor(i=0; i<UVSize; ++i) {\n"
2004.07.14 03:28:00
<- "598\t\t\tvr = uv+i*2*SizeOfBase+NoM;\n"
2004.07.14 03:28:00
<- "599\t\t\tvi = vr + 2*NoM;\n"
2004.07.14 03:28:00
<- "600\t\t\ttmpr = tmpi = 0.0;\n"
2004.07.14 03:28:00
<- "601\t\t\tfor(j=0; j<NoM; ++j) {\n"
2004.07.14 03:28:00
<- "602\t\t\t\ttmpr += bPr[j]*vr[j]-bPi[j]*vi[j];\n"
2004.07.14 03:28:00
<- "603\t\t\t\ttmpi += bPr[j]*vi[j]+bPi[j]*vr[j];\n"
2004.07.14 03:28:00
<- "604\t\t\t}\n"
2004.07.14 03:28:00
<- "605\t\t\tfprintf(out1, \"%d\\t%.16lf\\t%.16lf\\n"
   "\", i, tmpr, tmpi);\n"
2004.07.14 03:28:00
<- "606\t\t}\n"
2004.07.14 03:28:00
<- "607\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "608\t}\n"
2004.07.14 03:28:00
<- "609\t\n"
2004.07.14 03:28:00
<- "610\tdouble normSign(int nr){\n"
2004.07.14 03:28:00
<- "611\t\tint i,j,k;\t\n"
2004.07.14 03:28:00
<- "612\t\tdouble prour, proui, provr, provi;\n"
2004.07.14 03:28:00
<- "613\t\tdouble *ur, *ui, *vr, *vi;\n"
2004.07.14 03:28:00
<- "614\t\tdouble *br, *bi;\n"
2004.07.14 03:28:00
<- "615\t\tdouble normu, normv, norm;\n"
2004.07.14 03:28:00
<- "616\t\tdouble tmp;\n"
2004.07.14 03:28:00
<- "617\t\n"
2004.07.14 03:28:00
<- "618\t\tbr = boxPhi;\n"
2004.07.14 03:28:00
<- "619\t\tbi = boxPhi+NoM;\n"
2004.07.14 03:28:00
<- "620\t\n"
2004.07.14 03:28:00
<- "621\t\tur = z+nr*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "622\t\tvr = ur+NoM;\n"
2004.07.14 03:28:00
<- "623\t\tui = vr+NoM;\n"
2004.07.14 03:28:00
<- "624\t\tvi = ui+NoM;\n"
2004.07.14 03:28:00
<- "625\t\t\n"
2004.07.14 03:28:00
<- "626\t\tprour = 0.0;\n"
2004.07.14 03:28:00
<- "627\t\tproui = 0.0;\n"
2004.07.14 03:28:00
<- "628\t\tprovr = 0.0;\n"
2004.07.14 03:28:00
<- "629\t\tprovi = 0.0;\n"
2004.07.14 03:28:00
<- "630\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "631\t\t\tprour += ur[j]*br[j]+ui[j]*bi[j];\n"
2004.07.14 03:28:00
<- "632\t\t\tproui += ui[j]*br[j]-ur[j]*bi[j];\n"
2004.07.14 03:28:00
<- "633\t\t\tprovr += vr[j]*br[j]-vi[j]*bi[j]; //rzut v na psi^*\n"
2004.07.14 03:28:00
<- "634\t\t\tprovi += vi[j]*br[j]+vr[j]*bi[j];\n"
2004.07.14 03:28:00
<- "635\t\t}\n"
2004.07.14 03:28:00
<- "636\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "637\t\t\tur[j] -= prour*br[j]-proui*bi[j];\n"
2004.07.14 03:28:00
<- "638\t\t\tui[j] -= prour*bi[j]+proui*br[j];\n"
2004.07.14 03:28:00
<- "639\t\t\tvr[j] -= provr*br[j]+provi*bi[j];\n"
2004.07.14 03:28:00
<- "640\t\t\tvi[j] -= -provr*bi[j]+provi*br[j];\n"
2004.07.14 03:28:00
<- "641\t\t}\n"
2004.07.14 03:28:00
<- "642\t\tnormu = 0.0;\n"
2004.07.14 03:28:00
<- "643\t\tnormv = 0.0;\n"
2004.07.14 03:28:00
<- "644\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "645\t\t\tnormu += ur[j]*ur[j]+ui[j]*ui[j];\n"
2004.07.14 03:28:00
<- "646\t\t\tnormv += vr[j]*vr[j]+vi[j]*vi[j];\n"
2004.07.14 03:28:00
<- "647\t\t}\n"
2004.07.14 03:28:00
<- "648\t\ttmp = (normu-normv);\n"
2004.07.14 03:28:00
<- "649\t\tnorm = fabs(tmp);\n"
2004.07.14 03:28:00
<- "650\t\tif(norm < 1e-6) {\n"
2004.07.14 03:28:00
<- "651\t\t//\tprintf(\"n: %d\\n"
   "\", nr);\n"
2004.07.14 03:28:00
<- "652\t\t\treturn 0.0;\t\n"
2004.07.14 03:28:00
<- "653\t\t}\n"
2004.07.14 03:28:00
<- "654\t\tnorm = 1.0/sqrt(norm);\n"
2004.07.14 03:28:00
<- "655\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "656\t\t\tur[j] *= norm;\n"
2004.07.14 03:28:00
<- "657\t\t\tui[j] *= norm;\n"
2004.07.14 03:28:00
<- "658\t\t\tvr[j] *= norm;\n"
2004.07.14 03:28:00
<- "659\t\t\tvi[j] *= norm;\n"
2004.07.14 03:28:00
<- "660\t\t}\n"
2004.07.14 03:28:00
<- "661\t\n"
2004.07.14 03:28:00
<- "662\t\treturn tmp;\n"
2004.07.14 03:28:00
<- "663\t}\n"
2004.07.14 03:28:00
<- "664\t\n"
2004.07.14 03:28:00
<- "665\tvoid buildMatricesUVOrthoBasis(void){\n"
2004.07.14 03:28:00
<- "666\t\tint i, j, k;\t\n"
2004.07.14 03:28:00
<- "667\t\tdouble normur, normui, normuvr, normuvi;\n"
2004.07.14 03:28:00
<- "668\t\tdouble *ur, *ui, *vr, *vi; \n"
2004.07.14 03:28:00
<- "669\t\t\n"
2004.07.14 03:28:00
<- "670\t\tUM = (double *)malloc((size_t)4*UVSize*NoM*sizeof(double));\n"
2004.07.14 03:28:00
<- "671\t\tVM = (double *)malloc((size_t)4*UVSize*NoM*sizeof(double));\n"
2004.07.14 03:28:00
<- "672\t\n"
2004.07.14 03:28:00
<- "673\t\tfor (i=0; i<UVSize; ++i){\n"
2004.07.14 03:28:00
<- "674\t\t\tur = uv + i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "675\t\t\tvr = ur + NoM;\n"
2004.07.14 03:28:00
<- "676\t\t\tui = vr + NoM;\n"
2004.07.14 03:28:00
<- "677\t\t\tvi = ui + NoM;\n"
2004.07.14 03:28:00
<- "678\t\t\tfor (k=0; k<NoM; ++k){\n"
2004.07.14 03:28:00
<- "679\t\t\t\tUM[k*2*UVSize+i] = ur[k];\n"
2004.07.14 03:28:00
<- "680\t\t\t\tUM[(k+NoM)*2*UVSize+i+UVSize] = ur[k];\n"
2004.07.14 03:28:00
<- "681\t\t\t\tUM[(k+NoM)*2*UVSize+i] = ui[k];\n"
2004.07.14 03:28:00
<- "682\t\t\t\tUM[k*2*UVSize+i+UVSize] = -ui[k];\n"
2004.07.14 03:28:00
<- "683\t\t\t/*\t\n"
2004.07.14 03:28:00
<- "684\t\t\t\tUM[i*2*NoM+k] = ur[k];\n"
2004.07.14 03:28:00
<- "685\t\t\t\tUM[(i+UVSize)*2*NoM+k+NoM] = ur[k];\n"
2004.07.14 03:28:00
<- "686\t\t\t\tUM[i*2*NoM+k+NoM] = ui[k];\n"
2004.07.14 03:28:00
<- "687\t\t\t\tUM[(i+UVSize)*2*NoM+k] = -ui[k];\n"
2004.07.14 03:28:00
<- "688\t*/\n"
2004.07.14 03:28:00
<- "689\t\t\t\tVM[i*2*NoM+k] = vr[k];\n"
2004.07.14 03:28:00
<- "690\t\t\t\tVM[(i+UVSize)*2*NoM+k+NoM] = vr[k];\n"
2004.07.14 03:28:00
<- "691\t\t\t\tVM[i*2*NoM+k+NoM] = vi[k];\n"
2004.07.14 03:28:00
<- "692\t\t\t\tVM[(i+UVSize)*2*NoM+k] = -vi[k];\n"
2004.07.14 03:28:00
<- "693\t\n"
2004.07.14 03:28:00
<- "694\t\t\t}\n"
2004.07.14 03:28:00
<- "695\t\t}\n"
2004.07.14 03:28:00
<- "696\t\t\n"
2004.07.14 03:28:00
<- "697\t\tout1 = fopen(\"um.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "698\t\tfor (i=0; i<2*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "699\t\t\tfor (j=0; j<2*NoM; ++j){\n"
2004.07.14 03:28:00
<- "700\t\t\t\tfprintf(out1,\"%.8lf  \", UM[i*2*NoM+j]);\n"
2004.07.14 03:28:00
<- "701\t\t\t}\n"
2004.07.14 03:28:00
<- "702\t\t\tfprintf(out1, \"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "703\t\t}\n"
2004.07.14 03:28:00
<- "704\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "705\t\t\n"
2004.07.14 03:28:00
<- "706\t\tout1 = fopen(\"vm.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "707\t\tfor (i=0; i<2*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "708\t\t\tfor (j=0; j<2*NoM; ++j){\n"
2004.07.14 03:28:00
<- "709\t\t\t\tfprintf(out1,\"%.8lf  \", VM[i*2*NoM+j]);\n"
2004.07.14 03:28:00
<- "710\t\t\t}\n"
2004.07.14 03:28:00
<- "711\t\t\tfprintf(out1, \"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "712\t\t}\n"
2004.07.14 03:28:00
<- "713\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "714\t\n"
2004.07.14 03:28:00
<- "715\t}\n"
2004.07.14 03:28:00
<- "716\t\n"
2004.07.14 03:28:00
<- "717\tvoid buildMatricesUV(void){\n"
2004.07.14 03:28:00
<- "718\t\tint i, j, k;\t\n"
2004.07.14 03:28:00
<- "719\t\tdouble normur, normui, normuvr, normuvi;\n"
2004.07.14 03:28:00
<- "720\t\tdouble *ur1, *ui1, *vr1, *vi1; \n"
2004.07.14 03:28:00
<- "721\t\tdouble *ur2, *ui2, *vr2, *vi2; \n"
2004.07.14 03:28:00
<- "722\t\t\n"
2004.07.14 03:28:00
<- "723\t\tUM = (double *)malloc((size_t)4*UVSize*UVSize*sizeof(double));\n"
2004.07.14 03:28:00
<- "724\t\tVM = (double *)malloc((size_t)4*UVSize*UVSize*sizeof(double));\n"
2004.07.14 03:28:00
<- "725\t\n"
2004.07.14 03:28:00
<- "726\t\tfor (i=0; i<UVSize; ++i){\n"
2004.07.14 03:28:00
<- "727\t\t\tur2 = uv + i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "728\t\t\tvr2 = ur2 + NoM;\n"
2004.07.14 03:28:00
<- "729\t\t\tui2 = vr2 + NoM;\n"
2004.07.14 03:28:00
<- "730\t\t\tvi2 = ui2 + NoM;\n"
2004.07.14 03:28:00
<- "731\t\t\tfor (j=0; j<UVSize; ++j){\n"
2004.07.14 03:28:00
<- "732\t\t\t\tur1 = uv + j*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "733\t\t\t\tvr1 = ur1 + NoM;\n"
2004.07.14 03:28:00
<- "734\t\t\t\tui1 = vr1 + NoM;\n"
2004.07.14 03:28:00
<- "735\t\t\t\tvi1 = ui1 + NoM;\n"
2004.07.14 03:28:00
<- "736\t\t\t\tnormur = 0.0;\n"
2004.07.14 03:28:00
<- "737\t\t\t\tnormui = 0.0;\n"
2004.07.14 03:28:00
<- "738\t\t\t\tnormuvr= 0.0;\n"
2004.07.14 03:28:00
<- "739\t\t\t\tnormuvi= 0.0;\n"
2004.07.14 03:28:00
<- "740\t\t\t\tfor (k=0; k<NoM; ++k){\n"
2004.07.14 03:28:00
<- "741\t\t\t\t\tnormur += ur2[k]*ur1[k]+ui2[k]*ui1[k];\n"
2004.07.14 03:28:00
<- "742\t\t\t\t\tnormui += ur2[k]*ui1[k]-ui2[k]*ur1[k];\n"
2004.07.14 03:28:00
<- "743\t\n"
2004.07.14 03:28:00
<- "744\t\t\t\t\tnormuvr += vr2[k]*ur1[k]+vi2[k]*ui1[k];\n"
2004.07.14 03:28:00
<- "745\t\t\t\t\tnormuvi -= vr2[k]*ui1[k]+vi2[k]*ur1[k];\n"
2004.07.14 03:28:00
<- "746\t\t\t\t}\n"
2004.07.14 03:28:00
<- "747\t\t\t\t\n"
2004.07.14 03:28:00
<- "748\t\t\t\tUM[i*2*UVSize+j] = normur;\n"
2004.07.14 03:28:00
<- "749\t\t\t\tUM[(i+UVSize)*2*UVSize+j+UVSize] = normur;\n"
2004.07.14 03:28:00
<- "750\t\t\t\tUM[i*2*UVSize+j+UVSize] = -normui;\n"
2004.07.14 03:28:00
<- "751\t\t\t\tUM[(i+UVSize)*2*UVSize+j] = normui;\n"
2004.07.14 03:28:00
<- "752\t\t\t\t\n"
2004.07.14 03:28:00
<- "753\t\t\t\tVM[i*2*UVSize+j] = normuvr;\n"
2004.07.14 03:28:00
<- "754\t\t\t\tVM[(i+UVSize)*2*UVSize+j+UVSize] = normuvr;\n"
2004.07.14 03:28:00
<- "755\t\t\t\tVM[i*2*UVSize+j+UVSize] = -normuvi;\n"
2004.07.14 03:28:00
<- "756\t\t\t\tVM[(i+UVSize)*2*UVSize+j] = normuvi;\n"
2004.07.14 03:28:00
<- "757\t\n"
2004.07.14 03:28:00
<- "758\t\t\t}\n"
2004.07.14 03:28:00
<- "759\t\t}\n"
2004.07.14 03:28:00
<- "760\t\t\n"
2004.07.14 03:28:00
<- "761\t\tout1 = fopen(\"um.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "762\t\tfor (i=0; i<2*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "763\t\t\tfor (j=0; j<2*UVSize; ++j){\n"
2004.07.14 03:28:00
<- "764\t\t\t\tfprintf(out1,\"%.8lf  \", UM[i*2*UVSize+j]);\n"
2004.07.14 03:28:00
<- "765\t\t\t}\n"
2004.07.14 03:28:00
<- "766\t\t\tfprintf(out1, \"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "767\t\t}\n"
2004.07.14 03:28:00
<- "768\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "769\t\t\n"
2004.07.14 03:28:00
<- "770\t\tout1 = fopen(\"vm.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "771\t\tfor (i=0; i<2*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "772\t\t\tfor (j=0; j<2*UVSize; ++j){\n"
2004.07.14 03:28:00
<- "773\t\t\t\tfprintf(out1,\"%.8lf  \", VM[i*2*UVSize+j]);\n"
2004.07.14 03:28:00
<- "774\t\t\t}\n"
2004.07.14 03:28:00
<- "775\t\t\tfprintf(out1, \"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "776\t\t}\n"
2004.07.14 03:28:00
<- "777\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "778\t\n"
2004.07.14 03:28:00
<- "779\t}\n"
2004.07.14 03:28:00
<- "780\t\n"
2004.07.14 03:28:00
<- "781\tvoid invertUOrthoBasis(void){\n"
2004.07.14 03:28:00
<- "782\t\tint i, j, k;\t\n"
2004.07.14 03:28:00
<- "783\t\tdouble *test;\n"
2004.07.14 03:28:00
<- "784\t\tdouble tmp;\n"
2004.07.14 03:28:00
<- "785\t\n"
2004.07.14 03:28:00
<- "786\t\tint m = 2*UVSize, n = 2*NoM, lda = m;\n"
2004.07.14 03:28:00
<- "787\t\tint *ipiv;\n"
2004.07.14 03:28:00
<- "788\t\tint info;\n"
2004.07.14 03:28:00
<- "789\t\tdouble *work;\n"
2004.07.14 03:28:00
<- "790\t\tint lwork = m;\n"
2004.07.14 03:28:00
<- "791\t\n"
2004.07.14 03:28:00
<- "792\t\tipiv = (int *)malloc((size_t)m*sizeof(int));\n"
2004.07.14 03:28:00
<- "793\t\twork = (double *)malloc((size_t)lwork*sizeof(double));\n"
2004.07.14 03:28:00
<- "794\t\t\n"
2004.07.14 03:28:00
<- "795\t\n"
2004.07.14 03:28:00
<- "796\t\ttest = (double *)malloc((size_t)4*UVSize*NoM*sizeof(double));\n"
2004.07.14 03:28:00
<- "797\t\tfor (i=0; i<4*UVSize*NoM; ++i){\n"
2004.07.14 03:28:00
<- "798\t\t\ttest[i] = UM[i];\n"
2004.07.14 03:28:00
<- "799\t\t}\n"
2004.07.14 03:28:00
<- "800\t\n"
2004.07.14 03:28:00
<- "801\t\n"
2004.07.14 03:28:00
<- "802\t\tdgetrf_(&m, &n, UM, &lda, ipiv, &info);\n"
2004.07.14 03:28:00
<- "803\t\tprintf(\"dgetrf(): info=%d\\t\", info);\n"
2004.07.14 03:28:00
<- "804\t\n"
2004.07.14 03:28:00
<- "805\t\tdgetri_(&n, UM, &m, ipiv, work, &lwork, &info);\n"
2004.07.14 03:28:00
<- "806\t\tprintf(\"dgetri(): info=%d\\n"
   "\", info);\n"
2004.07.14 03:28:00
<- "807\t\t\n"
2004.07.14 03:28:00
<- "808\t\tout1 = fopen(\"invTest.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "809\t\tfor (i=0; i<2*NoM; ++i) {\n"
2004.07.14 03:28:00
<- "810\t\t\tfor (j=0; j<2*NoM; ++j) {\n"
2004.07.14 03:28:00
<- "811\t\t\t\ttmp = 0.0;\n"
2004.07.14 03:28:00
<- "812\t\t\t\tfor (k=0; k<2*UVSize; ++k) {\n"
2004.07.14 03:28:00
<- "813\t\t\t\t\ttmp += UM[i*2*UVSize+k]*test[k*2*NoM+j];\n"
2004.07.14 03:28:00
<- "814\t\t\t\t}\n"
2004.07.14 03:28:00
<- "815\t\t\t\tfprintf(out1, \"%.8lf\\t\", tmp);\n"
2004.07.14 03:28:00
<- "816\t\t\t}\n"
2004.07.14 03:28:00
<- "817\t\t\tfprintf(out1,\"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "818\t\t}\n"
2004.07.14 03:28:00
<- "819\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "820\t\n"
2004.07.14 03:28:00
<- "821\t\n"
2004.07.14 03:28:00
<- "822\t\tfree(ipiv);\n"
2004.07.14 03:28:00
<- "823\t\tfree(work);\n"
2004.07.14 03:28:00
<- "824\t}\n"
2004.07.14 03:28:00
<- "825\t\n"
2004.07.14 03:28:00
<- "826\t\n"
2004.07.14 03:28:00
<- "827\tvoid invertU(void){\n"
2004.07.14 03:28:00
<- "828\t\tint i, j;\t\n"
2004.07.14 03:28:00
<- "829\t\tdouble *b;\n"
2004.07.14 03:28:00
<- "830\t\tdouble *test;\n"
2004.07.14 03:28:00
<- "831\t\tdouble tmp;\n"
2004.07.14 03:28:00
<- "832\t\tint k;\n"
2004.07.14 03:28:00
<- "833\t/*\n"
2004.07.14 03:28:00
<- "834\t\ttest = (double *)malloc((size_t)4*UVSize*UVSize*sizeof(double));\n"
2004.07.14 03:28:00
<- "835\t\tfor (i=0; i<4*UVSize*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "836\t\t\ttest[i] = UM[i];\n"
2004.07.14 03:28:00
<- "837\t\t}\n"
2004.07.14 03:28:00
<- "838\t*/\tb = (double *)malloc((size_t)2*UVSize*sizeof(double));\n"
2004.07.14 03:28:00
<- "839\t\t\n"
2004.07.14 03:28:00
<- "840\t\tfor (i=0; i<2*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "841\t\t\tb[i] = 0.0;\n"
2004.07.14 03:28:00
<- "842\t\t}\n"
2004.07.14 03:28:00
<- "843\t/*Numerical recipes function gaussj has been changed from the original\n"
2004.07.14 03:28:00
<- "844\t * Now it is enough to pass beginings of arrays UM and b.*/\n"
2004.07.14 03:28:00
<- "845\t\tgaussj(UM, 2*UVSize, b, 1);\n"
2004.07.14 03:28:00
<- "846\t\n"
2004.07.14 03:28:00
<- "847\t/*\t\n"
2004.07.14 03:28:00
<- "848\t\tfor (i=0; i<2*UVSize; ++i) {\n"
2004.07.14 03:28:00
<- "849\t\t\tfor (j=0; j<2*UVSize; ++j) {\n"
2004.07.14 03:28:00
<- "850\t\t\t\ttmp = 0.0;\n"
2004.07.14 03:28:00
<- "851\t\t\t\tfor (k=0; k<2*UVSize; ++k) {\n"
2004.07.14 03:28:00
<- "852\t\t\t\t\ttmp += UM[k*2*UVSize+j]*test[i*2*UVSize+k];\n"
2004.07.14 03:28:00
<- "853\t\t\t\t}\n"
2004.07.14 03:28:00
<- "854\t\t\t\tprintf(\"%.8lf\\t\", tmp);\n"
2004.07.14 03:28:00
<- "855\t\t\t}\n"
2004.07.14 03:28:00
<- "856\t\t\tprintf(\"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "857\t\t}\n"
2004.07.14 03:28:00
<- "858\t*/\n"
2004.07.14 03:28:00
<- "859\t\t\n"
2004.07.14 03:28:00
<- "860\t\tfree(b);\n"
2004.07.14 03:28:00
<- "861\t}\n"
2004.07.14 03:28:00
<- "862\t\n"
2004.07.14 03:28:00
<- "863\tvoid createZM(void){\n"
2004.07.14 03:28:00
<- "864\t\tint i, j, k;\n"
2004.07.14 03:28:00
<- "865\t\tdouble tmpr, tmpi;\n"
2004.07.14 03:28:00
<- "866\t\t\n"
2004.07.14 03:28:00
<- "867\t\tZM = (double *)malloc((size_t)4*UVSize*UVSize*sizeof(double));\n"
2004.07.14 03:28:00
<- "868\t\n"
2004.07.14 03:28:00
<- "869\t\tfor (i=0; i<2*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "870\t\t\tfor (j=0; j<2*UVSize; ++j){\n"
2004.07.14 03:28:00
<- "871\t\t\t\ttmpr = 0.0;\n"
2004.07.14 03:28:00
<- "872\t\t\t\tfor (k=0; k<2*UVSize; ++k) {\n"
2004.07.14 03:28:00
<- "873\t\t\t\t\ttmpr += UM[i*2*UVSize + k]*VM[k*2*UVSize + j];\n"
2004.07.14 03:28:00
<- "874\t\t\t\t}\n"
2004.07.14 03:28:00
<- "875\t\t\t\tZM[i*2*UVSize+j] = tmpr;\t\n"
2004.07.14 03:28:00
<- "876\t\t\t}\n"
2004.07.14 03:28:00
<- "877\t\t}\n"
2004.07.14 03:28:00
<- "878\t/*\t\n"
2004.07.14 03:28:00
<- "879\t\tfor (i=0; i<UVSize; ++i){\n"
2004.07.14 03:28:00
<- "880\t\t\tfor (j=0; j<UVSize; ++j){\n"
2004.07.14 03:28:00
<- "881\t\t\t\ttmpr = 0.0;\n"
2004.07.14 03:28:00
<- "882\t\t\t\ttmpi = 0.0;\n"
2004.07.14 03:28:00
<- "883\t\t\t\tfor (k=0; k<UVSize; ++k) {\n"
2004.07.14 03:28:00
<- "884\t\t\t\t\ttmpr+= UM[i*2*UVSize + k]*VM[k*2*UVSize + j]-UM[i*2*UVSize +k+UVSize]*VM[(k+UVSize)*2*UVSize + j];\n"
2004.07.14 03:28:00
<- "885\t\t\t\t\ttmpi-= UM[i*2*UVSize + k+UVSize]*VM[k*2*UVSize + j]+UM[i*2*UVSize +k]*VM[(k+UVSize)*2*UVSize + j];\n"
2004.07.14 03:28:00
<- "886\t\t\t\t}\n"
2004.07.14 03:28:00
<- "887\t\t\t\tZM[i*2*UVSize+j] = tmpr;\t\n"
2004.07.14 03:28:00
<- "888\t\t\t\tZM[(i+UVSize)*2*UVSize+j+UVSize] = tmpr;\t\n"
2004.07.14 03:28:00
<- "889\t\t\t\tZM[i*2*UVSize+j+UVSize] = -tmpi;\t\n"
2004.07.14 03:28:00
<- "890\t\t\t\tZM[(i+UVSize)*2*UVSize+j] = tmpi;\t\n"
2004.07.14 03:28:00
<- "891\t\t\t}\n"
2004.07.14 03:28:00
<- "892\t\t}\n"
2004.07.14 03:28:00
<- "893\t*/\n"
2004.07.14 03:28:00
<- "894\t\tout1 = fopen(\"zm.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "895\t\tfor (i=0; i<2*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "896\t\t\tfor (j=0; j<2*UVSize; ++j){\n"
2004.07.14 03:28:00
<- "897\t\t\t\tfprintf(out1,\"%.8lf  \", ZM[i*2*UVSize+j]);\n"
2004.07.14 03:28:00
<- "898\t\t\t}\n"
2004.07.14 03:28:00
<- "899\t\t\tfprintf(out1, \"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "900\t\t}\n"
2004.07.14 03:28:00
<- "901\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "902\t\t\n"
2004.07.14 03:28:00
<- "903\t\tout1 = fopen(\"umo.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "904\t\tfor (i=0; i<2*UVSize; ++i){\n"
2004.07.14 03:28:00
<- "905\t\t\tfor (j=0; j<2*UVSize; ++j){\n"
2004.07.14 03:28:00
<- "906\t\t\t\tfprintf(out1,\"%.8lf  \", UM[i*2*UVSize+j]);\n"
2004.07.14 03:28:00
<- "907\t\t\t}\n"
2004.07.14 03:28:00
<- "908\t\t\tfprintf(out1, \"\\n"
   "\\n"
   "\");\n"
2004.07.14 03:28:00
<- "909\t\t}\n"
2004.07.14 03:28:00
<- "910\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "911\t\n"
2004.07.14 03:28:00
<- "912\t}\n"
2004.07.14 03:28:00
<- "913\t\n"
2004.07.14 03:28:00
<- "914\t\n"
2004.07.14 03:28:00
<- "915\tint diagonalize2(void){\n"
2004.07.14 03:28:00
<- "916\t\tchar balanc=\'N\', jobvl=\'N\';\n"
2004.07.14 03:28:00
<- "917\t\tchar jobvr=\'V\', sense=\'N\';\n"
2004.07.14 03:28:00
<- "918\t\tint il, iu, info, n=2*UVSize; //il>=1 iu=il+SizeOfBase-1\n"
2004.07.14 03:28:00
<- "919\t\tint lda=n, lwork=8*n;\n"
2004.07.14 03:28:00
<- "920\t\tdouble abnrm;\n"
2004.07.14 03:28:00
<- "921\t\tint *iwork;\n"
2004.07.14 03:28:00
<- "922\t\tdouble *scale,*rconde, *rcondv, *work;\n"
2004.07.14 03:28:00
<- "923\t\n"
2004.07.14 03:28:00
<- "924\t\tLambda = (double *)malloc((size_t)2*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "925\t\tBogoModes = (double *)malloc((size_t)n*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "926\t\n"
2004.07.14 03:28:00
<- "927\t\tiwork=(int *)malloc((size_t)2*n*sizeof(int));\n"
2004.07.14 03:28:00
<- "928\t\tscale=(double *)malloc((size_t)n*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "929\t\trconde=(double *)malloc((size_t)n*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "930\t\trcondv=(double *)malloc((size_t)n*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "931\t\twork=(double *)malloc((size_t)lwork*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "932\t\n"
2004.07.14 03:28:00
<- "933\t\tdgeevx_( &balanc, &jobvl, &jobvr, &sense, &n, ZM, &lda, Lambda, Lambda+n, BogoModes, &n, BogoModes, &n, &il, &iu, scale, &abnrm, rconde, rcondv, work, &lwork, iwork, &info );\n"
2004.07.14 03:28:00
<- "934\t\n"
2004.07.14 03:28:00
<- "935\t\tif(info){\n"
2004.07.14 03:28:00
<- "936\t\t\tfprintf(stderr,\"info = %d\\n"
   "\", info);\n"
2004.07.14 03:28:00
<- "937\t\t\tfprintf(stderr,\"diagonalize2(): Problem with diagonalization! \\n"
   "\");\n"
2004.07.14 03:28:00
<- "938\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "939\t\t}\n"
2004.07.14 03:28:00
<- "940\t\n"
2004.07.14 03:28:00
<- "941\t\tfree(scale);\n"
2004.07.14 03:28:00
<- "942\t\tfree(rconde);\n"
2004.07.14 03:28:00
<- "943\t\tfree(rcondv);\n"
2004.07.14 03:28:00
<- "944\t\tfree(iwork);\n"
2004.07.14 03:28:00
<- "945\t\tfree(work);\n"
2004.07.14 03:28:00
<- "946\t\treturn info;\n"
2004.07.14 03:28:00
<- "947\t}\n"
2004.07.14 03:28:00
<- "948\t\n"
2004.07.14 03:28:00
<- "949\tint findBogoModes(void){\n"
2004.07.14 03:28:00
<- "950\t\tchar jobz=\'V\', range=\'I\', uplo=\'U\';\n"
2004.07.14 03:28:00
<- "951\t\tint il=1, iu=2*M, info, m, n=2*UVSize; //il>=1 iu=il+SizeOfBase-1\n"
2004.07.14 03:28:00
<- "952\t\tint lda=n, ldz=n, lwork=8*n;\n"
2004.07.14 03:28:00
<- "953\t\tdouble abstol=acc, vl=50.0, vu=60.0;\n"
2004.07.14 03:28:00
<- "954\t\tint *ifail, *iwork;\n"
2004.07.14 03:28:00
<- "955\t\tdouble *work;\n"
2004.07.14 03:28:00
<- "956\t\n"
2004.07.14 03:28:00
<- "957\t\tifail=(int *)malloc((size_t)n*sizeof(int));\n"
2004.07.14 03:28:00
<- "958\t\tiwork=(int *)malloc((size_t)5*n*sizeof(int));\n"
2004.07.14 03:28:00
<- "959\t\twork=(double *)malloc((size_t)lwork*sizeof(double));\n"
2004.07.14 03:28:00
<- "960\t\tLambda = (double *)malloc((size_t)n*sizeof(double));\n"
2004.07.14 03:28:00
<- "961\t\tBogoModes = (double *)malloc((size_t)2*M*n*sizeof(double));\n"
2004.07.14 03:28:00
<- "962\t\n"
2004.07.14 03:28:00
<- "963\t\tdsyevx_( &jobz, &range, &uplo, &n, ZM, &lda, &vl, &vu,  &il,  &iu, &abstol,  &m,  Lambda,  BogoModes,  &ldz,  work, &lwork, iwork, ifail, &info );\n"
2004.07.14 03:28:00
<- "964\t\n"
2004.07.14 03:28:00
<- "965\t\tfree(ifail);\n"
2004.07.14 03:28:00
<- "966\t\tfree(iwork);\n"
2004.07.14 03:28:00
<- "967\t\tfree(work);\n"
2004.07.14 03:28:00
<- "968\t\treturn info;\n"
2004.07.14 03:28:00
<- "969\t}\n"
2004.07.14 03:28:00
<- "970\t\n"
2004.07.14 03:28:00
<- "971\tvoid fillBogoModesX(void){\n"
2004.07.14 03:28:00
<- "972\t\tint i,j,k,l, q;\n"
2004.07.14 03:28:00
<- "973\t\tdouble *vec;\n"
2004.07.14 03:28:00
<- "974\t\tchar buff[20];\n"
2004.07.14 03:28:00
<- "975\t \tdouble x;\n"
2004.07.14 03:28:00
<- "976\t\tdouble tmp;\n"
2004.07.14 03:28:00
<- "977\t\tdouble *bMr, *bMi;\n"
2004.07.14 03:28:00
<- "978\t\tdouble *ur, *ui;\n"
2004.07.14 03:28:00
<- "979\t\tint m;\n"
2004.07.14 03:28:00
<- "980\t\n"
2004.07.14 03:28:00
<- "981\t\tif(UVSize>=M) m = M;\n"
2004.07.14 03:28:00
<- "982\t\telse m = UVSize;\n"
2004.07.14 03:28:00
<- "983\t\t\n"
2004.07.14 03:28:00
<- "984\t\tBogoModesX = (double *)malloc((size_t)4*N*M*sizeof(double));\n"
2004.07.14 03:28:00
<- "985\t\tvec = (double *)malloc((size_t)2*N*sizeof(double));\n"
2004.07.14 03:28:00
<- "986\t\n"
2004.07.14 03:28:00
<- "987\t\tfor (i=0; i<2*m; i+=1) {\n"
2004.07.14 03:28:00
<- "988\t\t\tfor (k=0; k<2*N; ++k) {\n"
2004.07.14 03:28:00
<- "989\t\t\t\tvec[k] = 0.0;\n"
2004.07.14 03:28:00
<- "990\t\t\t}\n"
2004.07.14 03:28:00
<- "991\t\t\tbMr = BogoModes+i*2*UVSize;\n"
2004.07.14 03:28:00
<- "992\t\t\tbMi = bMr+UVSize;\n"
2004.07.14 03:28:00
<- "993\t\t\tfor (j=0; j<UVSize; ++j) {\n"
2004.07.14 03:28:00
<- "994\t\t\t\tur = uv+j*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "995\t\t\t\tui = ur+2*NoM;\n"
2004.07.14 03:28:00
<- "996\t\t\t\tfor (k=0; k<NoM; ++k) {\n"
2004.07.14 03:28:00
<- "997\t\t\t\t\tfor(l=0; l<N; ++l){\n"
2004.07.14 03:28:00
<- "998\t\t\t\t\t\tvec[l] += (bMr[j]*ur[k]-bMi[j]*ui[k])*BoxModes[k*N+l];\n"
2004.07.14 03:28:00
<- "999\t\t\t\t\t\tvec[l+N] += (bMr[j]*ui[k]+bMi[j]*ur[k])*BoxModes[k*N+l];\n"
2004.07.14 03:28:00
<- "1000\t\t\t\t\t}\n"
2004.07.14 03:28:00
<- "1001\t\t\t\t}\n"
2004.07.14 03:28:00
<- "1002\t\t\t}\n"
2004.07.14 03:28:00
<- "1003\t\t\ttmp = 0.0;\n"
2004.07.14 03:28:00
<- "1004\t\t\tfor (k=0; k<N; ++k) {\n"
2004.07.14 03:28:00
<- "1005\t\t\t\ttmp+= vec[k]*vec[k]+vec[k+N]*vec[k+N];\n"
2004.07.14 03:28:00
<- "1006\t\t\t}\n"
2004.07.14 03:28:00
<- "1007\t//\t\ttmp = 1.0/sqrt(tmp*dx);  //normalization also performed in BgControl.c\n"
2004.07.14 03:28:00
<- "1008\t\t\ttmp = 1.0/sqrt(tmp);\n"
2004.07.14 03:28:00
<- "1009\t\t\tfor (k=0; k<2*N; ++k) {\n"
2004.07.14 03:28:00
<- "1010\t\t\t\tBogoModesX[i*N + k] = vec[k]*tmp;\n"
2004.07.14 03:28:00
<- "1011\t\t\t}\n"
2004.07.14 03:28:00
<- "1012\t\t}\n"
2004.07.14 03:28:00
<- "1013\t\n"
2004.07.14 03:28:00
<- "1014\t\tfree(vec);\n"
2004.07.14 03:28:00
<- "1015\t}\n"
2004.07.14 03:28:00
<- "1016\t\n"
2004.07.14 03:28:00
<- "1017\tvoid outBogoModes(void){\n"
2004.07.14 03:28:00
<- "1018\t\tint i,j,k;\n"
2004.07.14 03:28:00
<- "1019\t\tchar buff[20];\n"
2004.07.14 03:28:00
<- "1020\t \tdouble x;\n"
2004.07.14 03:28:00
<- "1021\t\tdouble tmp;\n"
2004.07.14 03:28:00
<- "1022\t\tint m;\n"
2004.07.14 03:28:00
<- "1023\t\n"
2004.07.14 03:28:00
<- "1024\t\tif(UVSize>=M) m = M;\n"
2004.07.14 03:28:00
<- "1025\t\telse m = UVSize;\n"
2004.07.14 03:28:00
<- "1026\t\t\n"
2004.07.14 03:28:00
<- "1027\t\tfor (i=0; i<m; ++i) {\n"
2004.07.14 03:28:00
<- "1028\t\t\tsprintf(buff, \"mod%d.dat\", i+1);\n"
2004.07.14 03:28:00
<- "1029\t\t\tout1 = fopen(buff, \"wt\");\n"
2004.07.14 03:28:00
<- "1030\t\t\tx = xL;\n"
2004.07.14 03:28:00
<- "1031\t\t\tfor (k=0; k<N; ++k) {\n"
2004.07.14 03:28:00
<- "1032\t\t\t\ttmp = BogoModesX[i*N+k]*BogoModesX[i*N+k]+BogoModesX[i*N+k+N]*BogoModesX[i*N+k+N];\n"
2004.07.14 03:28:00
<- "1033\t\t\t\tfprintf(out1, \"%lf\\t%.12lf\\t%.12lf\\t%.12lf\\n"
   "\", x, BogoModesX[i*N+k], BogoModesX[i*N+k+N], tmp);\n"
2004.07.14 03:28:00
<- "1034\t\t\t\tx += dx;\n"
2004.07.14 03:28:00
<- "1035\t\t\t}\n"
2004.07.14 03:28:00
<- "1036\t\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1037\t\t}\n"
2004.07.14 03:28:00
<- "1038\t\n"
2004.07.14 03:28:00
<- "1039\t}\n"
2004.07.14 03:28:00
<- "1040\t\n"
2004.07.14 03:28:00
<- "1041\tvoid testH(void){\n"
2004.07.14 03:28:00
<- "1042\t\tint i, j, k;\n"
2004.07.14 03:28:00
<- "1043\t\tdouble *vec;\n"
2004.07.14 03:28:00
<- "1044\t\n"
2004.07.14 03:28:00
<- "1045\t\tvec = (double *)malloc((size_t)2*SizeOfBase*sizeof(double));\n"
2004.07.14 03:28:00
<- "1046\t\tfor (i=0; i<2*SizeOfBase; ++i) {\n"
2004.07.14 03:28:00
<- "1047\t\t\tvec[i] = 0.0;\n"
2004.07.14 03:28:00
<- "1048\t\t}\n"
2004.07.14 03:28:00
<- "1049\t\t\n"
2004.07.14 03:28:00
<- "1050\t\tfillHamiltonian(H);\n"
2004.07.14 03:28:00
<- "1051\t\tfor (i=0; i<SizeOfBase; ++i) {\n"
2004.07.14 03:28:00
<- "1052\t\t\tfor(j=0; j<NoM; ++j){\n"
2004.07.14 03:28:00
<- "1053\t\t\t\tvec[i] += H[i*2*SizeOfBase +j]*boxPhi[j];\n"
2004.07.14 03:28:00
<- "1054\t\t\t\tvec[i+SizeOfBase] += H[i*2*SizeOfBase +j]*boxPhi[j+NoM];\n"
2004.07.14 03:28:00
<- "1055\t\t\t}\n"
2004.07.14 03:28:00
<- "1056\t\t\tfor(j=NoM; j<SizeOfBase; ++j){\n"
2004.07.14 03:28:00
<- "1057\t\t\t\tvec[i] += H[i*2*SizeOfBase +j]*boxPhi[j-NoM];\n"
2004.07.14 03:28:00
<- "1058\t\t\t\tvec[i] += H[i*2*SizeOfBase +j]*boxPhi[j-NoM];\n"
2004.07.14 03:28:00
<- "1059\t\t\t}\n"
2004.07.14 03:28:00
<- "1060\t\t}\n"
2004.07.14 03:28:00
<- "1061\t\n"
2004.07.14 03:28:00
<- "1062\t\tout1 = fopen(\"Htest.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "1063\t\tfor (i=0; i<SizeOfBase; ++i) {\n"
2004.07.14 03:28:00
<- "1064\t\t\tfprintf(out1, \"%d\\t%lf\\n"
   "\", i, vec[i]);\n"
2004.07.14 03:28:00
<- "1065\t\t}\n"
2004.07.14 03:28:00
<- "1066\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1067\t\n"
2004.07.14 03:28:00
<- "1068\t\tfree(vec);\n"
2004.07.14 03:28:00
<- "1069\t}\n"
2004.07.14 03:28:00
<- "1070\t\n"
2004.07.14 03:28:00
<- "1071\tvoid outputs(void){\n"
2004.07.14 03:28:00
<- "1072\t\tint i, j, k;\n"
2004.07.14 03:28:00
<- "1073\t\tdouble tmp;\n"
2004.07.14 03:28:00
<- "1074\t\tdouble x;\n"
2004.07.14 03:28:00
<- "1075\t\tint SoB = 2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1076\t\tchar buff[20];\n"
2004.07.14 03:28:00
<- "1077\t\tint m;\n"
2004.07.14 03:28:00
<- "1078\t\n"
2004.07.14 03:28:00
<- "1079\t\tif(UVSize>=M) m = M;\n"
2004.07.14 03:28:00
<- "1080\t\telse m = UVSize;\n"
2004.07.14 03:28:00
<- "1081\t\n"
2004.07.14 03:28:00
<- "1082\t\toutBogoModes();\n"
2004.07.14 03:28:00
<- "1083\t\n"
2004.07.14 03:28:00
<- "1084\t\tout1 = fopen(\"uv.bin\", \"wb\");\n"
2004.07.14 03:28:00
<- "1085\t\tif (fwrite(uv, 8, SoB*UVSize, out1) != SoB*UVSize) {\n"
2004.07.14 03:28:00
<- "1086\t\t\tfprintf(stderr, \"bogo: Error writing uv! \\n"
   "\");\n"
2004.07.14 03:28:00
<- "1087\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "1088\t\t}\n"
2004.07.14 03:28:00
<- "1089\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1090\t\t\n"
2004.07.14 03:28:00
<- "1091\t\tout1 = fopen(\"lambda.bin\", \"wb\");\n"
2004.07.14 03:28:00
<- "1092\t\tfwrite(Lambda, 8, 2*m, out1);\n"
2004.07.14 03:28:00
<- "1093\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1094\t\t\n"
2004.07.14 03:28:00
<- "1095\t\tout1 = fopen(\"lambda.dat\", \"wt\");\n"
2004.07.14 03:28:00
<- "1096\t\tfor (i=0; i<2*m; ++i){\n"
2004.07.14 03:28:00
<- "1097\t\t\tfprintf(out1, \"%d\\t%.10lf\\t%.10lf\\n"
   "\", i, Lambda[i], Lambda[i+2*UVSize]);\n"
2004.07.14 03:28:00
<- "1098\t\t}\n"
2004.07.14 03:28:00
<- "1099\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1100\t\n"
2004.07.14 03:28:00
<- "1101\t\t/* This is to undo normalization*/\n"
2004.07.14 03:28:00
<- "1102\t\ttmp = 0.0;\n"
2004.07.14 03:28:00
<- "1103\t\tfor (i=0; i<N; ++i){\n"
2004.07.14 03:28:00
<- "1104\t\t\ttmp += phi[i]*phi[i];\n"
2004.07.14 03:28:00
<- "1105\t\t}\n"
2004.07.14 03:28:00
<- "1106\t\ttmp = 1.0/sqrt(tmp);\n"
2004.07.14 03:28:00
<- "1107\t\tfor (i=0; i<N; ++i){\n"
2004.07.14 03:28:00
<- "1108\t\t\t phi[i]*=tmp;\n"
2004.07.14 03:28:00
<- "1109\t\t}\n"
2004.07.14 03:28:00
<- "1110\t\tout1 = fopen(\"BDSModesX.bin\", \"wb\");\n"
2004.07.14 03:28:00
<- "1111\t\tif (fwrite(phi, 8, N, out1) != N) {\n"
2004.07.14 03:28:00
<- "1112\t\t\tfprintf(stderr, \"bogo: Error writing phi0! \\n"
   "\");\n"
2004.07.14 03:28:00
<- "1113\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "1114\t\t}\n"
2004.07.14 03:28:00
<- "1115\t\tif (fwrite(BogoModesX, 8, N*(M-1), out1) != N*(M-1)) {\n"
2004.07.14 03:28:00
<- "1116\t\t\tfprintf(stderr, \"bogo: Error writing BogoModes! \\n"
   "\");\n"
2004.07.14 03:28:00
<- "1117\t\t\texit(1);\n"
2004.07.14 03:28:00
<- "1118\t\t}\n"
2004.07.14 03:28:00
<- "1119\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1120\t/*\n"
2004.07.14 03:28:00
<- "1121\t\tout1 = fopen(\"BDSModes.bin\", \"wb\");\n"
2004.07.14 03:28:00
<- "1122\t\ttmp = 1.0;\n"
2004.07.14 03:28:00
<- "1123\t\tfwrite(&tmp, 8, 1, out1);\n"
2004.07.14 03:28:00
<- "1124\t\ttmp = 0.0;\n"
2004.07.14 03:28:00
<- "1125\t\tfor (i=0; i<M-1; ++i){\n"
2004.07.14 03:28:00
<- "1126\t\t\tfwrite(&tmp, 8, 1, out1);\n"
2004.07.14 03:28:00
<- "1127\t\t}\n"
2004.07.14 03:28:00
<- "1128\t\tfor (i=0; i<M-1; ++i){\n"
2004.07.14 03:28:00
<- "1129\t\t\tfwrite(&tmp, 8, 1, out1);\n"
2004.07.14 03:28:00
<- "1130\t\t\tfwrite(BogoModes+i*M, 8, M-1, out1);\n"
2004.07.14 03:28:00
<- "1131\t\t} \n"
2004.07.14 03:28:00
<- "1132\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1133\t*/\n"
2004.07.14 03:28:00
<- "1134\t}\n"
2004.07.14 03:28:00
<- "1135\t\n"
2004.07.14 03:28:00
<- "1136\tvoid fillBoxModes(void){\n"
2004.07.14 03:28:00
<- "1137\t\tint i,j;\n"
2004.07.14 03:28:00
<- "1138\t\tdouble x;\n"
2004.07.14 03:28:00
<- "1139\t\tdouble tmp1, tmp2, ar;\n"
2004.07.14 03:28:00
<- "1140\t\tdouble L=(xR-xL)*0.5, l;\n"
2004.07.14 03:28:00
<- "1141\t\n"
2004.07.14 03:28:00
<- "1142\t\n"
2004.07.14 03:28:00
<- "1143\t\tl=M_PI/L/2.0;\n"
2004.07.14 03:28:00
<- "1144\t\n"
2004.07.14 03:28:00
<- "1145\t\ttmp1 = sqrt(dx/L);\n"
2004.07.14 03:28:00
<- "1146\t\ttmp2 = sqrt(dx/L*0.5);\n"
2004.07.14 03:28:00
<- "1147\t\tfor(i=0; i<N; ++i){\n"
2004.07.14 03:28:00
<- "1148\t\t\tx=xL+i*dx;\n"
2004.07.14 03:28:00
<- "1149\t\t\tBoxModes[i] = tmp2;\n"
2004.07.14 03:28:00
<- "1150\t\t\tfor(j=1; j<NrOfBoxModes; j+=2){\n"
2004.07.14 03:28:00
<- "1151\t\t\t\tar = (double)(j+1)*l; \n"
2004.07.14 03:28:00
<- "1152\t\t\t\tBoxModes[j*N+i] = sin(ar*x)*tmp1;\n"
2004.07.14 03:28:00
<- "1153\t\t\t\tBoxModes[(j+1)*N+i] = cos(ar*x)*tmp1;\n"
2004.07.14 03:28:00
<- "1154\t\t\t}\n"
2004.07.14 03:28:00
<- "1155\t\t}\n"
2004.07.14 03:28:00
<- "1156\t}\n"
2004.07.14 03:28:00
<- "1157\t\n"
2004.07.14 03:28:00
<- "1158\tvoid outUVX(void){\n"
2004.07.14 03:28:00
<- "1159\t\tint i,j,k;\n"
2004.07.14 03:28:00
<- "1160\t\tdouble *vec;\n"
2004.07.14 03:28:00
<- "1161\t\tchar buff[20];\n"
2004.07.14 03:28:00
<- "1162\t \tdouble x;\n"
2004.07.14 03:28:00
<- "1163\t\tdouble tmp;\n"
2004.07.14 03:28:00
<- "1164\t\tdouble *bMr, *bMi;\n"
2004.07.14 03:28:00
<- "1165\t\t\n"
2004.07.14 03:28:00
<- "1166\t\tvec = (double *)malloc((size_t)2*N*sizeof(double));\n"
2004.07.14 03:28:00
<- "1167\t\n"
2004.07.14 03:28:00
<- "1168\t//\tfor (i=0; i<UVSize; ++i) {\n"
2004.07.14 03:28:00
<- "1169\t\tfor (i=0; i<2*NoM; ++i) {\n"
2004.07.14 03:28:00
<- "1170\t\t\tfor (k=0; k<2*N; ++k) {\n"
2004.07.14 03:28:00
<- "1171\t\t\t\tvec[k] = 0.0;\n"
2004.07.14 03:28:00
<- "1172\t\t\t}\n"
2004.07.14 03:28:00
<- "1173\t\t\tbMr = z+(indx[i]-1)*4*NoM;\n"
2004.07.14 03:28:00
<- "1174\t\t\tbMi = bMr+2*NoM;\n"
2004.07.14 03:28:00
<- "1175\t\t\tfor (j=0; j<NoM; ++j) {\n"
2004.07.14 03:28:00
<- "1176\t\t\t\tfor(k=0; k<N; ++k){\n"
2004.07.14 03:28:00
<- "1177\t\t\t\t\tvec[k] += bMr[j]*BoxModes[j*N+k];\n"
2004.07.14 03:28:00
<- "1178\t\t\t\t\tvec[k+N] += bMi[j]*BoxModes[j*N+k];\n"
2004.07.14 03:28:00
<- "1179\t\t\t\t}\n"
2004.07.14 03:28:00
<- "1180\t\t\t}\n"
2004.07.14 03:28:00
<- "1181\t\t\tfor(k=0; k<20; ++k) buff[k] = 0;\n"
2004.07.14 03:28:00
<- "1182\t\t\tsprintf(buff, \"u%d.dat\", i);\n"
2004.07.14 03:28:00
<- "1183\t\t\tout1 = fopen(buff, \"wt\");\n"
2004.07.14 03:28:00
<- "1184\t\t\tfor(k=0; k<N; ++k){\n"
2004.07.14 03:28:00
<- "1185\t\t\t\tx=xL+k*dx;\n"
2004.07.14 03:28:00
<- "1186\t\t\t\ttmp = vec[k]*vec[k]+vec[k+N]*vec[k+N];\n"
2004.07.14 03:28:00
<- "1187\t\t\t\tfprintf(out1, \"%lf\\t%.12lf\\t%.12lf\\t%.12lf\\n"
   "\", x, vec[k], vec[k+N], tmp);\n"
2004.07.14 03:28:00
<- "1188\t\t\t}\n"
2004.07.14 03:28:00
<- "1189\t\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1190\t\t\t\n"
2004.07.14 03:28:00
<- "1191\t\t}\n"
2004.07.14 03:28:00
<- "1192\t\t\n"
2004.07.14 03:28:00
<- "1193\t//\tfor (i=0; i<UVSize; ++i) {\n"
2004.07.14 03:28:00
<- "1194\t\tfor (i=0; i<2*NoM; ++i) {\n"
2004.07.14 03:28:00
<- "1195\t\t\tfor (k=0; k<2*N; ++k) {\n"
2004.07.14 03:28:00
<- "1196\t\t\t\tvec[k] = 0.0;\n"
2004.07.14 03:28:00
<- "1197\t\t\t}\n"
2004.07.14 03:28:00
<- "1198\t\t\tbMr = z+(indx[i]-1)*4*NoM+NoM;\n"
2004.07.14 03:28:00
<- "1199\t\t\tbMi = bMr+2*NoM;\n"
2004.07.14 03:28:00
<- "1200\t\t\tfor (j=0; j<NoM; ++j) {\n"
2004.07.14 03:28:00
<- "1201\t\t\t\tfor(k=0; k<N; ++k){\n"
2004.07.14 03:28:00
<- "1202\t\t\t\t\tvec[k] += bMr[j]*BoxModes[j*N+k];\n"
2004.07.14 03:28:00
<- "1203\t\t\t\t\tvec[k+N] += bMi[j]*BoxModes[j*N+k];\n"
2004.07.14 03:28:00
<- "1204\t\t\t\t}\n"
2004.07.14 03:28:00
<- "1205\t\t\t}\n"
2004.07.14 03:28:00
<- "1206\t\t\tfor(k=0; k<20; ++k) buff[k] = 0;\n"
2004.07.14 03:28:00
<- "1207\t\t\tsprintf(buff, \"v%d.dat\", i);\n"
2004.07.14 03:28:00
<- "1208\t\t\tout1 = fopen(buff, \"wt\");\n"
2004.07.14 03:28:00
<- "1209\t\t\tfor(k=0; k<N; ++k){\n"
2004.07.14 03:28:00
<- "1210\t\t\t\tx=xL+k*dx;\n"
2004.07.14 03:28:00
<- "1211\t\t\t\ttmp = vec[k]*vec[k]+vec[k+N]*vec[k+N];\n"
2004.07.14 03:28:00
<- "1212\t\t\t\tfprintf(out1, \"%lf\\t%.12lf\\t%.12lf\\t%.12lf\\n"
   "\", x, vec[k], vec[k+N], tmp);\n"
2004.07.14 03:28:00
<- "1213\t\t\t}\n"
2004.07.14 03:28:00
<- "1214\t\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1215\t\t\t\n"
2004.07.14 03:28:00
<- "1216\t\t}\n"
2004.07.14 03:28:00
<- "1217\t\n"
2004.07.14 03:28:00
<- "1218\t\tfree(vec);\n"
2004.07.14 03:28:00
<- "1219\t}\n"
2004.07.14 03:28:00
<- "1220\t\n"
2004.07.14 03:28:00
<- "1221\tint family(int nr){\n"
2004.07.14 03:28:00
<- "1222\t\tint j;\t\n"
2004.07.14 03:28:00
<- "1223\t\tdouble *ur, *ui, *vr, *vi;\n"
2004.07.14 03:28:00
<- "1224\t\tdouble normu, normv;\n"
2004.07.14 03:28:00
<- "1225\t\tdouble tmp, tmp1;\n"
2004.07.14 03:28:00
<- "1226\t\n"
2004.07.14 03:28:00
<- "1227\t\tur = z+nr*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1228\t\tvr = ur+NoM;\n"
2004.07.14 03:28:00
<- "1229\t\tui = vr+NoM;\n"
2004.07.14 03:28:00
<- "1230\t\tvi = ui+NoM;\n"
2004.07.14 03:28:00
<- "1231\t\t\n"
2004.07.14 03:28:00
<- "1232\t\tnormu = 0.0;\n"
2004.07.14 03:28:00
<- "1233\t\tnormv = 0.0;\n"
2004.07.14 03:28:00
<- "1234\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "1235\t\t\tnormu += ur[j]*ur[j]+ui[j]*ui[j];\n"
2004.07.14 03:28:00
<- "1236\t\t\tnormv += vr[j]*vr[j]+vi[j]*vi[j];\n"
2004.07.14 03:28:00
<- "1237\t\t}\n"
2004.07.14 03:28:00
<- "1238\t\ttmp = normu-normv;\n"
2004.07.14 03:28:00
<- "1239\t\t\n"
2004.07.14 03:28:00
<- "1240\t\tif(fabs(tmp) < 1e-12) {\n"
2004.07.14 03:28:00
<- "1241\t\t\treturn 0;\t\n"
2004.07.14 03:28:00
<- "1242\t\t}\n"
2004.07.14 03:28:00
<- "1243\t\telse{\n"
2004.07.14 03:28:00
<- "1244\t\t\ttmp1 = 1.0/sqrt(fabs(tmp));\n"
2004.07.14 03:28:00
<- "1245\t\t\tfor(j=0; j<NoM; j++){\n"
2004.07.14 03:28:00
<- "1246\t\t\t\tur[j]*=tmp1;\n"
2004.07.14 03:28:00
<- "1247\t\t\t\tui[j]*=tmp1;\n"
2004.07.14 03:28:00
<- "1248\t\t\t\tvr[j]*=tmp1;\n"
2004.07.14 03:28:00
<- "1249\t\t\t\tvi[j]*=tmp1;\n"
2004.07.14 03:28:00
<- "1250\t\t\t}\n"
2004.07.14 03:28:00
<- "1251\t\t\tif(tmp < 0) {\n"
2004.07.14 03:28:00
<- "1252\t\t\t\treturn -1;\t\n"
2004.07.14 03:28:00
<- "1253\t\t\t}\n"
2004.07.14 03:28:00
<- "1254\t\t\tif(tmp > 0) {\n"
2004.07.14 03:28:00
<- "1255\t\t\t\treturn 1;\t\n"
2004.07.14 03:28:00
<- "1256\t\t\t}\n"
2004.07.14 03:28:00
<- "1257\t\t}\n"
2004.07.14 03:28:00
<- "1258\t}\n"
2004.07.14 03:28:00
<- "1259\t\n"
2004.07.14 03:28:00
<- "1260\tvoid checkForSymmetry(void){\n"
2004.07.14 03:28:00
<- "1261\t\tint i, j;\n"
2004.07.14 03:28:00
<- "1262\t\tdouble sym, asym1, asym2;\n"
2004.07.14 03:28:00
<- "1263\t\n"
2004.07.14 03:28:00
<- "1264\t\tsym=0.0;\n"
2004.07.14 03:28:00
<- "1265\t\tasym1=0.0;\n"
2004.07.14 03:28:00
<- "1266\t\tfor(i=0; i<UVSize; ++i){\n"
2004.07.14 03:28:00
<- "1267\t\t\tfor(j=i+1; j<UVSize; ++j){\n"
2004.07.14 03:28:00
<- "1268\t\t\t\tsym += fabs(ZM[i*2*UVSize+j]-ZM[j*2*UVSize +i]);\n"
2004.07.14 03:28:00
<- "1269\t\t\t\tasym1 += fabs(ZM[i*2*UVSize+j+UVSize]-ZM[j*2*UVSize +i+UVSize]);\n"
2004.07.14 03:28:00
<- "1270\t\t\t}\n"
2004.07.14 03:28:00
<- "1271\t\t}\n"
2004.07.14 03:28:00
<- "1272\t\tasym2=0.0;\n"
2004.07.14 03:28:00
<- "1273\t\tfor(i=0; i<UVSize; ++i){\n"
2004.07.14 03:28:00
<- "1274\t\t\tasym2 += fabs(ZM[i*2*UVSize+i+UVSize]);\n"
2004.07.14 03:28:00
<- "1275\t\t}\n"
2004.07.14 03:28:00
<- "1276\t\n"
2004.07.14 03:28:00
<- "1277\t\tprintf(\"Symmetry of ZM: %.12lf\\t%.12lf\\t%.12lf\\n"
   "\", sym, asym1, asym2);\n"
2004.07.14 03:28:00
<- "1278\t\t\n"
2004.07.14 03:28:00
<- "1279\t}\n"
2004.07.14 03:28:00
<- "1280\t\n"
2004.07.14 03:28:00
<- "1281\tvoid reorganizeEigenvectors(void){\n"
2004.07.14 03:28:00
<- "1282\t\tint i,j;\n"
2004.07.14 03:28:00
<- "1283\t\tdouble *ur1, *ui1, *vr1, *vi1;\n"
2004.07.14 03:28:00
<- "1284\t\tdouble *ur2, *ui2, *vr2, *vi2;\n"
2004.07.14 03:28:00
<- "1285\t\tdouble *tmp1, *tmp2;\n"
2004.07.14 03:28:00
<- "1286\t\tdouble eps=1e-12;\n"
2004.07.14 03:28:00
<- "1287\t\n"
2004.07.14 03:28:00
<- "1288\t\ttmp1 = (double *)malloc((size_t)2*SizeOfBase*sizeof(double));\n"
2004.07.14 03:28:00
<- "1289\t\ttmp2 = (double *)malloc((size_t)2*SizeOfBase*sizeof(double));\n"
2004.07.14 03:28:00
<- "1290\t\n"
2004.07.14 03:28:00
<- "1291\t\tfor(i=0; i<2*SizeOfBase; ++i){\n"
2004.07.14 03:28:00
<- "1292\t\t\tif(fabs(wr[i]-wr[i+1])<eps && wi[i] && fabs(wi[i]+wi[i+1])<eps){\n"
2004.07.14 03:28:00
<- "1293\t\t\t\tur1 = z+i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1294\t\t\t\tvr1 = ur1+NoM;\n"
2004.07.14 03:28:00
<- "1295\t\t\t\tui1 = vr1+NoM;\n"
2004.07.14 03:28:00
<- "1296\t\t\t\tvi1 = ui1+NoM;\n"
2004.07.14 03:28:00
<- "1297\t\n"
2004.07.14 03:28:00
<- "1298\t\t\t\tur2 = vi1+NoM;\n"
2004.07.14 03:28:00
<- "1299\t\t\t\tvr2 = ur2+NoM;\n"
2004.07.14 03:28:00
<- "1300\t\t\t\tui2 = vr2+NoM;\n"
2004.07.14 03:28:00
<- "1301\t\t\t\tvi2 = ui2+NoM;\n"
2004.07.14 03:28:00
<- "1302\t\n"
2004.07.14 03:28:00
<- "1303\t\t\t\tfor(j=0; j<NoM; ++j){\n"
2004.07.14 03:28:00
<- "1304\t\t\t\t\ttmp1[j] = ur1[j]-ui2[j];\n"
2004.07.14 03:28:00
<- "1305\t\t\t\t\ttmp1[j+NoM] = vr1[j]-vi2[j];\n"
2004.07.14 03:28:00
<- "1306\t\t\t\t\ttmp1[j+2*NoM] = ui1[j]+ur2[j];\n"
2004.07.14 03:28:00
<- "1307\t\t\t\t\ttmp1[j+3*NoM] = vi1[j]+vr2[j];\n"
2004.07.14 03:28:00
<- "1308\t\t\t\t\t\n"
2004.07.14 03:28:00
<- "1309\t\t\t\t\ttmp2[j] = ur1[j]+ui2[j];\n"
2004.07.14 03:28:00
<- "1310\t\t\t\t\ttmp2[j+NoM] = vr1[j]+vi2[j];\n"
2004.07.14 03:28:00
<- "1311\t\t\t\t\ttmp2[j+2*NoM] = ui1[j]-ur2[j];\n"
2004.07.14 03:28:00
<- "1312\t\t\t\t\ttmp2[j+3*NoM] = vi1[j]-vr2[j];\n"
2004.07.14 03:28:00
<- "1313\t\t\t\t}\n"
2004.07.14 03:28:00
<- "1314\t\t\t\tfor(j=0; j<2*SizeOfBase; ++j){\n"
2004.07.14 03:28:00
<- "1315\t\t\t\t\tz[i*2*SizeOfBase +j] = tmp1[j];\n"
2004.07.14 03:28:00
<- "1316\t\t\t\t\tz[(i+1)*2*SizeOfBase +j] = tmp2[j];\n"
2004.07.14 03:28:00
<- "1317\t\t\t\t}\n"
2004.07.14 03:28:00
<- "1318\t\t\t\ti++;\n"
2004.07.14 03:28:00
<- "1319\t\t\t}\n"
2004.07.14 03:28:00
<- "1320\t\t}\n"
2004.07.14 03:28:00
<- "1321\t\n"
2004.07.14 03:28:00
<- "1322\t\tfree(tmp1);\n"
2004.07.14 03:28:00
<- "1323\t\tfree(tmp2);\n"
2004.07.14 03:28:00
<- "1324\t}\n"
2004.07.14 03:28:00
<- "1325\t\n"
2004.07.14 03:28:00
<- "1326\tvoid checkDum60(void){\n"
2004.07.14 03:28:00
<- "1327\t\tint i,j,k;\n"
2004.07.14 03:28:00
<- "1328\t\tdouble tmpr, tmpi;\n"
2004.07.14 03:28:00
<- "1329\t\tdouble *ur1, *ui1, *vr1, *vi1;\n"
2004.07.14 03:28:00
<- "1330\t\tdouble *ur2, *ui2, *vr2, *vi2;\n"
2004.07.14 03:28:00
<- "1331\t\n"
2004.07.14 03:28:00
<- "1332\t\tout1 = fopen(\"dum60a.dat\",\"wt\");\n"
2004.07.14 03:28:00
<- "1333\t\tfor(i=0; i<UVSize; ++i){\n"
2004.07.14 03:28:00
<- "1334\t\t\tur1 = uv+i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1335\t//\tfor(i=0; i<SizeOfBase; ++i){\n"
2004.07.14 03:28:00
<- "1336\t//\t\tur1 = z+i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1337\t\t\tvr1 = ur1 + NoM;\n"
2004.07.14 03:28:00
<- "1338\t\t\tui1 = vr1 + NoM;\n"
2004.07.14 03:28:00
<- "1339\t\t\tvi1 = ui1 + NoM;\n"
2004.07.14 03:28:00
<- "1340\t\t\tfor(j=i; j<UVSize; ++j){\n"
2004.07.14 03:28:00
<- "1341\t\t\t\tur2 = uv+j*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1342\t//\t\tfor(j=i; j<SizeOfBase; ++j){\n"
2004.07.14 03:28:00
<- "1343\t//\t\t\tur2 = z+j*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1344\t\t\t\tvr2 = ur2 + NoM;\n"
2004.07.14 03:28:00
<- "1345\t\t\t\tui2 = vr2 + NoM;\n"
2004.07.14 03:28:00
<- "1346\t\t\t\tvi2 = ui2 + NoM;\n"
2004.07.14 03:28:00
<- "1347\t\t\t\ttmpr=tmpi=0.0;\n"
2004.07.14 03:28:00
<- "1348\t\t\t\tfor(k=0; k<NoM; ++k){\n"
2004.07.14 03:28:00
<- "1349\t\t\t\t\ttmpr += ur1[k]*ur2[k] + ui1[k]*ui2[k] - vr1[k]*vr2[k] - vi1[k]*vi2[k]; \n"
2004.07.14 03:28:00
<- "1350\t\t\t\t\ttmpi += ur1[k]*ui2[k] - ui1[k]*ur2[k] - vr1[k]*vi2[k] + vi1[k]*vr2[k]; \n"
2004.07.14 03:28:00
<- "1351\t\t\t\t}\n"
2004.07.14 03:28:00
<- "1352\t\t\t\tfprintf(out1, \"%d\\t%d\\t%.16lf\\t%.16lf\\t%.16lf\\n"
   "\", i, j, tmpr, tmpi, tmpr*tmpr+tmpi*tmpi);\n"
2004.07.14 03:28:00
<- "1353\t\t\t}\n"
2004.07.14 03:28:00
<- "1354\t\t}\n"
2004.07.14 03:28:00
<- "1355\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1356\t\n"
2004.07.14 03:28:00
<- "1357\t\tout1 = fopen(\"dum60b.dat\",\"wt\");\n"
2004.07.14 03:28:00
<- "1358\t\tfor(i=0; i<UVSize; ++i){\n"
2004.07.14 03:28:00
<- "1359\t\t\tur1 = uv+i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1360\t//\tfor(i=0; i<SizeOfBase; ++i){\n"
2004.07.14 03:28:00
<- "1361\t//\t\tur1 = z+i*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1362\t\t\tvr1 = ur1 + NoM;\n"
2004.07.14 03:28:00
<- "1363\t\t\tui1 = vr1 + NoM;\n"
2004.07.14 03:28:00
<- "1364\t\t\tvi1 = ui1 + NoM;\n"
2004.07.14 03:28:00
<- "1365\t\t\tfor(j=i; j<UVSize; ++j){\n"
2004.07.14 03:28:00
<- "1366\t\t\t\tur2 = uv+j*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1367\t//\t\tfor(j=i; j<SizeOfBase; ++j){\n"
2004.07.14 03:28:00
<- "1368\t//\t\t\tur2 = z+j*2*SizeOfBase;\n"
2004.07.14 03:28:00
<- "1369\t\t\t\tvr2 = ur2 + NoM;\n"
2004.07.14 03:28:00
<- "1370\t\t\t\tui2 = vr2 + NoM;\n"
2004.07.14 03:28:00
<- "1371\t\t\t\tvi2 = ui2 + NoM;\n"
2004.07.14 03:28:00
<- "1372\t\t\t\ttmpr=tmpi=0.0;\n"
2004.07.14 03:28:00
<- "1373\t\t\t\tfor(k=0; k<NoM; ++k){\n"
2004.07.14 03:28:00
<- "1374\t\t\t\t\ttmpr += vr1[k]*ur2[k] - vi1[k]*ui2[k] - ur1[k]*vr2[k] + ui1[k]*vi2[k]; \n"
2004.07.14 03:28:00
<- "1375\t\t\t\t\ttmpi += -vr1[k]*ui2[k] - vi1[k]*ur2[k] + ur1[k]*vi2[k] + ui1[k]*vr2[k]; \n"
2004.07.14 03:28:00
<- "1376\t\t\t\t}\n"
2004.07.14 03:28:00
<- "1377\t\t\t\tfprintf(out1, \"%d\\t%d\\t%.16lf\\t%.16lf\\t%.16lf\\n"
   "\", i, j, tmpr, tmpi, tmpr*tmpr+tmpi*tmpi);\n"
2004.07.14 03:28:00
<- "1378\t\t\t}\n"
2004.07.14 03:28:00
<- "1379\t\t}\n"
2004.07.14 03:28:00
<- "1380\t\tfclose(out1);\n"
2004.07.14 03:28:00
<- "1381\t}\n"
2004.07.14 03:28:00
<- "(gdb) "
#  Reading file "BogolubovBox.c" from GDB...done.
#  Source "BogolubovBox.c" (from GDB) 1381 lines, 44648 characters
2004.07.14 03:28:00
-> "info source\n"
2004.07.14 03:28:00
<- "Current source file is BogolubovBox.c\n"
2004.07.14 03:28:00
<- "Compilation directory is /home/zbyszek/programs/shocks/src\n"
2004.07.14 03:28:00
<- "Located in /home/zbyszek/programs/shocks/src/BogolubovBox.c\n"
2004.07.14 03:28:00
<- "Contains 1381 lines.\n"
2004.07.14 03:28:00
<- "Source language is c.\n"
2004.07.14 03:28:00
<- "Compiled with DWARF 2 debugging format.\n"
2004.07.14 03:28:00
<- "Does not include preprocessor macro info.\n"
2004.07.14 03:28:00
<- "(gdb) "
2004.07.14 03:28:02
-> "next\n"
2004.07.14 03:28:02
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:279:5402:beg:0x8048c44\n"
2004.07.14 03:28:02
<- "(gdb) "
2004.07.14 03:28:02
-> "info breakpoints\n"
2004.07.14 03:28:02
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:02
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:02
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:02
<- "(gdb) "
#  File "/home/zbyszek/programs/shocks/src/BogolubovBox.c" 1381 lines, 44648 characters
2004.07.14 03:28:03
-> "next\n"
2004.07.14 03:28:03
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:280:5423:beg:0x8048c50\n"
2004.07.14 03:28:03
<- "(gdb) "
2004.07.14 03:28:03
-> "info breakpoints\n"
2004.07.14 03:28:03
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:03
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:03
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:03
<- "(gdb) "
2004.07.14 03:28:03
-> "next\n"
2004.07.14 03:28:03
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:282:5478:beg:0x8048c65\n"
2004.07.14 03:28:03
<- "(gdb) "
2004.07.14 03:28:03
-> "info breakpoints\n"
2004.07.14 03:28:03
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:03
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:03
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:03
<- "(gdb) "
2004.07.14 03:28:04
-> "next\n"
2004.07.14 03:28:04
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:283:5506:beg:0x8048c7e\n"
2004.07.14 03:28:04
<- "(gdb) "
2004.07.14 03:28:04
-> "info breakpoints\n"
2004.07.14 03:28:04
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:04
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:04
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:04
<- "(gdb) "
2004.07.14 03:28:04
-> "next\n"
2004.07.14 03:28:04
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:287:5607:beg:0x8048cd5\n"
2004.07.14 03:28:04
<- "(gdb) "
2004.07.14 03:28:04
-> "info breakpoints\n"
2004.07.14 03:28:04
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:04
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:04
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:04
<- "(gdb) "
2004.07.14 03:28:05
-> "next\n"
2004.07.14 03:28:05
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:291:5704:beg:0x8048d20\n"
2004.07.14 03:28:05
<- "(gdb) "
2004.07.14 03:28:05
-> "info breakpoints\n"
2004.07.14 03:28:05
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:05
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:05
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:05
<- "(gdb) "
2004.07.14 03:28:05
-> "next\n"
2004.07.14 03:28:05
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:295:5802:beg:0x8048d6a\n"
2004.07.14 03:28:05
<- "(gdb) "
2004.07.14 03:28:05
-> "info breakpoints\n"
2004.07.14 03:28:05
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:05
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:05
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:05
<- "(gdb) "
2004.07.14 03:28:05
-> "next\n"
2004.07.14 03:28:05
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:297:5817:beg:0x8048d77\n"
2004.07.14 03:28:05
<- "(gdb) "
2004.07.14 03:28:05
-> "info breakpoints\n"
2004.07.14 03:28:05
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:05
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:05
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:05
<- "(gdb) "
2004.07.14 03:28:05
-> "next\n"
2004.07.14 03:28:05
<- "GP energy = 52.820173\n"
2004.07.14 03:28:05
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:299:5852:beg:0x8048d8d\n"
2004.07.14 03:28:05
<- "(gdb) "
2004.07.14 03:28:05
-> "info breakpoints\n"
2004.07.14 03:28:05
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:05
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:05
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:05
<- "(gdb) "
2004.07.14 03:28:06
-> "next\n"
2004.07.14 03:28:06
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:300:5862:beg:0x8048d92\n"
2004.07.14 03:28:06
<- "(gdb) "
2004.07.14 03:28:06
-> "info breakpoints\n"
2004.07.14 03:28:06
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:06
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:06
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:06
<- "(gdb) "
2004.07.14 03:28:06
-> "next\n"
2004.07.14 03:28:06
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:301:5885:beg:0x8048da7\n"
2004.07.14 03:28:06
<- "(gdb) "
2004.07.14 03:28:06
-> "info breakpoints\n"
2004.07.14 03:28:06
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:06
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:06
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:06
<- "(gdb) "
2004.07.14 03:28:06
-> "next\n"
2004.07.14 03:28:06
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:300:5862:beg:0x8048dd4\n"
2004.07.14 03:28:06
<- "(gdb) "
2004.07.14 03:28:06
-> "info breakpoints\n"
2004.07.14 03:28:06
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:06
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:06
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:06
<- "(gdb) "
2004.07.14 03:28:06
-> "next\n"
2004.07.14 03:28:06
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:301:5885:beg:0x8048da7\n"
2004.07.14 03:28:06
<- "(gdb) "
2004.07.14 03:28:06
-> "info breakpoints\n"
2004.07.14 03:28:06
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:06
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:06
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:06
<- "(gdb) "
2004.07.14 03:28:07
-> "next\n"
2004.07.14 03:28:07
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:300:5862:beg:0x8048dd4\n"
2004.07.14 03:28:07
<- "(gdb) "
2004.07.14 03:28:07
-> "info breakpoints\n"
2004.07.14 03:28:07
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:07
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:07
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:07
<- "(gdb) "
2004.07.14 03:28:07
-> "next\n"
2004.07.14 03:28:07
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:301:5885:beg:0x8048da7\n"
2004.07.14 03:28:07
<- "(gdb) "
2004.07.14 03:28:07
-> "info breakpoints\n"
2004.07.14 03:28:07
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:07
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:07
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:07
<- "(gdb) "
2004.07.14 03:28:12
-> "until\n"
2004.07.14 03:28:12
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:300:5862:beg:0x8048dd4\n"
2004.07.14 03:28:12
<- "(gdb) "
2004.07.14 03:28:12
-> "info breakpoints\n"
2004.07.14 03:28:12
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:12
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:12
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:12
<- "(gdb) "
2004.07.14 03:28:13
-> "until\n"
2004.07.14 03:28:15
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:303:5912:beg:0x8048ddb\n"
2004.07.14 03:28:15
<- "(gdb) "
2004.07.14 03:28:15
-> "info breakpoints\n"
2004.07.14 03:28:15
<- "Num Type           Disp Enb Address    What\n"
2004.07.14 03:28:15
<- "2   breakpoint     keep y   0x08048a05 in main at BogolubovBox.c:91\n"
2004.07.14 03:28:15
<- "\tbreakpoint already hit 1 time\n"
2004.07.14 03:28:15
<- "(gdb) "
2004.07.14 03:28:31
-> "frame\n"
2004.07.14 03:28:31
<- "#0  allocateAndFillUpArrays () at BogolubovBox.c:303\n"
2004.07.14 03:28:31
<- "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:303:5912:beg:0x8048ddb\n"
2004.07.14 03:28:31
<- "(gdb) "
2004.07.14 03:28:31
-> "display phi\n"
2004.07.14 03:28:31
<- "(gdb) "
2004.07.14 03:28:31
-> "display\n"
2004.07.14 03:28:31
<- "1: phi = (double *) 0x4068c600\n"
2004.07.14 03:28:31
<- "(gdb) "
#  Creating display...
!  Nothing to plot.
2004.07.14 03:28:31
-> "undisplay 1\n"
#  Creating display...done.
2004.07.14 03:28:31
<- "(gdb) "
2004.07.14 03:28:43
-> "frame\n"
2004.07.14 03:28:43
<- "#0  allocateAndFillUpArrays () at BogolubovBox.c:303\n"
   "\032\032/home/zbyszek/programs/shocks/src/BogolubovBox.c:303:5912:beg:0x8048ddb\n"
2004.07.14 03:28:43
<- "(gdb) "
2004.07.14 03:28:43
-> "display address@hidden"
2004.07.14 03:28:43
<- "(gdb) "
2004.07.14 03:28:43
-> "display\n"
2004.07.14 03:28:43
<- "2: phi[0] @ 20 = {0.063663936122987763, 0.063663936336536123, 0.063663936950764638, 0.063663937888723723, 0.063663939030006209, 0.06366394022316621, 0.063663941301954477, 0.063663942105007074, 0.06366394249764426, 0.063663942393125394, 0.063663941769833143, 0.063663940681879014, 0.063663939259716681, 0.063663937699201437, 0.063663936238511867, 0.063663935125107193, 0.063663934576178277, 0.063663934739237135, 0.063663935658942336, 0.06366393725738155}\n"
2004.07.14 03:28:43
<- "(gdb) "
#  Creating display...
#  Creating display...done.
>> "set term xlib\n"
   "set parametric\n"
   "set urange [0:1]\n"
   "set vrange [0:1]\n"
   "set trange [0:1]"
>> "set noborder\n"
   "plot \"/tmp/dddzCG9Xm\" title \"phi[0] @ 20\""
2004.07.14 03:28:43
-> "output &(phi[0] @ 20)\n"
#  Display 2: phi[0] @ 20 (enabled, scope allocateAndFillUpArrays)
2004.07.14 03:28:43
<- "(double (*)[20]) 0x4068c600(gdb) "
#  Display 2: phi[0] @ 20 (enabled, scope allocateAndFillUpArrays, address 0x4068c600)
<< "G0 29360142"
<< "QD"
<< "A90"
<< "A0"
<< "W0001"
<< "L-002"
<< "M05610326"
<< "V05880326"
<< "M39340326"
<< "V39070326"
<< "J0002"
<< "QF"
<< "T05210326 0.0636639"
<< "W0001"
<< "L-002"
<< "M05610736"
<< "V05880736"
<< "M39340736"
<< "V39070736"
<< "J0002"
<< "T05210736 0.0636639"
<< "W0001"
<< "L-002"
<< "M05611146"
<< "V05881146"
<< "M39341146"
<< "V39071146"
<< "J0002"
<< "T05211146 0.0636639"
<< "W0001"
<< "L-002"
<< "M05611556"
<< "V05881556"
<< "M39341556"
<< "V39071556"
<< "J0002"
<< "T05211556 0.0636639"
<< "W0001"
<< "L-002"
<< "M05611966"
<< "V05881966"
<< "M39341966"
<< "V39071966"
<< "J0002"
<< "T05211966 0.0636639"
<< "W0001"
<< "L-002"
<< "M05612376"
<< "V05882376"
<< "M39342376"
<< "V39072376"
<< "J0002"
<< "T05212376 0.0636639"
<< "W0001"
<< "L-002"
<< "M05612786"
<< "V05882786"
<< "M39342786"
<< "V39072786"
<< "J0002"
<< "T05212786 0.0636639"
<< "W0001"
<< "L-002"
<< "M05613196"
<< "V05883196"
<< "M39343196"
<< "V39073196"
<< "J0002"
<< "T05213196 0.0636639"
<< "W0001"
<< "L-002"
<< "M05613606"
<< "V05883606"
<< "M39343606"
<< "V39073606"
<< "J0002"
<< "T05213606 0.0636639"
<< "W0001"
<< "L-002"
<< "M05614016"
<< "V05884016"
<< "M39344016"
<< "V39074016"
<< "J0002"
<< "T05214016 0.0636639"
<< "W0001"
<< "L-002"
<< "A0"
<< "M05610326"
<< "V05610366"
<< "M05614016"
<< "V05613976"
<< "J0001"
<< "T05610163 0"
<< "W0001"
<< "L-002"
<< "M08980326"
<< "V08980366"
<< "M08984016"
<< "V08983976"
<< "J0001"
<< "T08980163 2"
<< "W0001"
<< "L-002"
<< "M12360326"
<< "V12360366"
<< "M12364016"
<< "V12363976"
<< "J0001"
<< "T12360163 4"
<< "W0001"
<< "L-002"
<< "M15730326"
<< "V15730366"
<< "M15734016"
<< "V15733976"
<< "J0001"
<< "T15730163 6"
<< "W0001"
<< "L-002"
<< "M19100326"
<< "V19100366"
<< "M19104016"
<< "V19103976"
<< "J0001"
<< "T19100163 8"
<< "W0001"
<< "L-002"
<< "M22480326"
<< "V22480366"
<< "M22484016"
<< "V22483976"
<< "J0001"
<< "T22480163 10"
<< "W0001"
<< "L-002"
<< "M25850326"
<< "V25850366"
<< "M25854016"
<< "V25853976"
<< "J0001"
<< "T25850163 12"
<< "W0001"
<< "L-002"
<< "M29220326"
<< "V29220366"
<< "M29224016"
<< "V29223976"
<< "J0001"
<< "T29220163 14"
<< "W0001"
<< "L-002"
<< "M32590326"
<< "V32590366"
<< "M32594016"
<< "V32593976"
<< "J0001"
<< "T32590163 16"
<< "W0001"
<< "L-002"
<< "M35970326"
<< "V35970366"
<< "M35974016"
<< "V35973976"
<< "J0001"
<< "T35970163 18"
<< "W0001"
<< "L-002"
<< "M39340326"
<< "V39340366"
<< "M39344016"
<< "V39343976"
<< "J0001"
<< "T39340163 20"
<< "W0001"
<< "L-002"
<< "A0"
<< "P-2 13 20"
<< "P-2 13 20"
<< "W0001"
<< "L0000"
<< "L-002"
<< "J0002"
<< "T36273895phi[0] @ 20"
<< "L0000"
<< "P0 561 1196"
<< "P0 730 1284"
<< "P0 898 1536"
<< "P0 1067 1920"
<< "P0 1236 2388"
<< "P0 1404 2877"
<< "P0 1573 3320"
<< "P0 1742 3649"
<< "P0 1910 3810"
<< "P0 2079 3767"
<< "P0 2248 3512"
<< "P0 2416 3066"
<< "P0 2585 2482"
<< "P0 2753 1843"
<< "P0 2922 1244"
<< "P0 3091 787"
<< "P0 3259 562"
<< "P0 3428 629"
<< "P0 3597 1006"
<< "P0 3765 1662"
<< "P0 3760 3895"
<< "P-2 13 20"
<< "E"
2004.07.14 03:28:45
-> "output fopen\n"
2004.07.14 03:28:45
<- "{<text variable, no debug info>} 0x4060d990 <fopen>"
2004.07.14 03:28:45
<- "(gdb) "
2004.07.14 03:28:47
-> "output fread\n"
2004.07.14 03:28:47
<- "{<text variable, no debug info>} 0x4057c7c0 <fread>"
2004.07.14 03:28:47
<- "(gdb) "
2004.07.14 03:29:01
-> "output gp.bin\n"
2004.07.14 03:29:01
<- "No symbol \"gp\" in current context."
2004.07.14 03:29:01
<- "\n"
2004.07.14 03:29:01
<- "(gdb) "
!  X error
X Error of failed request:  BadWindow (invalid Window parameter)
  Major opcode of failed request:  10 (X_UnmapWindow)
  Resource id in failed request:  0x0
  Serial number of failed request:  43008
  Current serial number in output stream:  43011
!  X error
#  Saving options in "/home/zbyszek/.ddd/init"...
#  Searching "ddd/Ddd"...
#  Trying "/home/zbyszek/.ddd/ddd/Ddd"
#  Trying "/usr/share/ddd-3.3.9/ddd/Ddd"
#  Searching "ddd/Ddd"..."/usr/share/ddd-3.3.9/ddd/Ddd".
#  Saving options in "/home/zbyszek/.ddd/init"...done.
#  Saving history in "/home/zbyszek/.ddd/history"...
#  Saving history in "/home/zbyszek/.ddd/history"...done.
#  Thanks for using DDD 3.3.9!  (We apologize for the inconvenience.)
X Error of failed request:  BadWindow (invalid Window parameter)
  Major opcode of failed request:  10 (X_UnmapWindow)
  Resource id in failed request:  0x0
  Serial number of failed request:  43010
  Current serial number in output stream:  43355