source: molecuilder/src/defs.hpp@ c75363

Last change on this file since c75363 was c75363, checked in by Frederik Heber <heber@…>, 17 years ago

HUGE REWRITE to allow for adaptive increase of the bond order, first working commit

Lots of code was thrown out:
-BottomUp, TopDown and GetAtomicFragments
-TEFactors are now used as "CreationCounters", i.e. the count how often a fragment as been created (ideal would be only once)
-ReduceToUniqueOnes and stuff all thrown out, since they are out-dated since use of hash table
Other changes:
-CreateListofUniqueFragments renamed to PowerSetGenerator
-PowerSetGenerator goes not over all reachable roots, but one given by calling function FragmentBOSSANOVA
-FragmentBOSSANOVA loops over all possible root sites and hands this over to PowerSetGenerator
-by the virtue of the hash table it is not important anymore whether created keysets are unique or not, as this is recognized in log(n). Hence, the label < label is not important anymore (and not applicable in an adaptive scheme with old, parsed keysets and unknown labels) (THIS IS HOWEVER NOT DONE YET!)
The setup is then as follows:

  1. FragmentMolecule
    1. parses adjacency, keysets and orderatsite files
    2. performs DFS to find connected subgraphs (to leave this in was a design decision: might be useful later)
    3. a RootStack is created for every subgraph (here, later we implement the "update 10 sites with highest energy contribution", and that's why this consciously not done in the following loop)
    4. in a loop over all subgraphs

d1. calls FragmentBOSSANOVA with this RootStack and within the subgraph molecule structure
d2. creates molecule (fragment)s from the returned keysets (StoreFragmentFromKeySet)

  1. combines the generated molecule lists from all subgraphs
  2. saves to disk: fragment configs, adjacency, orderatsite, keyset files
  1. FragmentBOSSANOVA
    1. constructs a complete keyset of the molecule
    2. In a loop over all possible roots from the given rootstack

b1. increases order of root site
b2. calls PowerSetGenerator with this order, the complete keyset and the rootkeynr
b3. for all consecutive lower levels PowerSetGenerator is called with the suborder, the higher order keyset as the restricted one and each site in the set as the root)
b4. these are merged into a fragment list of keysets

  1. All fragment lists (for all orders, i.e. from all destination fields) are merged into one list for return
  1. PowerSetGenerator
    1. initialises UniqueFragments structure
    2. fills edge list via BFS
    3. creates the fragment by calling recursive function SPFragmentGenerator with UniqueFragments structure, 0 as root distance, the edge set, its dimension and the current suborder
    4. Free'ing structure
  2. SPFragmentGenerator (nothing much has changed here)
    1. loops over every possible combination (2dimension of edge set)

a1. inserts current set, if there's still space left

a11. yes: calls SPFragmentGenerator with structure, created new edge list and size respective to root distance+1
a12. no: stores fragment into keyset list by calling InsertFragmentIntoGraph

a2. removes all items added into the snake stack (in UniqueFragments structure) added during level (root distance) and current set

  • Property mode set to 100644
File size: 2.3 KB
Line 
1/** \file defs.hpp
2 *
3 * Constant definitons and enumerations.
4 *
5 */
6
7#ifndef DEFS_HPP_
8#define DEFS_HPP_
9
10using namespace std;
11
12#define MYEPSILON 1e-13 //!< machine epsilon precision
13#define NDIM 3 //!< number of spatial dimensions
14#define MAX_ELEMENTS 128 //!< maximum number of elements for certain lookup tables
15#define AtomicLengthToAngstroem 0.52917721 //!< conversion factor from atomic length/bohrradius to angstroem
16#define BONDTHRESHOLD 0.5 //!< CSD threshold in bond check which is the width of the interval whose center is the sum of the covalent radii
17
18#define VERSIONSTRING "v1.0"
19
20//enum BondOrderScheme { NoScheme, BottomUp, TopDown, ANOVA, Combined }; //!< Fragmentation scheme used in BOSS
21
22enum EdgeType { Undetermined, TreeEdge, BackEdge }; //!< edge type in a graph after Depth-First-Search analysis.
23
24enum Shading { white, lightgray, darkgray, black }; //!< color in Breadth-First-Search analysis
25
26//enum CutCyclicBond { KeepBond, SaturateBond }; //!< Saturation scheme either atom- or bondwise
27
28// Specifting whether a value in the parameter file must be specified or is optional
29enum necessity { optional, //!< parameter is optional, if not given sensible value is chosen
30 critical //!< parameter must be given or programme won't initiate
31 };
32
33// Specifying the status of the on command line given config file
34enum ConfigStatus { absent, empty, present };
35
36// Specifying the cast type to be read of a parameter, see ParseForParameter()
37enum value_type { string_type, double_type, int_type, row_int, row_double, grid, lower_trigrid, upper_trigrid};
38
39// maximum length of any char array
40#define MAXSTRINGSIZE 255
41
42// various standard filenames
43#define DEFAULTCONFIG "main_pcp_linux"
44#define KEYSETFILE "KeySets.dat"
45#define ADJACENCYFILE "Adjacency.dat"
46#define ORDERATSITEFILE "OrderAtSite.dat"
47#define FRAGMENTPREFIX "BondFragment"
48#define STANDARDCONFIG "unknown.conf"
49#define STANDARDELEMENTSDB "elements.db"
50#define STANDARDVALENCEDB "valence.db"
51#define STANDARDORBITALDB "orbitals.db"
52#define STANDARDHBONDDISTANCEDB "Hbonddistance.db"
53#define STANDARDHBONDANGLEDB "Hbondangle.db"
54
55#define UPDATECOUNT 10 //!< update ten sites per BOSSANOVA interval
56
57#endif /*DEFS_HPP_*/
Note: See TracBrowser for help on using the repository browser.