[Trilinos-Users] Nonlinear maps & VectorToMatrixMarketFile
Klosiewicz Przemyslaw
Przemyslaw.Klosiewicz at ua.ac.be
Sun May 8 10:34:48 MDT 2011
Any ideas?
Greetings,
Przemyslaw
On 29 Apr 2011, at 19:18, Klosiewicz Przemyslaw wrote:
> Hi all,
>
> Up till now I've been using EpetraExt::VectorToMatrixMarketFile(...) and similar to quickly dump distributed vectors to files. However, as soon as "custom-defined" maps for these vectors are used the family of the above IO functions generates results suggesting the requirement of a default linear map.
> As an example, consider an 'interleaving' map:
>
> int myGlobalElements[2];
> if (comm.MyPID() == 0) {
> myGlobalElements[0] = 0;
> myGlobalElements[1] = 2;
> } else {
> myGlobalElements[0] = 1;
> myGlobalElements[1] = 3;
> }
> Epetra_Map map(4, 2, myGlobalElements, 0, comm);
>
> and a vector:
>
> Epetra_Vector v(map);
> for (int idx = 0; idx < map.NumMyElements(); idx++) {
> v[idx] = comm.MyPID();
> }
>
> to create the vector [0, 1, 0, 1] when executed on 2 nodes. A quick cout << v indicates that all is fine.
> Writing the file with VectorToMatrixMarketFile("v.mtx", v) though, gives:
>
> %%MatrixMarket matrix array real general
> 4 1
> 0.0000000000000000e+00
> 0.0000000000000000e+00
> 1.0000000000000000e+00
> 1.0000000000000000e+00
>
> The right behavior can of course be obtained by using a suitable Epetra_Import but I was expecting the IO functions not to depend on linear ordering.
> Is this a bug, a feature or am I missing something crucial?
>
> Thanks,
> Przemyslaw
>
>
> --
> Przemyslaw Klosiewicz
> PhD Student @ CoMP Research Group
> Department of Mathematics and Computer Science, University of Antwerp
> Middelheimlaan 1, 2020 Antwerp, Belgium
>
> Office (G2.07): T (+32) 3 265 32 80, F (+32) 3 265 37 77, E Przemyslaw.Klosiewicz at ua.ac.be
> Home: T (+32) 497 24 18 82, E przemek.klosiewicz at gmail.com, Skype przemek_klosiewicz
>
>
>
>
>
>
> _______________________________________________
> Trilinos-Users mailing list
> Trilinos-Users at software.sandia.gov
> http://software.sandia.gov/mailman/listinfo/trilinos-users
>
--
Przemyslaw Klosiewicz
PhD Student @ CoMP Research Group
Department of Mathematics and Computer Science, University of Antwerp
Middelheimlaan 1, 2020 Antwerp, Belgium
Office (G2.07): T (+32) 3 265 32 80, F (+32) 3 265 37 77, E Przemyslaw.Klosiewicz at ua.ac.be
Home: T (+32) 497 24 18 82, E przemek.klosiewicz at gmail.com, Skype przemek_klosiewicz
More information about the Trilinos-Users
mailing list