/[escript]/trunk/escript/src/Data.h
ViewVC logotype

Diff of /trunk/escript/src/Data.h

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

trunk/esys2/escript/src/Data/Data.h revision 97 by jgs, Tue Dec 14 05:39:33 2004 UTC trunk/escript/src/Data/Data.h revision 286 by gross, Thu Dec 1 06:40:20 2005 UTC
# Line 1  Line 1 
1  // $Id$  // $Id$
2  /*=============================================================================  /*
   
3   ******************************************************************************   ******************************************************************************
4   *                                                                            *   *                                                                            *
5   *       COPYRIGHT ACcESS 2004 -  All Rights Reserved                         *   *       COPYRIGHT ACcESS 2004 -  All Rights Reserved                         *
# Line 12  Line 11 
11   * person has a software license agreement with ACcESS.                       *   * person has a software license agreement with ACcESS.                       *
12   *                                                                            *   *                                                                            *
13   ******************************************************************************   ******************************************************************************
14    */
15    
16  ******************************************************************************/  /** \file Data.h */
17    
18  #ifndef DATA_H  #ifndef DATA_H
19  #define DATA_H  #define DATA_H
20    
21  #include "escript/Data/DataAbstract.h"  #include "escript/Data/DataAbstract.h"
22  #include "escript/Data/DataTagged.h"  #include "escript/Data/DataTagged.h"
23    #include "escript/Data/DataAlgorithm.h"
24  #include "escript/Data/FunctionSpace.h"  #include "escript/Data/FunctionSpace.h"
25  #include "escript/Data/BinaryOp.h"  #include "escript/Data/BinaryOp.h"
26  #include "escript/Data/UnaryOp.h"  #include "escript/Data/UnaryOp.h"
27    #include "escript/Data/DataException.h"
28    
29  extern "C" {  extern "C" {
30  #include "escript/Data/DataC.h"  #include "escript/Data/DataC.h"
# Line 39  extern "C" { Line 41  extern "C" {
41  #include <boost/python/tuple.hpp>  #include <boost/python/tuple.hpp>
42  #include <boost/python/numeric.hpp>  #include <boost/python/numeric.hpp>
43    
44    namespace escript {
45    
46    //
47    // Forward declaration for various implementations of Data.
48    class DataEmpty;
49    class DataConstant;
50    class DataTagged;
51    class DataExpanded;
52    
53  /**  /**
54     \brief     \brief
55     Data is essentially a factory class which creates the appropriate Data     Data creates the appropriate Data object for the given construction
56     object for the given construction arguments. It retains control over     arguments.
    the object created for the lifetime of the object.  
    The type of Data object referred to may change during the lifetime of  
    the Data object.  
57    
58     Description:     Description:
59     Data is essentially a factory class which creates the appropriate Data     Data is essentially a factory class which creates the appropriate Data
# Line 54  extern "C" { Line 62  extern "C" {
62     The type of Data object referred to may change during the lifetime of     The type of Data object referred to may change during the lifetime of
63     the Data object.     the Data object.
64  */  */
   
 namespace escript {  
   
 //  
 // Forward declaration for various implimentations of Data.  
 class DataEmpty;  
 class DataConstant;  
 class DataTagged;  
 class DataExpanded;  
   
65  class Data {  class Data {
66    
67    public:    public:
68    
69      // These typedefs allow function names to be cast to pointers
70      // to functions of the appropriate type when calling unaryOp etc.
71    typedef double (*UnaryDFunPtr)(double);    typedef double (*UnaryDFunPtr)(double);
72    typedef double (*BinaryDFunPtr)(double,double);    typedef double (*BinaryDFunPtr)(double,double);
73    
# Line 206  class Data { Line 206  class Data {
206         const boost::python::tuple& shape=boost::python::make_tuple(),         const boost::python::tuple& shape=boost::python::make_tuple(),
207         const FunctionSpace& what=FunctionSpace(),         const FunctionSpace& what=FunctionSpace(),
208         bool expanded=false);         bool expanded=false);
209      /**
210         \brief
211         Destructor
212      */
213      ~Data();
214    
215    /**    /**
216       \brief       \brief
# Line 220  class Data { Line 225  class Data {
225    
226    /**    /**
227       \brief       \brief
228         Return the values of all data-points as a single python numarray object.
229      */
230      const boost::python::numeric::array
231      convertToNumArray();
232    
233      /**
234         \brief
235         Return the values of all data-points for the given sample as a single python numarray object.
236      */
237      const boost::python::numeric::array
238      convertToNumArrayFromSampleNo(int sampleNo);
239    
240      /**
241         \brief
242         Return the value of the specified data-point as a single python numarray object.
243      */
244      const boost::python::numeric::array
245      convertToNumArrayFromDPNo(int sampleNo,
246                                int dataPointNo);
247    
248      /**
249         \brief
250         Fills the expanded Data object from values of a python numarray object.
251      */
252      void
253      fillFromNumArray(const boost::python::numeric::array);
254    
255      /**
256         \brief
257         Return the tag number associated with the given data-point.
258    
259         The data-point number here corresponds to the data-point number in the
260         numarray returned by convertToNumArray.
261      */
262      int
263      getTagNumber(int dpno);
264    
265      /**
266         \brief
267       Return the C wrapper for the Data object.       Return the C wrapper for the Data object.
268    */    */
269    escriptDataC    escriptDataC
# Line 397  class Data { Line 441  class Data {
441    
442    /**    /**
443       \brief       \brief
444         Assign the given value to the data-points referenced by the given
445         reference number.
446    
447         The value supplied is a python numarray object.  The data from this numarray
448         is unpacked into a DataArray, and this is used to set the corresponding
449         data-points in the underlying Data object.
450    
451         If the underlying Data object cannot be accessed via reference numbers, an
452         exception will be thrown.
453    
454         \param ref - Input - reference number.
455         \param value - Input - value to assign to data-points associated with
456                                the given reference number.
457      */
458      void
459      setRefValue(int ref,
460                  const boost::python::numeric::array& value);
461    
462      /**
463         \brief
464         Return the values associated with the data-points referenced by the given
465         reference number.
466    
467         The value supplied is a python numarray object. The data from the corresponding
468         data-points in this Data object are packed into the given numarray object.
469    
470         If the underlying Data object cannot be accessed via reference numbers, an
471         exception will be thrown.
472    
473         \param ref - Input - reference number.
474         \param value - Output - object to receive values from data-points
475                                 associated with the given reference number.
476      */
477      void
478      getRefValue(int ref,
479                  boost::python::numeric::array& value);
480    
481      /**
482         \brief
483       Return a view into the data for the data point specified.       Return a view into the data for the data point specified.
484       NOTE: Construction of the DataArrayView is a relatively expensive       NOTE: Construction of the DataArrayView is a relatively expensive
485       operation.       operation.
# Line 422  class Data { Line 505  class Data {
505       \brief       \brief
506       Return the data point shape as a tuple of integers.       Return the data point shape as a tuple of integers.
507    */    */
508    boost::python::tuple    const boost::python::tuple
509    getShapeTuple() const;    getShapeTuple() const;
510    
511    /**    /**
# Line 459  class Data { Line 542  class Data {
542       cannot be converted to a DataTagged object.       cannot be converted to a DataTagged object.
543       \param tagKey - Input - Integer key.       \param tagKey - Input - Integer key.
544       \param value - Input - Value to associate with given key.       \param value - Input - Value to associate with given key.
      Note: removed for now - this version not needed, and breaks escript.cpp  
545    */    */
   /*  
546    void    void
547    setTaggedValue(int tagKey,    setTaggedValueFromCPP(int tagKey,
548                   const DataArrayView& value);                          const DataArrayView& value);
   */  
549    
550    /**    /**
551      \brief      \brief
# Line 483  class Data { Line 563  class Data {
563       \brief       \brief
564       Interpolates this onto the given functionspace and returns       Interpolates this onto the given functionspace and returns
565       the result as a Data object.       the result as a Data object.
566         *
567    */    */
568    Data    Data
569    interpolate(const FunctionSpace& functionspace) const;    interpolate(const FunctionSpace& functionspace) const;
# Line 491  class Data { Line 572  class Data {
572       \brief       \brief
573       Calculates the gradient of the data at the data points of functionspace.       Calculates the gradient of the data at the data points of functionspace.
574       If functionspace is not present the function space of Function(getDomain()) is used.       If functionspace is not present the function space of Function(getDomain()) is used.
575         *
576    */    */
577    Data    Data
578    gradOn(const FunctionSpace& functionspace) const;    gradOn(const FunctionSpace& functionspace) const;
# Line 501  class Data { Line 583  class Data {
583    /**    /**
584       \brief       \brief
585       Calculate the integral over the function space domain.       Calculate the integral over the function space domain.
586         *
587    */    */
588    boost::python::numeric::array    boost::python::numeric::array
589    integrate() const;    integrate() const;
# Line 508  class Data { Line 591  class Data {
591    /**    /**
592       \brief       \brief
593       Return a Data with a 1 for +ive values and a 0 for 0 or -ive values.       Return a Data with a 1 for +ive values and a 0 for 0 or -ive values.
594         *
595    */    */
596    Data    Data
597    wherePositive() const;    wherePositive() const;
# Line 515  class Data { Line 599  class Data {
599    /**    /**
600       \brief       \brief
601       Return a Data with a 1 for -ive values and a 0 for +ive or 0 values.       Return a Data with a 1 for -ive values and a 0 for +ive or 0 values.
602         *
603    */    */
604    Data    Data
605    whereNegative() const;    whereNegative() const;
# Line 522  class Data { Line 607  class Data {
607    /**    /**
608       \brief       \brief
609       Return a Data with a 1 for +ive or 0 values and a 0 for -ive values.       Return a Data with a 1 for +ive or 0 values and a 0 for -ive values.
610         *
611    */    */
612    Data    Data
613    whereNonNegative() const;    whereNonNegative() const;
# Line 529  class Data { Line 615  class Data {
615    /**    /**
616       \brief       \brief
617       Return a Data with a 1 for -ive or 0 values and a 0 for +ive values.       Return a Data with a 1 for -ive or 0 values and a 0 for +ive values.
618         *
619    */    */
620    Data    Data
621    whereNonPositive() const;    whereNonPositive() const;
# Line 536  class Data { Line 623  class Data {
623    /**    /**
624       \brief       \brief
625       Return a Data with a 1 for 0 values and a 0 for +ive or -ive values.       Return a Data with a 1 for 0 values and a 0 for +ive or -ive values.
626         *
627    */    */
628    Data    Data
629    whereZero() const;    whereZero() const;
# Line 543  class Data { Line 631  class Data {
631    /**    /**
632       \brief       \brief
633       Return a Data with a 0 for 0 values and a 1 for +ive or -ive values.       Return a Data with a 0 for 0 values and a 1 for +ive or -ive values.
634         *
635    */    */
636    Data    Data
637    whereNonZero() const;    whereNonZero() const;
638    
639    /**    /**
640       \brief       \brief
641         Return the maximum absolute value of this Data object.
642         *
643      */
644      double
645      Lsup() const;
646    
647      /**
648         \brief
649         Return the minimum absolute value of this Data object.
650         *
651      */
652      double
653      Linf() const;
654    
655      /**
656         \brief
657         Return the maximum value of this Data object.
658         *
659      */
660      double
661      sup() const;
662    
663      /**
664         \brief
665         Return the minimum value of this Data object.
666         *
667      */
668      double
669      inf() const;
670    
671      /**
672         \brief
673         Return the absolute value of each data point of this Data object.
674         *
675      */
676      Data
677      abs() const;
678    
679      /**
680         \brief
681         Return the maximum value of each data point of this Data object.
682         *
683      */
684      Data
685      maxval() const;
686    
687      /**
688         \brief
689         Return the minimum value of each data point of this Data object.
690         *
691      */
692      Data
693      minval() const;
694    
695      /**
696         \brief
697         Return the (sample number, data-point number) of the data point with
698         the minimum value in this Data object.
699      */
700      const boost::python::tuple
701      mindp() const;
702    
703      void
704      calc_mindp(int& SampleNo,
705                 int& DataPointNo) const;
706    
707      /**
708         \brief
709         Return the sign of each data point of this Data object.
710         -1 for negative values, zero for zero values, 1 for positive values.
711         *
712      */
713      Data
714      sign() const;
715    
716      /**
717         \brief
718         Transpose each data point of this Data object around the given axis.
719         --* not implemented yet *--
720         *
721      */
722      Data
723      transpose(int axis) const;
724    
725      /**
726         \brief
727         Calculate the trace of each data point of this Data object.
728         sum(A[i,i,i,i])
729         *
730      */
731      Data
732      trace() const;
733    
734      /**
735         \brief
736       Return the sin of each data point of this Data object.       Return the sin of each data point of this Data object.
737         *
738    */    */
739    Data    Data
740    sin() const;    sin() const;
# Line 557  class Data { Line 742  class Data {
742    /**    /**
743       \brief       \brief
744       Return the cos of each data point of this Data object.       Return the cos of each data point of this Data object.
745         *
746    */    */
747    Data    Data
748    cos() const;    cos() const;
# Line 564  class Data { Line 750  class Data {
750    /**    /**
751       \brief       \brief
752       Return the tan of each data point of this Data object.       Return the tan of each data point of this Data object.
753         *
754    */    */
755    Data    Data
756    tan() const;    tan() const;
757    
758    /**    /**
759       \brief       \brief
760       Return the log to base 10 of each data point of this Data object.       Return the asin of each data point of this Data object.
761         *
762    */    */
763    Data    Data
764    log() const;    asin() const;
765    
766    /**    /**
767       \brief       \brief
768       Return the natural log of each data point of this Data object.       Return the acos of each data point of this Data object.
769         *
770    */    */
771    Data    Data
772    ln() const;    acos() const;
773    
774    /**    /**
775       \brief       \brief
776       Return the maximum absolute value of this Data object.       Return the atan of each data point of this Data object.
777         *
778    */    */
779    double    Data
780    Lsup() const;    atan() const;
781    
782    /**    /**
783       \brief       \brief
784       Return the maximum value of this Data object.       Return the sinh of each data point of this Data object.
785         *
786    */    */
787    double    Data
788    sup() const;    sinh() const;
789    
790    /**    /**
791       \brief       \brief
792       Return the minimum value of this Data object.       Return the cosh of each data point of this Data object.
793         *
794    */    */
795    double    Data
796    inf() const;    cosh() const;
797    
798    /**    /**
799       \brief       \brief
800       Return the absolute value of each data point of this Data object.       Return the tanh of each data point of this Data object.
801         *
802    */    */
803    Data    Data
804    abs() const;    tanh() const;
805    
806    /**    /**
807       \brief       \brief
808       Return the maximum value of each data point of this Data object.       Return the asinh of each data point of this Data object.
809         *
810    */    */
811    Data    Data
812    maxval() const;    asinh() const;
813    
814    /**    /**
815       \brief       \brief
816       Return the minimum value of each data point of this Data object.       Return the acosh of each data point of this Data object.
817         *
818    */    */
819    Data    Data
820    minval() const;    acosh() const;
821    
822    /**    /**
823       \brief       \brief
824       Return the length of each data point of this Data object.       Return the atanh of each data point of this Data object.
825       sqrt(sum(A[i,j,k,l]^2))       *
826    */    */
827    Data    Data
828    length() const;    atanh() const;
829    
830    /**    /**
831       \brief       \brief
832       Return the sign of each data point of this Data object.       Return the log to base 10 of each data point of this Data object.
833       -1 for negative values, zero for zero values, 1 for positive values.       *
834    */    */
835    Data    Data
836    sign() const;    log10() const;
837    
838    /**    /**
839      \transpose       \brief
840      Transpose each data point of this Data object around the given axis.       Return the natural log of each data point of this Data object.
841         *
842    */    */
843    Data    Data
844    transpose(int axis) const;    log() const;
845    
846    /**    /**
847      \trace       \brief
848      Calculate the trace of each data point of this Data object.       Return the exponential function of each data point of this Data object.
849      sum(A[i,i,i,i])       *
850    */    */
851    Data    Data
852    trace() const;    exp() const;
853    
854    /**    /**
855      \exp       \brief
856      Return the exponential function of each data point of this Data object.       Return the square root of each data point of this Data object.
857         *
858    */    */
859    Data    Data
860    exp() const;    sqrt() const;
861    
862      /**
863         \brief
864         Return the negation of each data point of this Data object.
865         *
866      */
867      Data
868      neg() const;
869    
870    /**    /**
871      \sqrt       \brief
872      Return the square root of each data point of this Data object.       Return the identity of each data point of this Data object.
873         Simply returns this object unmodified.
874         *
875    */    */
876    Data    Data
877    sqrt() const;    pos() const;
878    
879    /**    /**
880       \brief       \brief
881       Return the given power of each data point of this Data object.       Return the given power of each data point of this Data object.
882    
883         \param right Input - the power to raise the object to.
884         *
885    */    */
886    Data    Data
887    powD(const Data& right) const;    powD(const Data& right) const;
888    
889      /**
890         \brief
891         Return the given power of each data point of this boost python object.
892        
893         \param right Input - the power to raise the object to.
894         *
895       */
896    Data    Data
897    powO(const boost::python::object& right) const;    powO(const boost::python::object& right) const;
898    
899    /**    /**
900      \brief       \brief
901      Return the negation of each data point of this Data object.       writes the object to a file in the DX file format
902    */    */
903    Data    void
904    neg() const;    saveDX(std::string fileName) const;
905    
906    /**    /**
907      \brief       \brief
908      Return the identity of each data point of this Data object.       writes the object to a file in the VTK file format
     Simply returns this object unmodified.  
909    */    */
910    Data    void
911    pos() const;    saveVTK(std::string fileName) const;
912    
913    /**    /**
914       \brief       \brief
915       Overloaded operator +=       Overloaded operator +=
916       \param right - Input - The right hand side.       \param right - Input - The right hand side.
917         *
918    */    */
919    Data& operator+=(const Data& right);    Data& operator+=(const Data& right);
920    Data& operator+=(const boost::python::object& right);    Data& operator+=(const boost::python::object& right);
# Line 706  class Data { Line 923  class Data {
923       \brief       \brief
924       Overloaded operator -=       Overloaded operator -=
925       \param right - Input - The right hand side.       \param right - Input - The right hand side.
926         *
927    */    */
928    Data& operator-=(const Data& right);    Data& operator-=(const Data& right);
929    Data& operator-=(const boost::python::object& right);    Data& operator-=(const boost::python::object& right);
# Line 714  class Data { Line 932  class Data {
932       \brief       \brief
933       Overloaded operator *=       Overloaded operator *=
934       \param right - Input - The right hand side.       \param right - Input - The right hand side.
935         *
936    */    */
937    Data& operator*=(const Data& right);    Data& operator*=(const Data& right);
938    Data& operator*=(const boost::python::object& right);    Data& operator*=(const boost::python::object& right);
# Line 722  class Data { Line 941  class Data {
941       \brief       \brief
942       Overloaded operator /=       Overloaded operator /=
943       \param right - Input - The right hand side.       \param right - Input - The right hand side.
944         *
945    */    */
946    Data& operator/=(const Data& right);    Data& operator/=(const Data& right);
947    Data& operator/=(const boost::python::object& right);    Data& operator/=(const boost::python::object& right);
# Line 755  class Data { Line 975  class Data {
975       \brief       \brief
976       Copies slice from value into this Data object.       Copies slice from value into this Data object.
977    
      \description  
978       Implements the [] set operator in python.       Implements the [] set operator in python.
979       Calls setSlice.       Calls setSlice.
980    
# Line 788  class Data { Line 1007  class Data {
1007       Return a Data object containing the specified slice of       Return a Data object containing the specified slice of
1008       this Data object.       this Data object.
1009       \param region - Input - Region to copy.       \param region - Input - Region to copy.
1010         *
1011    */    */
1012    Data    Data
1013    getSlice(const DataArrayView::RegionType& region) const;    getSlice(const DataArrayView::RegionType& region) const;
# Line 798  class Data { Line 1018  class Data {
1018       Data object.       Data object.
1019       \param value - Input - Data to copy from.       \param value - Input - Data to copy from.
1020       \param region - Input - Region to copy.       \param region - Input - Region to copy.
1021         *
1022    */    */
1023    void    void
1024    setSlice(const Data& value,    setSlice(const Data& value,
1025             const DataArrayView::RegionType& region);             const DataArrayView::RegionType& region);
1026    
1027      /**
1028         \brief
1029         Archive the current Data object to the given file.
1030         \param fileName - Input - file to archive to.
1031      */
1032      void
1033      archiveData(const std::string fileName);
1034    
1035      /**
1036         \brief
1037         Extract the Data object archived in the given file, overwriting
1038         the current Data object.
1039         Note - the current object must be of type DataEmpty.
1040         \param fileName - Input - file to extract from.
1041         \param fspace - Input - a suitable FunctionSpace descibing the data.
1042      */
1043      void
1044      extractData(const std::string fileName,
1045                  const FunctionSpace& fspace);
1046    
1047   protected:   protected:
1048    
1049   private:   private:
# Line 823  class Data { Line 1064  class Data {
1064    /**    /**
1065       \brief       \brief
1066       Perform the specified reduction algorithm on every element of every data point in       Perform the specified reduction algorithm on every element of every data point in
1067       this Data object and return the single double value result.       this Data object according to the given function and return the single value result.
1068    */    */
1069    template <class UnaryFunction>    template <class BinaryFunction>
1070    inline    inline
1071    double    double
1072    algorithm(UnaryFunction operation) const;    algorithm(BinaryFunction operation,
1073                double initial_value) const;
1074    
1075      /**
1076         \brief
1077         Reduce each data-point in this Data object using the given operation. Return a Data
1078         object with the same number of data-points, but with each data-point containing only
1079         one value - the result of the reduction operation on the corresponding data-point in
1080         this Data object
1081      */
1082      template <class BinaryFunction>
1083      inline
1084      Data
1085      dp_algorithm(BinaryFunction operation,
1086                   double initial_value) const;
1087    
1088    /**    /**
1089       \brief       \brief
# Line 895  class Data { Line 1150  class Data {
1150    // pointer to the actual data object    // pointer to the actual data object
1151    boost::shared_ptr<DataAbstract> m_data;    boost::shared_ptr<DataAbstract> m_data;
1152    
1153      //
1154      // pointer to the internal profiling data
1155      struct profDataEntry *profData;
1156    
1157  };  };
1158    
1159  template <class IValueType>  template <class IValueType>
# Line 1190  Data::unaryOp(UnaryFunction operation) Line 1449  Data::unaryOp(UnaryFunction operation)
1449    object (*this) is a rank n Data object, and returned object is a scalar.    object (*this) is a rank n Data object, and returned object is a scalar.
1450    Calls escript::algorithm.    Calls escript::algorithm.
1451  */  */
1452  template <class UnaryFunction>  template <class BinaryFunction>
1453  inline  inline
1454  double  double
1455  Data::algorithm(UnaryFunction operation) const  Data::algorithm(BinaryFunction operation, double initial_value) const
1456  {  {
1457    if (isExpanded()) {    if (isExpanded()) {
1458      DataExpanded* leftC=dynamic_cast<DataExpanded*>(m_data.get());      DataExpanded* leftC=dynamic_cast<DataExpanded*>(m_data.get());
1459      EsysAssert((leftC!=0), "Programming error - casting to DataExpanded.");      EsysAssert((leftC!=0), "Programming error - casting to DataExpanded.");
1460      return escript::algorithm(*leftC,operation);      return escript::algorithm(*leftC,operation,initial_value);
1461    } else if (isTagged()) {    } else if (isTagged()) {
1462      DataTagged* leftC=dynamic_cast<DataTagged*>(m_data.get());      DataTagged* leftC=dynamic_cast<DataTagged*>(m_data.get());
1463      EsysAssert((leftC!=0), "Programming error - casting to DataTagged.");      EsysAssert((leftC!=0), "Programming error - casting to DataTagged.");
1464      return escript::algorithm(*leftC,operation);      return escript::algorithm(*leftC,operation,initial_value);
1465    } else if (isConstant()) {    } else if (isConstant()) {
1466      DataConstant* leftC=dynamic_cast<DataConstant*>(m_data.get());      DataConstant* leftC=dynamic_cast<DataConstant*>(m_data.get());
1467      EsysAssert((leftC!=0), "Programming error - casting to DataConstant.");      EsysAssert((leftC!=0), "Programming error - casting to DataConstant.");
1468      return escript::algorithm(*leftC,operation);      return escript::algorithm(*leftC,operation,initial_value);
1469    }    }
1470    return 0;    return 0;
1471  }  }
# Line 1219  Data::algorithm(UnaryFunction operation) Line 1478  Data::algorithm(UnaryFunction operation)
1478    rank 0 Data object.    rank 0 Data object.
1479    Calls escript::dp_algorithm.    Calls escript::dp_algorithm.
1480  */  */
1481  template <class UnaryFunction>  template <class BinaryFunction>
1482  inline  inline
1483  Data  Data
1484  dp_algorithm(const Data& data,  Data::dp_algorithm(BinaryFunction operation, double initial_value) const
              UnaryFunction operation)  
1485  {  {
1486    Data result(0,DataArrayView::ShapeType(),data.getFunctionSpace(),data.isExpanded());    Data result(0,DataArrayView::ShapeType(),getFunctionSpace(),isExpanded());
1487    if (data.isExpanded()) {    if (isExpanded()) {
1488      DataExpanded* dataE=dynamic_cast<DataExpanded*>(data.m_data.get());      DataExpanded* dataE=dynamic_cast<DataExpanded*>(m_data.get());
1489      DataExpanded* resultE=dynamic_cast<DataExpanded*>(result.m_data.get());      DataExpanded* resultE=dynamic_cast<DataExpanded*>(result.m_data.get());
1490      EsysAssert((dataE!=0), "Programming error - casting data to DataExpanded.");      EsysAssert((dataE!=0), "Programming error - casting data to DataExpanded.");
1491      EsysAssert((resultE!=0), "Programming error - casting result to DataExpanded.");      EsysAssert((resultE!=0), "Programming error - casting result to DataExpanded.");
1492      escript::dp_algorithm(*dataE,*resultE,operation);      escript::dp_algorithm(*dataE,*resultE,operation,initial_value);
1493    } else if (data.isTagged()) {    } else if (isTagged()) {
1494      DataTagged* dataT=dynamic_cast<DataTagged*>(data.m_data.get());      DataTagged* dataT=dynamic_cast<DataTagged*>(m_data.get());
1495      DataTagged* resultT=dynamic_cast<DataTagged*>(result.m_data.get());      DataTagged* resultT=dynamic_cast<DataTagged*>(result.m_data.get());
1496      EsysAssert((dataT!=0), "Programming error - casting data to DataTagged.");      EsysAssert((dataT!=0), "Programming error - casting data to DataTagged.");
1497      EsysAssert((resultT!=0), "Programming error - casting result to DataTagged.");      EsysAssert((resultT!=0), "Programming error - casting result to DataTagged.");
1498      escript::dp_algorithm(*dataT,*resultT,operation);      escript::dp_algorithm(*dataT,*resultT,operation,initial_value);
1499    } else if (data.isConstant()) {    } else if (isConstant()) {
1500      DataConstant* dataC=dynamic_cast<DataConstant*>(data.m_data.get());      DataConstant* dataC=dynamic_cast<DataConstant*>(m_data.get());
1501      DataConstant* resultC=dynamic_cast<DataConstant*>(result.m_data.get());      DataConstant* resultC=dynamic_cast<DataConstant*>(result.m_data.get());
1502      EsysAssert((dataC!=0), "Programming error - casting data to DataConstant.");      EsysAssert((dataC!=0), "Programming error - casting data to DataConstant.");
1503      EsysAssert((resultC!=0), "Programming error - casting result to DataConstant.");      EsysAssert((resultC!=0), "Programming error - casting result to DataConstant.");
1504      escript::dp_algorithm(*dataC,*resultC,operation);      escript::dp_algorithm(*dataC,*resultC,operation,initial_value);
1505    }    }
1506    return result;    return result;
1507  }  }

Legend:
Removed from v.97  
changed lines
  Added in v.286

  ViewVC Help
Powered by ViewVC 1.1.26