tsi/TODO at master · msf/tsi · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
TSI TODO LIST!


DOING DSS CLEANUP.
    - make thread save:
        srchnod
        krige
        ksol


---------------------------------------------------------------
TO BE DONE:
- implement mask
    -DSS
        - srchnod()
        - krige() ?
        - main loop()
    -TSI
        Normalizar sistema de coordenadas 3d <-> index para DSS e TSI.
        .. problemas com sistema de coordenadas 3d -> index
        ..todo o lado
- convert lookup tables for x, y z to "index" table.
- significantly reduce lookup table sizes, implies:
    - search radius must be significantly smaller than grid size.
	- calculate search radius relative to x, y, and z
	- set lookup table size from search radius for x, y, z


GRID HEAP:
- remove grid_heap usage from cmgrid.
- remove covtable grids from grid_heap
- make a "dump_grid()" for BestAI and nextBAI, nextBCM at the end of any iteration.
- make a "free_grids()" in grid_heap to remove grids from memory -> seismic_grid

- forward "grid_heap" knowledge inside DSS.
- simplify grid_heap algorithm.

- reduce memory use in 3 GRIDS!!
	- preload random path in order.
	- keep _only_ a "chunk" (set of 50 million points..) of order grid in memory.
		- based on current position from order that is being used.
		- keep base position index for chunk[0]
		- create a get_order_pos(index)
	- resort corr grid and secondary data grid based on random path order.
		- keep only a chunk of corr grid and secondary grid.
	   	- base pointers for chunks
		- get_bai_pos(index)
		- get_corr_pos(index)

---------------------------------------------------------------
DONE
- remove bug that ocours when running in:
	- grid bigger than 33million points +
		linux 64 bits executable or
		any windows executable or
		linux 32 bits and optimizations other than: -03 -ffast-math -fomit-frame-pointer
- camada de abstracao para syscalls nao portaveis:
	time: gettimeofday()<
- printlayers()
- error reporting no "read_harddata_file()"
- mersenne twister como RNG
- removi tsi_io_mpi
- bug fix importante na DSS (bug nao esta presente na tsi-5.4 mas tava no "current")
- small fixes no tsi_io.c (melhorar bug reporting)
- pequenas mudanças de ficheiros para algumas funcoes.
- correcoes no tsi.conf