View Javadoc
1   // ******************************************************************************
2   //
3   // Title:       Force Field X.
4   // Description: Force Field X - Software for Molecular Biophysics.
5   // Copyright:   Copyright (c) Michael J. Schnieders 2001-2025.
6   //
7   // This file is part of Force Field X.
8   //
9   // Force Field X is free software; you can redistribute it and/or modify it
10  // under the terms of the GNU General Public License version 3 as published by
11  // the Free Software Foundation.
12  //
13  // Force Field X is distributed in the hope that it will be useful, but WITHOUT
14  // ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
15  // FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
16  // details.
17  //
18  // You should have received a copy of the GNU General Public License along with
19  // Force Field X; if not, write to the Free Software Foundation, Inc., 59 Temple
20  // Place, Suite 330, Boston, MA 02111-1307 USA
21  //
22  // Linking this library statically or dynamically with other modules is making a
23  // combined work based on this library. Thus, the terms and conditions of the
24  // GNU General Public License cover the whole combination.
25  //
26  // As a special exception, the copyright holders of this library give you
27  // permission to link this library with independent modules to produce an
28  // executable, regardless of the license terms of these independent modules, and
29  // to copy and distribute the resulting executable under terms of your choice,
30  // provided that you also meet, for each linked independent module, the terms
31  // and conditions of the license of that module. An independent module is a
32  // module which is not derived from or based on this library. If you modify this
33  // library, you may extend this exception to your version of the library, but
34  // you are not obligated to do so. If you do not wish to do so, delete this
35  // exception statement from your version.
36  //
37  // ******************************************************************************
38  package ffx.xray.solvent;
39  
40  import edu.rit.pj.ParallelTeam;
41  import ffx.crystal.Crystal;
42  import ffx.potential.bonded.Atom;
43  import ffx.potential.nonbonded.RowRegion;
44  
45  import java.util.logging.Level;
46  import java.util.logging.Logger;
47  
48  import static java.util.Arrays.fill;
49  
50  /**
51   * This class implements a spatial decomposition based on partitioning a grid into rows. The
52   * over-ridden "selectAtoms" method selects atoms that are not in the asymmetric unit, but are
53   * within the supplied cutoff radius.
54   *
55   * @author Michael J. Schnieders
56   * @since 1.0
57   */
58  public class BulkSolventRowRegion extends RowRegion {
59  
60    /** Constant <code>logger</code> */
61    protected static final Logger logger = Logger.getLogger(BulkSolventRowRegion.class.getName());
62  
63    private final BulkSolventList bulkSolventList;
64    private final int gZ;
65    private final int gY;
66  
67    /**
68     * Constructor for BulkSolventDensityRegion.
69     *
70     * @param gX a int.
71     * @param gY a int.
72     * @param gZ a int.
73     * @param grid an array of double.
74     * @param nSymm a int.
75     * @param threadCount a int.
76     * @param crystal a {@link ffx.crystal.Crystal} object.
77     * @param atoms an array of {@link ffx.potential.bonded.Atom} objects.
78     * @param coordinates an array of double.
79     * @param cutoff a double.
80     * @param parallelTeam a {@link edu.rit.pj.ParallelTeam} object.
81     */
82    public BulkSolventRowRegion(int gX, int gY, int gZ, double[] grid, int nSymm,
83        int threadCount, Crystal crystal, Atom[] atoms, double[][][] coordinates, double cutoff,
84        ParallelTeam parallelTeam) {
85      super(gX, gY, gZ, grid, nSymm, threadCount, atoms, coordinates);
86      this.gZ = gZ;
87      this.gY = gY;
88      // Asymmetric unit atoms never selected by this class.
89      fill(select[0], false);
90      bulkSolventList = new BulkSolventList(crystal, atoms, cutoff, parallelTeam);
91    }
92  
93    /** {@inheritDoc} */
94    @Override
95    public void run() {
96      try {
97        execute(0, (gZ * gY) - 1, rowLoop[getThreadIndex()]);
98      } catch (Exception e) {
99        String message = " Exception in BulkSolventRowRegion.";
100       logger.log(Level.SEVERE, message, e);
101     }
102   }
103 
104   /** {@inheritDoc} */
105   @Override
106   public void selectAtoms() {
107     bulkSolventList.buildList(coordinates, select, false);
108   }
109 }