/[escript]/branches/domexper/dudley/src/Util.c
ViewVC logotype

Diff of /branches/domexper/dudley/src/Util.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3195 by jfenwick, Fri Aug 27 05:26:25 2010 UTC revision 3196 by jfenwick, Wed Sep 22 01:18:52 2010 UTC
# Line 86  void Dudley_Util_AddScatter(dim_t len,in Line 86  void Dudley_Util_AddScatter(dim_t len,in
86    
87  /*          A(1:A1,1:A2)=B(1:A1,1:B2)*C(1:B2,1:A2) */  /*          A(1:A1,1:A2)=B(1:A1,1:B2)*C(1:B2,1:A2) */
88    
89  void Dudley_Util_SmallMatMult(dim_t A1,dim_t A2, double* A, dim_t B2, double*B, double* C) {  void Dudley_Util_SmallMatMult(dim_t A1,dim_t A2, double* A, dim_t B2, const double*B, const double* C) {
90      dim_t i,j,s;      dim_t i,j,s;
91      register double rtmp;      register double rtmp;
92         for (i=0;i<A1;i++) {         for (i=0;i<A1;i++) {
# Line 102  void Dudley_Util_SmallMatMult(dim_t A1,d Line 102  void Dudley_Util_SmallMatMult(dim_t A1,d
102    
103  /*        A(1:A1,1:A2,i)=B(1:A1,1:B2,i)*C(1:B2,1:A2,i) i=1,len */  /*        A(1:A1,1:A2,i)=B(1:A1,1:B2,i)*C(1:B2,1:A2,i) i=1,len */
104    
105  void Dudley_Util_SmallMatSetMult(dim_t len,dim_t A1,dim_t A2, double* A, dim_t B2, double*B, double* C) {  void Dudley_Util_SmallMatSetMult(dim_t len,dim_t A1,dim_t A2, double* A, dim_t B2, const double*B, const double* C) {
106      dim_t q,i,j,s;      dim_t q,i,j,s;
107      register double rtmp;      register double rtmp;
108      for (q=0;q<len;q++) {      for (q=0;q<len;q++) {
# Line 119  void Dudley_Util_SmallMatSetMult(dim_t l Line 119  void Dudley_Util_SmallMatSetMult(dim_t l
119    
120  /*        A(1:A1,1:A2,i)=B(1:A1,1:B2,i)*C(1:B2,1:A2) i=1,len */  /*        A(1:A1,1:A2,i)=B(1:A1,1:B2,i)*C(1:B2,1:A2) i=1,len */
121    
122  void Dudley_Util_SmallMatSetMult1(dim_t len,dim_t A1,dim_t A2, double* A, dim_t B2, double*B, double* C) {  void Dudley_Util_SmallMatSetMult1(dim_t len,dim_t A1,dim_t A2, double* A, dim_t B2, const double*B, const double* C) {
123      dim_t q,i,j,s;      dim_t q,i,j,s;
124      register double rtmp;      register double rtmp;
125      for (q=0;q<len;q++) {      for (q=0;q<len;q++) {

Legend:
Removed from v.3195  
changed lines
  Added in v.3196

  ViewVC Help
Powered by ViewVC 1.1.26