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.openmm;
39  
40  import com.sun.jna.ptr.IntByReference;
41  import com.sun.jna.ptr.PointerByReference;
42  
43  import java.nio.IntBuffer;
44  
45  import static edu.uiowa.jopenmm.OpenMMLibrary.OpenMM_Discrete3DFunction_create;
46  import static edu.uiowa.jopenmm.OpenMMLibrary.OpenMM_Discrete3DFunction_destroy;
47  import static edu.uiowa.jopenmm.OpenMMLibrary.OpenMM_Discrete3DFunction_getFunctionParameters;
48  import static edu.uiowa.jopenmm.OpenMMLibrary.OpenMM_Discrete3DFunction_setFunctionParameters;
49  
50  /**
51   * This is a TabulatedFunction that computes a discrete three dimensional function f(x,y,z).
52   * To evaluate it, x, y, and z are each rounded to the nearest integer and the table element with those
53   * indices is returned. If any index is outside the range [0, size), the result is undefined.
54   */
55  public class Discrete3DFunction extends TabulatedFunction {
56  
57    /**
58     * Create a Discrete3DFunction f(x,y,z) based on a set of tabulated values.
59     *
60     * @param xsize  the number of table elements along the x direction
61     * @param ysize  the number of table elements along the y direction
62     * @param zsize  the number of table elements along the z direction
63     * @param values the tabulated values of the function f(x,y,z), ordered so that
64     *               values[i+xsize*j+xsize*ysize*k] = f(i,j,k). This must be of length xsize*ysize*zsize.
65     */
66    public Discrete3DFunction(int xsize, int ysize, int zsize, PointerByReference values) {
67      super(OpenMM_Discrete3DFunction_create(xsize, ysize, zsize, values));
68    }
69  
70    /**
71     * Destroy the discrete 3D function.
72     */
73    @Override
74    public void destroy() {
75      if (pointer != null) {
76        OpenMM_Discrete3DFunction_destroy(pointer);
77        pointer = null;
78      }
79    }
80  
81    /**
82     * Get the parameters for the tabulated function.
83     *
84     * @param xsize  the number of table elements along the x direction
85     * @param ysize  the number of table elements along the y direction
86     * @param zsize  the number of table elements along the z direction
87     * @param values the tabulated values of the function f(x,y,z), ordered so that
88     *               values[i+xsize*j+xsize*ysize*k] = f(i,j,k). This must be of length xsize*ysize*zsize.
89     */
90    public void getFunctionParameters(IntByReference xsize, IntByReference ysize, IntByReference zsize, PointerByReference values) {
91      OpenMM_Discrete3DFunction_getFunctionParameters(pointer, xsize, ysize, zsize, values);
92    }
93  
94    /**
95     * Get the parameters for the tabulated function.
96     *
97     * @param xsize  the number of table elements along the x direction
98     * @param ysize  the number of table elements along the y direction
99     * @param zsize  the number of table elements along the z direction
100    * @param values the tabulated values of the function f(x,y,z), ordered so that
101    *               values[i+xsize*j+xsize*ysize*k] = f(i,j,k). This must be of length xsize*ysize*zsize.
102    */
103   public void getFunctionParameters(IntBuffer xsize, IntBuffer ysize, IntBuffer zsize, PointerByReference values) {
104     OpenMM_Discrete3DFunction_getFunctionParameters(pointer, xsize, ysize, zsize, values);
105   }
106 
107   /**
108    * Set the parameters for the tabulated function.
109    *
110    * @param xsize  the number of table elements along the x direction
111    * @param ysize  the number of table elements along the y direction
112    * @param zsize  the number of table elements along the z direction
113    * @param values the tabulated values of the function f(x,y,z), ordered so that
114    *               values[i+xsize*j+xsize*ysize*k] = f(i,j,k). This must be of length xsize*ysize*zsize.
115    */
116   public void setFunctionParameters(int xsize, int ysize, int zsize, PointerByReference values) {
117     OpenMM_Discrete3DFunction_setFunctionParameters(pointer, xsize, ysize, zsize, values);
118   }
119 
120 }