Changeset ac6d04 for src/solver/solver_regular.cpp
- Timestamp:
- Apr 10, 2012, 1:55:49 PM (14 years ago)
- Children:
- a40eea
- Parents:
- d24c2f
- File:
-
- 1 edited
-
src/solver/solver_regular.cpp (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
src/solver/solver_regular.cpp
rd24c2f rac6d04 38 38 #endif 39 39 40 this->Realloc(rhs.Global(). SizeGlobal().Product());40 this->Realloc(rhs.Global().GlobalSize().Product()); 41 41 42 42 for (grid_iter = rhs.Iterators().Local().Begin(); grid_iter != rhs.Iterators().Local().End(); ++grid_iter) { 43 43 44 mat_index = rhs.GlobalLinearIndex(*grid_iter + rhs.Global(). BeginLocal());44 mat_index = rhs.GlobalLinearIndex(*grid_iter + rhs.Global().LocalBegin()); 45 45 46 46 assert(mat_index >= 0 && mat_index<this->Size()); … … 56 56 for (stencil_iter = A.begin(); stencil_iter != A.end(); ++stencil_iter) { 57 57 58 mat_index2 = rhs.GlobalLinearIndex(*grid_iter + rhs.Global(). BeginLocal() + stencil_iter->Disp());58 mat_index2 = rhs.GlobalLinearIndex(*grid_iter + rhs.Global().LocalBegin() + stencil_iter->Disp()); 59 59 60 60 assert(mat_index2 >= 0 && mat_index2<this->Size()); … … 69 69 for (grid_iter = rhs.Iterators().Boundary1()[i].Begin(); grid_iter != rhs.Iterators().Boundary1()[i].End(); ++grid_iter) { 70 70 71 mat_index = rhs.GlobalLinearIndex(*grid_iter + rhs.Global(). BeginLocal());71 mat_index = rhs.GlobalLinearIndex(*grid_iter + rhs.Global().LocalBegin()); 72 72 73 73 assert(mat_index >= 0 && mat_index<this->Size()); … … 84 84 for (grid_iter = rhs.Iterators().Boundary2()[i].Begin(); grid_iter != rhs.Iterators().Boundary2()[i].End(); ++grid_iter) { 85 85 86 mat_index = rhs.GlobalLinearIndex(*grid_iter + rhs.Global(). BeginLocal());86 mat_index = rhs.GlobalLinearIndex(*grid_iter + rhs.Global().LocalBegin()); 87 87 88 88 assert(mat_index >= 0 && mat_index<this->Size()); … … 110 110 111 111 for (Grid::iterator iter = sol.Iterators().CompleteGrid().Begin(); iter != sol.Iterators().CompleteGrid().End(); ++iter) { 112 index = sol.GlobalLinearIndex(sol.Global(). BeginLocal() + *iter - offset);112 index = sol.GlobalLinearIndex(sol.Global().LocalBegin() + *iter - offset); 113 113 sol(*iter) = this->Sol(index); 114 114 }
Note:
See TracChangeset
for help on using the changeset viewer.
