forked from pygbe/pygbe
-
Notifications
You must be signed in to change notification settings - Fork 1
/
README_input_format
82 lines (53 loc) · 3.69 KB
/
README_input_format
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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
In this file, the format for the input files is detailed.
To run, PyGBe needs two input files: a "config" file and a "parameters" file, which are in the ./bem_pycuda/input_files directory.
Config file format:
------------------
This file has two parts:
FILE: indicates in that line a mesh file will be specified.
This code interfaces with the MSMS format, which generates a 'filename.vert' and 'filename.face' files.
Each 'filename' has the information for the mesh of one surface. To account for more than one surface (ie. for Stern layers, solvent filled cavities, several proteins), more than one FILE line is needed.
After 'filename', the user must specify what kind of surface 'filename' is. It can be:
- stern_layer: surface interfacer stern layer with solvent
- dielectric_interface: surface interfaces low dielectric (inside protein) and high dielectric (outside protein)
- internal_cavity: surface is an internal cavity (a result of '-all_components' flag in MSMS). This is important to specify because by default MSMS changes the vertex ordering for internal cavities.
- dirichlet_surface: surface of specified potential. The value of this potential is read from a text file which has to be specified next to 'dirichlet_surface'. (See .phi0 file in example)
- neumann_surface: surface of specified potential. The value of this potential is read from a text file which has to be specified next to 'neumann_surface'. (See .phi0 file in example)
FIELD: specifies physical parameters of the different regions.
Parameters:
- LorY: 1: Laplace,
2: Yukawa
- E?: 1: calculate the energy in this region
0: don't calculate energy in this region
Note: if region is surrounded by a dirichlet or neumann surface,
surface energy will be calculated.
- Dielec: Dielectric constant
- kappa: reciprocal of Debye length
- charges?: 0: No charges inside this region
1: There are charges in this region
- coulomb?: 0: don't calculate coulomb energy in this region
1: calculate coulomb energy in this region
- charge_file: location of the '.pqr' file with the location of the charges
- Nparent: Number of 'parent' surfaces (surface containing this region)
Of course, this is either 1 of 0 (if it corresponds to the infinite region)
- parent: file of the parent surface mesh, according to their position in the FILE list,
starting from 0 (eg. if the mesh file for the parent is the third one specified
in the FILE section, parent=2)
- Nchild: number of child surfaces (surfaces completely contained in this region).
- children: position of the mesh files for the children surface in the FILE section
Parameters file format:
----------------------
- Precision: double or float. (float not supported yet!).
- K: number of Gauss points per element (1, 3, 4, and 7 are supported).
- Nk: number of Gauss points per triangle edge for semi-analytical integration.
- threshold: defines region near singularity where semi-analytical technique is used.
if sqrt(2*Area)/r > threshold, integration is done semi-analytically.
- BSZ: CUDA block size.
- restart: number of iterations for GMRES to do restart.
- tolerance: GMRES tolerance.
- max_iter: maximum number of GMRES iterations.
- P: order of expansion in treecode.
- eps: epsilon machine.
- NCRIT: maximum number of boundary elements per twig box of tree structure.
- theta: multipole acceptance criterion of treecode.
- GPU: 0: don't use GPU.
1: use GPU.