Changeset dac8deb in tol for tol_pkg/FunRn2R/BaseC2.Cache.tol

Timestamp:

Dec 15, 2009, 5:39:55 PM (17 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

8fd1737

Parents:

8642d7e

Message:

Enabling cache for integration methods

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@1813 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol_pkg/FunRn2R/BaseC2.Cache.tol (modified) (4 diffs)

tol_pkg/FunRn2R/BaseC2.Cache.tol

@@ -17 +17 @@
   Set _.partial.derivative.1 = Copy(Empty);
   Set _.partial.derivative.2 = Copy(Empty);
+  Set _.partial.integrate.1 = Copy(Empty);
 
   //Returns the cached value of the k-th basic function for an 
@@ -37 +38 @@
     Real j_ = Min(i,j);
     SubCol(_.partial.derivative.2[i_][j_],[[k]]) 
+  };
+  //Returns the cached value of the first partial integration of k-th basic 
+  //function for an  array of points given as the rows of a matrix 
+  VMatrix partial.integrate.1.v(Real k, VMatrix X, Real i)
+  { 
+    SubCol(_.partial.integrate.1[i],[[k]]) 
   };
 
@@ -62 +69 @@
   };
 
+  //Evaluates the first partial integration of all basic functions, one in 
+  //each returned column, for an  array of points given as the rows of a 
+  //matrix 
+  VMatrix partial.integrate.1.m(VMatrix X, Real i)
+  { 
+    _.partial.integrate.1[i]
+  };
   //Some bases are more efficient using cache
-  Real _build.cache(VMatrix X, Real order)
+  Real _build.cache(VMatrix X, Real deriv.order, Real integ.order)
   {
-    VMatrix If(order>=0, { _.function.r := 
-      _build.function.m(X)
-    });
-    Set If(order>=1, { _.partial.derivative.1 := 
+    VMatrix _.function := _build.function.m(X);
+    Set If(deriv.order>=1, { _.partial.derivative.1 := 
       For(1, _.n, VMatrix(Real i)
       { _build.partial.derivative.1.m(X,i) })
     });
-    Set If(order>=2, { _.partial.derivative.2 := 
+    Set If(deriv.order>=2, { _.partial.derivative.2 := 
       For(1, _.n, Set(Real i)
       { 
@@ -81 +93 @@
       })
     });
+    Set If(integ.order>=1, { _.partial.integrate.1 := 
+      For(1, _.n, VMatrix(Real i)
+      { _build.partial.integrate.1.m(X,i) })
+    });
     True
   }