extendedNeighborhoodLattices2D.hh

00001 /* This file is part of the Palabos library.
00002  *
00003  * Copyright (C) 2011 FlowKit Sarl
00004  * Avenue de Chailly 23
00005  * 1012 Lausanne, Switzerland
00006  * E-mail contact: contact@flowkit.com
00007  *
00008  * The most recent release of Palabos can be downloaded at 
00009  * <http://www.palabos.org/>
00010  *
00011  * The library Palabos is free software: you can redistribute it and/or
00012  * modify it under the terms of the GNU Affero General Public License as
00013  * published by the Free Software Foundation, either version 3 of the
00014  * License, or (at your option) any later version.
00015  *
00016  * The library is distributed in the hope that it will be useful,
00017  * but WITHOUT ANY WARRANTY; without even the implied warranty of
00018  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00019  * GNU Affero General Public License for more details.
00020  *
00021  * You should have received a copy of the GNU Affero General Public License
00022  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
00023 */
00024 
00028 #ifndef EXTENDED_NEIGHBORHOOD_LATTICES_2D_HH
00029 #define EXTENDED_NEIGHBORHOOD_LATTICES_2D_HH
00030 
00031 #include "latticeBoltzmann/extendedNeighborhoodLattices2D.h"
00032 
00033 namespace plb {
00034 
00035 namespace descriptors {
00036     // D2Q37 ////////////////////////////////////////////////////////////
00037 
00038     template<typename T>
00039     const T D2Q37Constants<T>::invD = (T)1 / (T) d;
00040 
00041     template<typename T>
00042     const int D2Q37Constants<T>::vicinity = 3;
00043 
00044     template<typename T>
00045     const int D2Q37Constants<T>::c
00046         [D2Q37Constants<T>::q][D2Q37Constants<T>::d] =
00047         {
00048             { 0, 0},
00049             {-1, 3}, {-1, 2}, {-1, 1}, {-2, 2}, {-2, 1}, {-3, 1}, {-1, 0}, {-2, 0}, {-3, 0},
00050             {-3,-1}, {-2,-1}, {-1,-1}, {-2,-2}, {-1,-2}, {-1,-3}, { 0,-1}, { 0,-2}, { 0,-3},
00051             
00052             { 1,-3}, { 1,-2}, { 1,-1}, { 2,-2}, { 2,-1}, { 3,-1}, { 1, 0}, { 2, 0}, { 3, 0},
00053             { 3, 1}, { 2, 1}, { 1, 1}, { 2, 2}, { 1, 2}, { 1, 3}, { 0, 1}, { 0, 2}, { 0, 3}
00054         };
00055 
00056     template<typename T>
00057     const int D2Q37Constants<T>::cNormSqr[D2Q37Constants<T>::q] =
00058         { 
00059             0, 
00060             10, 5, 2, 8, 5, 10, 1, 4, 9,
00061             10, 5, 2, 8, 5, 10, 1, 4, 9,
00062             10, 5, 2, 8, 5, 10, 1, 4, 9,
00063             10, 5, 2, 8, 5, 10, 1, 4, 9
00064         };
00065 
00066     template<typename T>
00067     const T D2Q37Constants<T>::t[D2Q37Constants<T>::q] =
00068         {
00069             0.23315066913235250228650,
00070             0.00028341425299419821740, 0.00535304900051377523273, 0.05766785988879488203006, 0.00101193759267357547541, 0.00535304900051377523273, 0.00028341425299419821740, 0.10730609154221900241246, 0.01420821615845075026469, 0.00024530102775771734547,
00071             0.00028341425299419821740, 0.00535304900051377523273, 0.05766785988879488203006, 0.00101193759267357547541, 0.00535304900051377523273, 0.00028341425299419821740, 0.10730609154221900241246, 0.01420821615845075026469, 0.00024530102775771734547,
00072             0.00028341425299419821740, 0.00535304900051377523273, 0.05766785988879488203006, 0.00101193759267357547541, 0.00535304900051377523273, 0.00028341425299419821740, 0.10730609154221900241246, 0.01420821615845075026469, 0.00024530102775771734547,
00073             0.00028341425299419821740, 0.00535304900051377523273, 0.05766785988879488203006, 0.00101193759267357547541, 0.00535304900051377523273, 0.00028341425299419821740, 0.10730609154221900241246, 0.01420821615845075026469, 0.00024530102775771734547
00074         };
00075 
00076     template<typename T>
00077     const T D2Q37Constants<T>::cs2 = 1./ (1.19697977039307435897239*1.19697977039307435897239);
00078 
00079     template<typename T>
00080     const T D2Q37Constants<T>::invCs2 = 1.19697977039307435897239*1.19697977039307435897239;
00081 
00082     template<typename T>
00083     const char D2Q37Descriptor<T>::name[] = "D2Q37";
00084 
00085     template<typename T>
00086     const char ForcedD2Q37Descriptor<T>::name[] = "ForcedD2Q37";
00087 
00088 }  // namespace descriptors
00089 
00090 }  // namespace plb
00091 
00092 #endif