Add dOptionsRealUnits()

jedbrown · May 25, 2011 · 0f22606 · 0f22606
1 parent b77a3f3
commit 0f22606
Show file tree

Hide file tree

Showing 2 changed files with 21 additions and 1 deletion.
diff --git a/include/dohpunits.h b/include/dohpunits.h
@@ -27,6 +27,6 @@ dScalar dUnitDimensionalize(dUnit,dScalar);
 dScalar dUnitNonDimensionalize(dUnit,dScalar);
 dScalar dUnitDimensionalizeSI(dUnit,dScalar);
 dScalar dUnitNonDimensionalizeSI(dUnit,dScalar);
-
+dErr dOptionsRealUnits(const char opt[],const char text[],const char man[],dUnit unit,PetscReal defaultv,PetscReal *value,PetscBool *set);
 dEXTERN_C_END
 #endif
diff --git a/src/util/dunits.c b/src/util/dunits.c
@@ -277,3 +277,23 @@ dErr dUnitsDestroy(dUnits *unp)
   err = PetscHeaderDestroy(unp);dCHK(err);
   dFunctionReturn(0);
 }
+dErr dOptionsRealUnits(const char opt[],const char text[],const char man[],dUnit unit,PetscReal defaultv,PetscReal *value,PetscBool *set)
+{
+  char text2[512];
+  PetscReal default2,value2;
+  dBool set2;
+  dErr err;
+
+  dFunctionBegin;
+  if (unit) {
+    default2 = dUnitDimensionalize(unit,defaultv);
+    err = PetscSNPrintf(text2,sizeof text2,"%s (%s)",text,dUnitName(unit));dCHK(err);
+    err = PetscOptionsReal(opt,text2,man,default2,&value2,&set2);dCHK(err);
+    if (set2) *value = dUnitNonDimensionalize(unit,value2);
+    if (set) *set = set2;
+  } else {
+    err = PetscSNPrintf(text2,sizeof text2,"%s (nondimensional)",text);dCHK(err);
+    err = PetscOptionsReal(opt,text2,man,defaultv,value,set);dCHK(err);
+  }
+  dFunctionReturn(0);
+}