/[escript]/branches/doubleplusgood/paso/src/PasoUtil.h
ViewVC logotype

Diff of /branches/doubleplusgood/paso/src/PasoUtil.h

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

revision 1388 by trankine, Fri Jan 11 07:45:58 2008 UTC revision 4154 by jfenwick, Tue Jan 22 09:30:23 2013 UTC
# Line 1  Line 1 
1    
2  /* $Id$ */  /*****************************************************************************
3    *
4    * Copyright (c) 2003-2013 by University of Queensland
5    * http://www.uq.edu.au
6    *
7    * Primary Business: Queensland, Australia
8    * Licensed under the Open Software License version 3.0
9    * http://www.opensource.org/licenses/osl-3.0.php
10    *
11    * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12    * Development since 2012 by School of Earth Sciences
13    *
14    *****************************************************************************/
15    
 /*******************************************************  
  *  
  *           Copyright 2003-2007 by ACceSS MNRF  
  *       Copyright 2007 by University of Queensland  
  *  
  *                http://esscc.uq.edu.au  
  *        Primary Business: Queensland, Australia  
  *  Licensed under the Open Software License version 3.0  
  *     http://www.opensource.org/licenses/osl-3.0.php  
  *  
  *******************************************************/  
16    
17  #ifndef INC_PASO_UTIL  #ifndef INC_PASO_UTIL
18  #define INC_PASO_UTIL  #define INC_PASO_UTIL
19    
20  /**************************************************************/  /************************************************************************************/
21    
22  /*   Some utility routines: */  /*   Some utility routines: */
23    
24  /**************************************************************/  /************************************************************************************/
25    
26  /*   Copyrights by ACcESS Australia, 2003,2004,2005 */  /*   Copyrights by ACcESS Australia, 2003,2004,2005 */
27  /*   author: gross@access.edu.au */  /*   author: l.gross@uq.edu.au */
28    
29  /**************************************************************/  /************************************************************************************/
30    
31  #include "Common.h"  #include "Common.h"
32    #include "esysUtils/Esys_MPI.h"
33    #ifdef _OPENMP
34    #include <omp.h>
35    #endif
36    
37  /**************************************************************/  /************************************************************************************/
38    
39  index_t Paso_Util_cumsum(dim_t,index_t*);  index_t Paso_Util_cumsum(dim_t,index_t*);
 void Paso_copyDouble(dim_t n,double* source, double* target);  
40  bool_t Paso_Util_isAny(dim_t N,index_t* array,index_t value);  bool_t Paso_Util_isAny(dim_t N,index_t* array,index_t value);
41    void Paso_zeroes(const dim_t n, double* x);
42    void Paso_Update(const dim_t n, const double a, double* x, const double b, const double* y);
43    void Paso_LinearCombination(const dim_t n, double*z, const double a,const double* x, const double b, const double* y);
44    double Paso_InnerProduct(const dim_t n,const double* x, const double* y, Esys_MPIInfo* mpiinfo);
45    double Paso_l2(const dim_t n, const double* x, Esys_MPIInfo* mpiinfo);
46    void ApplyGivensRotations(const dim_t n,double* v,const double* c,const double* s);
47    void Paso_Copy(const dim_t n, double* out, const double* in);
48    bool_t Paso_fileExists( const char* filename );
49    double Paso_lsup(const dim_t n, const double* x, Esys_MPIInfo* mpiinfo);
50    index_t Paso_Util_cumsum_maskedTrue(dim_t N,index_t* array, bool_t* mask);
51    index_t Paso_Util_cumsum_maskedFalse(dim_t N,index_t* array, bool_t* mask);
52    index_t Paso_Util_arg_max(dim_t n, dim_t* lambda);
53    index_t Paso_Util_iMax(const dim_t N,const index_t* array);
54    dim_t Paso_Util_numPositives(const dim_t N, const double *x);
55    
56    #define Paso_Scale(n, x, a) Paso_Update(n, a, x, 0, x);
57    #define Paso_AXPY(n, x, a, y) Paso_Update(n, 1., x, a,  y);
58    #define Paso_copyShortDouble(n, source, target)  memcpy(target,source,sizeof(double)*(size_t)n)
59    
 #endif /* #ifndef INC_PASO_UTIL */  
60    
61  /*  #endif /* #ifndef INC_PASO_UTIL */
  * $Log$  
  * Revision 1.2  2005/09/15 03:44:39  jgs  
  * Merge of development branch dev-02 back to main trunk on 2005-09-15  
  *  
  * Revision 1.1.2.1  2005/09/05 06:29:48  gross  
  * These files have been extracted from finley to define a stand alone libray for iterative  
  * linear solvers on the ALTIX. main entry through Paso_solve. this version compiles but  
  * has not been tested yet.  
  *  
  *  
  */  

Legend:
Removed from v.1388  
changed lines
  Added in v.4154

  ViewVC Help
Powered by ViewVC 1.1.26