Single Halo

For this test, we take the grid from the null result test, and add a blob around the center of one of the tiles. The code used to create the halo files is in single_halo.c. Again, you'll need to modify readfiles.c to read this and the velocity field file in.

I'm going to skip the diatribe about what the programs do, go back to the null result test if things are unclear.

Create Grid

cable:~/packages/voboz1.1.1 rpwagner$ gcc single_halo.c
cable:~/packages/voboz1.1.1 rpwagner$ a.out
cable:~/packages/voboz1.1.1 rpwagner$ ls -l single_halo*
-rw-r--r--   1 rpwagner  rpwagner  1132 Jan 14 13:52 single_halo.raw
-rw-r--r--   1 rpwagner  rpwagner  1132 Jan 14 13:52 single_halo_vel.raw
cable:~/packages/voboz1.1.1 rpwagner$ 

vozinit

The blob (halo) is inside of the 8 particles in the tile, so it should be well blocked from the guard particles.

cable:~/packages/voboz1.1.1 rpwagner$ ./vozinit single_halo.raw 0.3 1 2 goo

cable:~/packages/voboz1.1.1 rpwagner$ cat scrgoo 
#!/bin/bash -f
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 0 0 0
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 0 0 1
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 0 1 0
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 0 1 1
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 1 0 0
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 1 0 1
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 1 1 0
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 1 1 1
./voztie 2 goo
cable:~/packages/voboz1.1.1 rpwagner$ 

voz1b1 and voztie

cable:~/packages/voboz1.1.1 rpwagner$ ./scrgoo 
<< blah, blah, blah >>

We've got the results from each tile,

cable:~/packages/voboz1.1.1 rpwagner$ ls -1 part.goo.0*
part.goo.00.00.00
part.goo.00.00.01
part.goo.00.01.00
part.goo.00.01.01
part.goo.01.00.00
part.goo.01.00.01
part.goo.01.01.00
part.goo.01.01.01
cable:~/packages/voboz1.1.1 rpwagner$ 

and the data files from voztie.

cable:~/packages/voboz1.1.1 rpwagner$ ls -1 *goo.dat
adjgoo.dat
volgoo.dat
cable:~/packages/voboz1.1.1 rpwagner$ 

jovoz

cable:~/packages/voboz1.1.1 rpwagner$ ./jovoz adjgoo.dat volgoo.dat goozone.dat goozone.txt 0.01
About to jump ...
Post-jump ...
1 initial haloes found
Found zone adjacencies
Volumes range from 1.401526e-02 to 1.468750e+00.
Maxprob = 104.796516.
about to sort ...
cable:~/packages/voboz1.1.1 rpwagner$ 

Again, some files.

cable:~/packages/voboz1.1.1 rpwagner$ ls -1 goozone.*
goozone.dat
goozone.txt

And now for something new. The new clump creates a zone with a different volume.

cable:~/packages/voboz1.1.1 rpwagner$ head goozone.txt 
1
0       94      1.047965e+02    1.468750e+00
cable:~/packages/voboz1.1.1 rpwagner$}}}

=== {{{boz}}} ===

{{{boz}}} gives us even more new stuff, like finding a single halo.

{{{
cable:~/packages/voboz1.1.1 rpwagner$ ./boz 1.0 0.3 1.0 single_halo.raw single_halo_vel.raw goozone.dat goohalos.dat goohalos.txt 1.2
a = 1.000000 => Multiplying H0 by 1.000000e+00
cell, potconst = 0.007812,1.707547e-04
np = 94
np: 94, x: 0.000000,0.875000; y: 0.000000,0.875000; z: 0.000000,0.875000
vx: 0.000000,0.000000; vy: 0.000000,0.000000; vz: 0.000000,0.000000
np = 94, nhaloes = 1
Done! Outputting ...
Nhunb = 1
cable:~/packages/voboz1.1.1 rpwagner$
}}}

Here's the info on the halo found. Interestingly, all of the particles are included.

{{{
cable:~/packages/voboz1.1.1 rpwagner$ cat goohalos.txt 
1       1
0       94      94      1.047965e+02 1.401526e-02 -inf -inf     70      0       0       1 1 0.000000 0.000000 0.000000 0.000000e+00 0.000000e+00 0.000000e+00
cable:~/packages/voboz1.1.1 rpwagner$ 
}}}