What I am very wondering is that if rVec is zero-vector (no any value is added), it works very well, but solution is also zero. However, after any value is added, it does not work as the error was shown
--------------------------------------------------- Exception on processing: Iterative method reported convergence failure in step 1 with residual 31.7698 Aborting! ---------------------------------------------------- Definitely, I have the maximum iteration number more than one.! ??????? Seems to me that it's a mystery.. What is wrong? On Tue, September 30, 2008 9:18 pm, Wolfgang Bangerth wrote: > >> Is that wrong? > > No, this looks alright. > > What I would do is print, right before you call the solver, the sizes of > all > blocks of rVec as well as the locally available elements. I.e., do > something > like this: > > std::cout << get_this_mpi_process(mpi_communicator) << " " > << rVec.size() << " " > << rVec.n_blocks() << " " > << rVec.block(0).size() << " " > << rVec.block(1).size() << " " > << rVec.block(0).local_range().first() << " " > << rVec.block(0).local_range().second() << " " > << rVec.block(1).local_range().first() << " " > << rVec.block(1).local_range().second() << " " > << std::endl > > If one of these numbers looks suspicious, you'll need to investigate what > is > going on. > > Best > W. > > -- > ------------------------------------------------------------------------- > Wolfgang Bangerth email: [EMAIL PROTECTED] > www: http://www.math.tamu.edu/~bangerth/ > _______________________________________________
