mirror of
				https://github.com/paboyle/Grid.git
				synced 2025-11-04 05:54:32 +00:00 
			
		
		
		
	Zero changes, literally
This commit is contained in:
		@@ -75,7 +75,7 @@ template<class vtype,class ltype,class rtype, int N> accelerator_inline void sub
 | 
			
		||||
	sub(&ret->_internal[c1][c2],&lhs->_internal,&rhs->_internal[c1][c2]);
 | 
			
		||||
      } else {
 | 
			
		||||
	// Fails -- need unary minus. Catalogue other unops?
 | 
			
		||||
	ret->_internal[c1][c2]=zero;
 | 
			
		||||
	ret->_internal[c1][c2]=Zero();
 | 
			
		||||
	ret->_internal[c1][c2]=ret->_internal[c1][c2]-rhs->_internal[c1][c2];
 | 
			
		||||
 | 
			
		||||
      }
 | 
			
		||||
 
 | 
			
		||||
@@ -272,7 +272,7 @@ public:
 | 
			
		||||
  static auto traceIndex(const iVector<vtype,N> arg) ->  iScalar<RemoveCRV(arg._internal[0])>
 | 
			
		||||
  {
 | 
			
		||||
    iScalar<RemoveCRV(arg._internal[0])> ret;
 | 
			
		||||
    ret._internal=zero;
 | 
			
		||||
    ret._internal=Zero();
 | 
			
		||||
    for(int i=0;i<N;i++){
 | 
			
		||||
      ret._internal = ret._internal+ arg._internal[i];
 | 
			
		||||
    }
 | 
			
		||||
@@ -302,7 +302,7 @@ public:
 | 
			
		||||
  static auto transposeIndex(const iMatrix<vtype,N> arg)  ->  iMatrix<vtype,N> 
 | 
			
		||||
  {
 | 
			
		||||
    iMatrix<vtype,N> ret;
 | 
			
		||||
    ret=zero;
 | 
			
		||||
    ret=Zero();
 | 
			
		||||
    for(int i=0;i<N;i++){
 | 
			
		||||
      for(int j=0;j<N;j++){
 | 
			
		||||
	ret._internal[i][j] = ret._internal[i][j]+arg._internal[i][j];
 | 
			
		||||
 
 | 
			
		||||
		Reference in New Issue
	
	Block a user