aboutsummaryrefslogtreecommitdiff
path: root/src/macro/DXDB_guts.h
diff options
context:
space:
mode:
Diffstat (limited to 'src/macro/DXDB_guts.h')
-rw-r--r--src/macro/DXDB_guts.h39
1 files changed, 39 insertions, 0 deletions
diff --git a/src/macro/DXDB_guts.h b/src/macro/DXDB_guts.h
new file mode 100644
index 0000000..7d76ddc
--- /dev/null
+++ b/src/macro/DXDB_guts.h
@@ -0,0 +1,39 @@
+/*@@
+ @header DXDB_guts.h
+ @date Jun 98
+ @author Gabrielle Allen
+ @desc
+
+ Macro to calculate the first derivatives of the
+ shift with respect to x
+
+ The macro is defined in terms of standard variables in
+ @seefile DXDB_declare.h
+ @enddesc
+@@*/
+
+#ifndef DXDB_GUTS
+#define DXDB_GUTS
+
+#ifdef FCODE
+
+ DXDB_OO2DX = 1D0/(2D0*DXDB_DX)
+
+ DXDB_DXDBX = DXDB_OO2DX*(DXDB_BX_IP - DXDB_BX_IM)
+ DXDB_DXDBY = DXDB_OO2DX*(DXDB_BY_IP - DXDB_BY_IM)
+ DXDB_DXDBZ = DXDB_OO2DX*(DXDB_BZ_IP - DXDB_BZ_IM)
+
+#endif
+
+#ifdef CCODE
+
+ DXDB_OO2DX = 1/(2*GH->dx0);
+
+ DXDB_DXDBX = DXDB_OO2DX*(DXDB_BX_IP - DXDB_BX_IM);
+ DXDB_DXDBY = DXDB_OO2DX*(DXDB_BY_IP - DXDB_BY_IM);
+ DXDB_DXDBZ = DXDB_OO2DX*(DXDB_BZ_IP - DXDB_BZ_IM);
+
+#endif
+
+#endif
+