Changeset 2291 in Sophya for trunk/SophyaLib/Samba/lambdaBuilder.h

Timestamp:

Dec 5, 2002, 4:43:59 PM (23 years ago)

Author:

lemeur

Message:

protections de dimensions

File:

• trunk/SophyaLib/Samba/lambdaBuilder.h (modified) (5 diffs)

trunk/SophyaLib/Samba/lambdaBuilder.h

@@ -14 +14 @@
 
  public :
-Legendre();
 Legendre(r_8 x, int_4 lmax);
- inline r_8 getPl(int_4 l) 
+
+inline r_8 getPl(int_4 l) const
 {
   if (l>lmax_) 
@@ -25 +25 @@
 }
 
+
  private : 
-  void array_init(int_4 lmax);
+
+Legendre() {;}
+
+void array_init(int_4 lmax);
 
   r_8 x_;
@@ -45 +49 @@
 
  /*! return the value of the coefficient \f$  \lambda_l^m \f$ */
-inline double lamlm(int l, int m) const { return lambda_(l,m); }
+inline double lamlm(int l, int m) const 
+  {
+    lambda_.CheckRelativeIndices(l,m);
+    return lambda_(l,m); 
+  }
 
-inline double lamlm(int index) const { return lambda_(index); }
+inline double lamlm(int index) const {
+ return lambda_(index); 
+}
 
 
@@ -92 +102 @@
 inline void lam_wx(int l, int m, r_8& w, r_8& x) const
    {
+  if (m > l ) 
+    {
+      throw RangeCheckError("LambdaWXBuilder::lam_wx : l < m !" );
+    }
      w=lamWlm_(l,m);
      x=lamXlm_(l,m);
@@ -117 +131 @@
 inline void lam_pm(int l, int m, r_8& lambda_plus, r_8& lambda_moins) const
    {
+  if (m > l ) 
+    {
+      throw RangeCheckError("LambdaPMBuilder::lam_pm : l < m !" );
+    }
       lambda_plus = lamPlm_(l,m);
      lambda_moins = lamMlm_(l,m);