diff --git a/examples/main.c b/examples/main.c
index 01365a419d24e659575f05eb6e73429a9e876d7f..8589798ade0c578a0e0559947db36eeed6437814 100644
--- a/examples/main.c
+++ b/examples/main.c
@@ -232,7 +232,7 @@ int main(int argc, char *argv[]) {
/* Initialize unit system */
struct UnitSystem us;
- initUnitSystem(&us, ¶ms);
+ units_init(&us, ¶ms);
if (myrank == 0) {
message("Unit system: U_M = %e g.", us.UnitMass_in_cgs);
message("Unit system: U_L = %e cm.", us.UnitLength_in_cgs);
diff --git a/src/common_io.c b/src/common_io.c
index 6183effe9ce392ab930c581cbd118f025bbce773..4a6dd1db5629552820dec0e6ef5204c059add603 100644
--- a/src/common_io.c
+++ b/src/common_io.c
@@ -282,15 +282,15 @@ void writeUnitSystem(hid_t h_file, struct UnitSystem* us) {
if (h_grpunit < 0) error("Error while creating Unit System group");
writeAttribute_d(h_grpunit, "Unit mass in cgs (U_M)",
- getBaseUnit(us, UNIT_MASS));
+ units_get_base_unit(us, UNIT_MASS));
writeAttribute_d(h_grpunit, "Unit length in cgs (U_L)",
- getBaseUnit(us, UNIT_LENGTH));
+ units_get_base_unit(us, UNIT_LENGTH));
writeAttribute_d(h_grpunit, "Unit time in cgs (U_t)",
- getBaseUnit(us, UNIT_TIME));
+ units_get_base_unit(us, UNIT_TIME));
writeAttribute_d(h_grpunit, "Unit current in cgs (U_I)",
- getBaseUnit(us, UNIT_CURRENT));
+ units_get_base_unit(us, UNIT_CURRENT));
writeAttribute_d(h_grpunit, "Unit temperature in cgs (U_T)",
- getBaseUnit(us, UNIT_TEMPERATURE));
+ units_get_base_unit(us, UNIT_TEMPERATURE));
H5Gclose(h_grpunit);
}
diff --git a/src/parallel_io.c b/src/parallel_io.c
index 0076c225e1c5361287280f8a567c8062aefd914e..e9d17d9c7dcc14b60bd721fae1a5b06a51140632 100644
--- a/src/parallel_io.c
+++ b/src/parallel_io.c
@@ -274,11 +274,11 @@ void writeArrayBackEnd(hid_t grp, char* fileName, FILE* xmfFile,
type);
/* Write unit conversion factors for this data set */
- conversionString(buffer, us, convFactor);
+ units_conversion_string(buffer, us, convFactor);
writeAttribute_d(h_data, "CGS conversion factor",
- conversionFactor(us, convFactor));
- writeAttribute_f(h_data, "h-scale exponent", hFactor(us, convFactor));
- writeAttribute_f(h_data, "a-scale exponent", aFactor(us, convFactor));
+ units_conversion_factor(us, convFactor));
+ writeAttribute_f(h_data, "h-scale exponent", units_h_factor(us, convFactor));
+ writeAttribute_f(h_data, "a-scale exponent", units_a_factor(us, convFactor));
writeAttribute_s(h_data, "Conversion factor", buffer);
/* Free and close everything */
diff --git a/src/serial_io.c b/src/serial_io.c
index 40bd2b1c8921f4acbfa0950984d6915ebd3d241e..b4a0269c58b56675c833f2ed4a3eea74042ee8ba 100644
--- a/src/serial_io.c
+++ b/src/serial_io.c
@@ -237,11 +237,11 @@ void prepareArray(hid_t grp, char* fileName, FILE* xmfFile,
writeXMFline(xmfFile, fileName, partTypeGroupName, name, N_total, dim, type);
/* Write unit conversion factors for this data set */
- conversionString(buffer, us, convFactor);
+ units_conversion_string(buffer, us, convFactor);
writeAttribute_d(h_data, "CGS conversion factor",
- conversionFactor(us, convFactor));
- writeAttribute_f(h_data, "h-scale exponent", hFactor(us, convFactor));
- writeAttribute_f(h_data, "a-scale exponent", aFactor(us, convFactor));
+ units_conversion_factor(us, convFactor));
+ writeAttribute_f(h_data, "h-scale exponent", units_h_factor(us, convFactor));
+ writeAttribute_f(h_data, "a-scale exponent", units_a_factor(us, convFactor));
writeAttribute_s(h_data, "Conversion factor", buffer);
H5Pclose(h_prop);
diff --git a/src/single_io.c b/src/single_io.c
index 801428433ef5170082b68dec425e52f845bb41ae..65f7e83edcf2fcca465c9b2f092e1ccfe7a55160 100644
--- a/src/single_io.c
+++ b/src/single_io.c
@@ -241,11 +241,11 @@ void writeArrayBackEnd(hid_t grp, char* fileName, FILE* xmfFile,
writeXMFline(xmfFile, fileName, partTypeGroupName, name, N, dim, type);
/* Write unit conversion factors for this data set */
- conversionString(buffer, us, convFactor);
+ units_conversion_string(buffer, us, convFactor);
writeAttribute_d(h_data, "CGS conversion factor",
- conversionFactor(us, convFactor));
- writeAttribute_f(h_data, "h-scale exponent", hFactor(us, convFactor));
- writeAttribute_f(h_data, "a-scale exponent", aFactor(us, convFactor));
+ units_conversion_factor(us, convFactor));
+ writeAttribute_f(h_data, "h-scale exponent", units_h_factor(us, convFactor));
+ writeAttribute_f(h_data, "a-scale exponent", units_a_factor(us, convFactor));
writeAttribute_s(h_data, "Conversion factor", buffer);
/* Free and close everything */
diff --git a/src/units.c b/src/units.c
index 46710b9d131b880057d13c31d42c4e0aaf3ec352..184dbe8a0df000008dba1d7003558d83b1f08cad 100644
--- a/src/units.c
+++ b/src/units.c
@@ -48,7 +48,7 @@
* @param us The UnitSystem to initialize.
* @param params The parsed parameter file.
*/
-void initUnitSystem(struct UnitSystem* us, const struct swift_params* params) {
+void units_init(struct UnitSystem* us, const struct swift_params* params) {
us->UnitMass_in_cgs =
parser_get_param_double(params, "UnitSystem:UnitMass_in_cgs");
@@ -68,7 +68,8 @@ void initUnitSystem(struct UnitSystem* us, const struct swift_params* params) {
* @param us The UnitSystem used
* @param baseUnit The base unit
*/
-double getBaseUnit(const struct UnitSystem* us, enum BaseUnits baseUnit) {
+double units_get_base_unit(const struct UnitSystem* us,
+ enum BaseUnits baseUnit) {
switch (baseUnit) {
case UNIT_MASS:
return us->UnitMass_in_cgs;
@@ -90,7 +91,7 @@ double getBaseUnit(const struct UnitSystem* us, enum BaseUnits baseUnit) {
* @brief Returns the base unit symbol
* @param baseUnit The base unit
*/
-const char* getBaseUnitSymbol(enum BaseUnits baseUnit) {
+const char* units_get_base_unit_symbol(enum BaseUnits baseUnit) {
switch (baseUnit) {
case UNIT_MASS:
return "U_M";
@@ -112,7 +113,7 @@ const char* getBaseUnitSymbol(enum BaseUnits baseUnit) {
* @brief Returns the base unit symbol in the cgs system
* @param baseUnit The base unit
*/
-const char* getBaseUnitCGSSymbol(enum BaseUnits baseUnit) {
+const char* units_get_base_unit_CGS_symbol(enum BaseUnits baseUnit) {
switch (baseUnit) {
case UNIT_MASS:
return "g";
@@ -130,8 +131,8 @@ const char* getBaseUnitCGSSymbol(enum BaseUnits baseUnit) {
return "";
}
-void getBaseUnitExponantsArray(float baseUnitsExp[5],
- enum UnitConversionFactor unit) {
+void units_get_base_unit_exponants_array(float baseUnitsExp[5],
+ enum UnitConversionFactor unit) {
switch (unit) {
case UNIT_CONV_NO_UNITS:
break;
@@ -272,13 +273,13 @@ void getBaseUnitExponantsArray(float baseUnitsExp[5],
* @param us The system of units in use
* @param unit The unit to convert
*/
-double conversionFactor(const struct UnitSystem* us,
- enum UnitConversionFactor unit) {
+double units_conversion_factor(const struct UnitSystem* us,
+ enum UnitConversionFactor unit) {
float baseUnitsExp[5] = {0.f};
- getBaseUnitExponantsArray(baseUnitsExp, unit);
+ units_get_base_unit_exponants_array(baseUnitsExp, unit);
- return generalConversionFactor(us, baseUnitsExp);
+ return units_general_conversion_factor(us, baseUnitsExp);
}
/**
@@ -286,12 +287,13 @@ double conversionFactor(const struct UnitSystem* us,
* @param us The system of units in use
* @param unit The unit to convert
*/
-float hFactor(const struct UnitSystem* us, enum UnitConversionFactor unit) {
+float units_h_factor(const struct UnitSystem* us,
+ enum UnitConversionFactor unit) {
float baseUnitsExp[5] = {0.f};
- getBaseUnitExponantsArray(baseUnitsExp, unit);
+ units_get_base_unit_exponants_array(baseUnitsExp, unit);
- return generalhFactor(us, baseUnitsExp);
+ return units_general_h_factor(us, baseUnitsExp);
}
/**
@@ -299,25 +301,26 @@ float hFactor(const struct UnitSystem* us, enum UnitConversionFactor unit) {
* @param us The system of units in use
* @param unit The unit to convert
*/
-float aFactor(const struct UnitSystem* us, enum UnitConversionFactor unit) {
+float units_a_factor(const struct UnitSystem* us,
+ enum UnitConversionFactor unit) {
float baseUnitsExp[5] = {0.f};
- getBaseUnitExponantsArray(baseUnitsExp, unit);
+ units_get_base_unit_exponants_array(baseUnitsExp, unit);
- return generalaFactor(us, baseUnitsExp);
+ return units_general_a_factor(us, baseUnitsExp);
}
/**
* @brief Returns a string containing the exponents of the base units making up
* the conversion factors
*/
-void conversionString(char* buffer, const struct UnitSystem* us,
- enum UnitConversionFactor unit) {
+void units_conversion_string(char* buffer, const struct UnitSystem* us,
+ enum UnitConversionFactor unit) {
float baseUnitsExp[5] = {0.f};
- getBaseUnitExponantsArray(baseUnitsExp, unit);
+ units_get_base_unit_exponants_array(baseUnitsExp, unit);
- generalConversionString(buffer, us, baseUnitsExp);
+ units_general_conversion_string(buffer, us, baseUnitsExp);
}
/**
@@ -327,14 +330,14 @@ void conversionString(char* buffer, const struct UnitSystem* us,
* @param baseUnitsExponants The exponent of each base units required to form
* the desired quantity. See conversionFactor() for a working example
*/
-double generalConversionFactor(const struct UnitSystem* us,
- float baseUnitsExponants[5]) {
+double units_general_conversion_factor(const struct UnitSystem* us,
+ float baseUnitsExponants[5]) {
double factor = 1.;
int i;
for (i = 0; i < 5; ++i)
if (baseUnitsExponants[i] != 0)
- factor *= pow(getBaseUnit(us, i), baseUnitsExponants[i]);
+ factor *= pow(units_get_base_unit(us, i), baseUnitsExponants[i]);
return factor;
}
@@ -345,7 +348,8 @@ double generalConversionFactor(const struct UnitSystem* us,
* @param baseUnitsExponants The exponent of each base units required to form
* the desired quantity. See conversionFactor() for a working example
*/
-float generalhFactor(const struct UnitSystem* us, float baseUnitsExponants[5]) {
+float units_general_h_factor(const struct UnitSystem* us,
+ float baseUnitsExponants[5]) {
float factor_exp = 0.f;
factor_exp += -baseUnitsExponants[UNIT_MASS];
@@ -362,7 +366,8 @@ float generalhFactor(const struct UnitSystem* us, float baseUnitsExponants[5]) {
* @param baseUnitsExponants The exponent of each base units required to form
* the desired quantity. See conversionFactor() for a working example
*/
-float generalaFactor(const struct UnitSystem* us, float baseUnitsExponants[5]) {
+float units_general_a_factor(const struct UnitSystem* us,
+ float baseUnitsExponants[5]) {
float factor_exp = 0.f;
factor_exp += baseUnitsExponants[UNIT_LENGTH];
@@ -379,11 +384,11 @@ float generalaFactor(const struct UnitSystem* us, float baseUnitsExponants[5]) {
* @param baseUnitsExponants The exponent of each base units required to form
* the desired quantity. See conversionFactor() for a working example
*/
-void generalConversionString(char* buffer, const struct UnitSystem* us,
- float baseUnitsExponants[5]) {
+void units_general_conversion_string(char* buffer, const struct UnitSystem* us,
+ float baseUnitsExponants[5]) {
char temp[14];
- double a_exp = generalaFactor(us, baseUnitsExponants);
- double h_exp = generalhFactor(us, baseUnitsExponants);
+ double a_exp = units_general_a_factor(us, baseUnitsExponants);
+ double h_exp = units_general_h_factor(us, baseUnitsExponants);
int i;
/* Check whether we are unitless or not */
@@ -423,12 +428,13 @@ void generalConversionString(char* buffer, const struct UnitSystem* us,
if (baseUnitsExponants[i] == 0.)
sprintf(temp, " ");
else if (baseUnitsExponants[i] == 1.)
- sprintf(temp, "%s ", getBaseUnitSymbol(i));
+ sprintf(temp, "%s ", units_get_base_unit_symbol(i));
else if (remainder(baseUnitsExponants[i], 1.) == 0)
- sprintf(temp, "%s^%d ", getBaseUnitSymbol(i),
+ sprintf(temp, "%s^%d ", units_get_base_unit_symbol(i),
(int)baseUnitsExponants[i]);
else
- sprintf(temp, "%s^%7.4f ", getBaseUnitSymbol(i), baseUnitsExponants[i]);
+ sprintf(temp, "%s^%7.4f ", units_get_base_unit_symbol(i),
+ baseUnitsExponants[i]);
strncat(buffer, temp, 12);
}
@@ -440,12 +446,12 @@ void generalConversionString(char* buffer, const struct UnitSystem* us,
if (baseUnitsExponants[i] == 0.)
continue;
else if (baseUnitsExponants[i] == 1.)
- sprintf(temp, "%s ", getBaseUnitCGSSymbol(i));
+ sprintf(temp, "%s ", units_get_base_unit_CGS_symbol(i));
else if (remainder(baseUnitsExponants[i], 1.) == 0)
- sprintf(temp, "%s^%d ", getBaseUnitCGSSymbol(i),
+ sprintf(temp, "%s^%d ", units_get_base_unit_CGS_symbol(i),
(int)baseUnitsExponants[i]);
else
- sprintf(temp, "%s^%7.4f ", getBaseUnitCGSSymbol(i),
+ sprintf(temp, "%s^%7.4f ", units_get_base_unit_CGS_symbol(i),
baseUnitsExponants[i]);
strncat(buffer, temp, 12);
}
diff --git a/src/units.h b/src/units.h
index f96824190ace6c50cbbe218e658d9b740ae3bb32..3e349dc16787cd4052a3e9205b21dce3c3732448 100644
--- a/src/units.h
+++ b/src/units.h
@@ -92,21 +92,25 @@ enum UnitConversionFactor {
UNIT_CONV_TEMPERATURE
};
-void initUnitSystem(struct UnitSystem*, const struct swift_params*);
-double getBaseUnit(const struct UnitSystem*, enum BaseUnits);
-const char* getBaseUnitSymbol(enum BaseUnits);
-const char* getBaseUnitCGSSymbol(enum BaseUnits);
-double generalConversionFactor(const struct UnitSystem* us,
- float baseUnitsExponants[5]);
-double conversionFactor(const struct UnitSystem* us,
- enum UnitConversionFactor unit);
-float generalhFactor(const struct UnitSystem* us, float baseUnitsExponants[5]);
-float hFactor(const struct UnitSystem* us, enum UnitConversionFactor unit);
-float generalaFactor(const struct UnitSystem* us, float baseUnitsExponants[5]);
-float aFactor(const struct UnitSystem* us, enum UnitConversionFactor unit);
-void generalConversionString(char* buffer, const struct UnitSystem* us,
+void units_init(struct UnitSystem*, const struct swift_params*);
+double units_get_base_unit(const struct UnitSystem*, enum BaseUnits);
+const char* units_get_base_unit_symbol(enum BaseUnits);
+const char* units_get_base_unit_CGS_symbol(enum BaseUnits);
+double units_general_conversion_factor(const struct UnitSystem* us,
+ float baseUnitsExponants[5]);
+double units_conversion_factor(const struct UnitSystem* us,
+ enum UnitConversionFactor unit);
+float units_general_h_factor(const struct UnitSystem* us,
float baseUnitsExponants[5]);
-void conversionString(char* buffer, const struct UnitSystem* us,
- enum UnitConversionFactor unit);
+float units_h_factor(const struct UnitSystem* us,
+ enum UnitConversionFactor unit);
+float units_general_a_factor(const struct UnitSystem* us,
+ float baseUnitsExponants[5]);
+float units_a_factor(const struct UnitSystem* us,
+ enum UnitConversionFactor unit);
+void units_general_conversion_string(char* buffer, const struct UnitSystem* us,
+ float baseUnitsExponants[5]);
+void units_conversion_string(char* buffer, const struct UnitSystem* us,
+ enum UnitConversionFactor unit);
#endif /* SWIFT_UNITS_H */