Changeset 668 in Sophya for trunk/SophyaLib/SkyT/gaussfilt.cc

Timestamp:

Nov 29, 1999, 3:16:10 PM (26 years ago)

Author:

ansari

Message:

Ajout de classes deleguees PPersist et correction integration - Sophie 29/11/99

File:

• trunk/SophyaLib/SkyT/gaussfilt.cc (modified) (4 diffs)

trunk/SophyaLib/SkyT/gaussfilt.cc

@@ -1 +1 @@
 //--------------------------------------------------------------------------
 // File and Version Information:
-//      $Id: gaussfilt.cc,v 1.3 1999-11-21 23:25:45 ansari Exp $
+//      $Id: gaussfilt.cc,v 1.4 1999-11-29 14:16:04 ansari Exp $
 //
 // Description:
@@ -22 +22 @@
 // Constructor --
 //----------------
+GaussianFilter::GaussianFilter()
+  : SpectralResponse()
+{
+  setParams(100, 10., 1.);
+}
+
 GaussianFilter::GaussianFilter(double nu0, double s, double a, double numin, double numax)
   : SpectralResponse(numin, numax)
 {
-  if (s < 1.e-19) s = 1.e-19;
-  _s = s;
-  _nu0 = nu0;
-  _a = a;
+  setParams(nu0, s, a);
 }
 
@@ -66 +69 @@
 }
 
+
+
+void 
+GaussianFilter::setParams(double nu0, double s, double a)
+{
+  if (s < 1.e-19) s = 1.e-19;
+  _s = s;
+  _nu0 = nu0;
+  _a = a;
+}
+
+
+
 void 
 GaussianFilter::Print(ostream& os)  const
@@ -76 +92 @@
 
 }
+
+void
+ObjFileIO<GaussianFilter>::WriteSelf(POutPersist& s) const
+{
+  if(dobj == NULL)
+    {
+      cout << " ObjFileIO<GaussianFilter>::WriteSelf:: dobj= null " << endl;
+      return;
+    }
+  int_4 version, nothing;
+  version = 1;
+  nothing = 0;   // Reserved for future use 
+  s.PutI4(version);
+  s.PutI4(nothing);
+
+  s.PutR8(dobj->minFreq());
+  s.PutR8(dobj->maxFreq());
+  s.PutR8(dobj->giveNorm());
+  s.PutR8(dobj->giveNu0());
+  s.PutR8(dobj->giveDNu());
+}
+
+void
+ObjFileIO<GaussianFilter>::ReadSelf(PInPersist& s)  
+{
+  int_4 version, nothing;
+  version = 1;
+  nothing = 0;   // Reserved for future use 
+  s.GetI4(version);
+  s.GetI4(nothing);
+
+  if(dobj == NULL)
+    {
+      dobj= new GaussianFilter();
+      ownobj= true;
+    }
+
+  r_8 a, nu0, dnu, numin, numax;
+  s.GetR8(numin);
+  s.GetR8(numax);
+  s.GetR8(a);
+  s.GetR8(nu0);
+  s.GetR8(dnu);
+  dobj->setMinMaxFreq(numin, numax);
+  dobj->setParams(nu0, dnu, a);
+  //  cout << " Norm - Nu0 - DNu - minFreq - maxFreq " << endl; 
+  // cout << _a << "-" << _nu0  << "-" << _s << "-" << _numin << "-" << _numax << endl;
+}
+
+
+#ifdef __CXX_PRAGMA_TEMPLATES__
+#pragma define_template ObjFileIO<GaussianFilter>
+#endif
+
+#if defined(ANSI_TEMPLATES) || defined(GNU_TEMPLATES)
+template class ObjFileIO<GaussianFilter>;
+#endif