Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042 #ifndef PLAYA_BLOCKEDVECTORSPACEBASEIMPL_HPP
00043 #define PLAYA_BLOCKEDVECTORSPACEBASEIMPL_HPP
00044
00045 #include "PlayaBlockVectorSpaceDecl.hpp"
00046 #include "PlayaVectorSpaceDecl.hpp"
00047
00048 namespace Playa
00049 {
00050 template <class Scalar> inline
00051 int BlockVectorSpaceBase<Scalar>::dim() const
00052 {
00053 int rtn = 0;
00054 for (int i=0; i<this->numBlocks(); i++)
00055 {
00056 rtn += this->getBlock(i).dim();
00057 }
00058 return rtn;
00059 }
00060
00061 template <class Scalar> inline
00062 int BlockVectorSpaceBase<Scalar>::numLocalElements() const
00063 {
00064 int rtn = 0;
00065 for (int i=0; i<this->numBlocks(); i++)
00066 {
00067 rtn += this->getBlock(i).numLocalElements();
00068 }
00069 return rtn;
00070 }
00071
00072 template <class Scalar> inline
00073 bool BlockVectorSpaceBase<Scalar>::isCompatible(const VectorSpaceBase<Scalar>* other) const
00074 {
00075 if (this->numBlocks() != other->numBlocks()) return false;
00076
00077 const BlockVectorSpaceBase<Scalar>* bs
00078 = dynamic_cast<const BlockVectorSpaceBase<Scalar>* >(other);
00079
00080 if (bs == 0) return false;
00081
00082 for (int i=0; i<this->numBlocks(); i++)
00083 {
00084 if (! (this->getBlock(i).isCompatible(bs->getBlock(i))))
00085 return false;
00086 }
00087 return true;
00088 }
00089
00090
00091 template <class Scalar> inline
00092 std::string BlockVectorSpaceBase<Scalar>::description() const
00093 {
00094 std::ostringstream rtn;
00095 rtn << "BlockVS[";
00096 for (int b=0; b<this->numBlocks(); b++)
00097 {
00098 if (b > 0) rtn << ", ";
00099 rtn << this->getBlock(b).description();
00100 }
00101 rtn << "]";
00102 return rtn.str();
00103 }
00104
00105
00106
00107
00108 }
00109
00110 #endif