Changeset 2976 in Sophya

Timestamp:

Jun 20, 2006, 6:05:40 PM (19 years ago)

Author:

ansari

Message:

ajout commande resol2szidx ds skymapmodule.cc , Reza 20/6/2006

File:

• trunk/SophyaPI/ProgPI/skymapmodule.cc (modified) (4 diffs)

trunk/SophyaPI/ProgPI/skymapmodule.cc

@@ -46 +46 @@
   virtual int   Execute(string& keyw, vector<string>& args, string& toks);
 
+  void Resol2SizeIndex(vector<string>& tokens);
   void TypeMap(vector<string>& tokens);
   void Map2Double(vector<string>& tokens);
@@ -99 +100 @@
 mpiac->RegisterCommand(kw, usage, this, hgrp);
 
+kw = "resol2szidx";
+usage = "Compute SizeIndex value (=nside for HEALPix) for a";
+usage += "\n given resolution, (resol in arcminutes)";
+usage += "\n Usage: resol2szidx resol";
+mpiac->RegisterCommand(kw, usage, this, hgrp);
+
 kw = "typemap";
 usage = "Imprime le type de map";
@@ -274 +281 @@
 if(kw == "settypemap") {
   SetTypeMap(tokens);
+} else if(kw == "resol2szidx") {
+  if(tokens.size()<1) {
+    cout<<"Usage: resol2szidx resol"<<endl;
+    return(0);
+  }
+  Resol2SizeIndex(tokens);
 } else if(kw == "typemap") {
   if(tokens.size()<1) {
@@ -415 +428 @@
 
 return(0);
+}
+
+/* --Methode-- */
+void skymapmoduleExecutor::Resol2SizeIndex(vector<string>& tokens)
+{
+ double resamin = atof(tokens[0].c_str());
+ double resrad = Angle(resamin,Angle::ArcMin).ToRadian();
+ int_4 mh = SphereHEALPix<r_4>::ResolToSizeIndex(resrad);
+ int_4 mt = SphereThetaPhi<r_4>::ResolToSizeIndex(resrad);
+ cout << "Resol2SizeIndex: Resol= " << resamin << " arcmin ->" 
+      << resrad << " radian  NPix ~= " << 4.*M_PI/(resrad*resrad) 
+      << " \n  ===> nside=m_HEALPix= " << mh << "  m_ThetaPhi= " << mt << endl;
+ return;
 }