From: kronbichler Date: Wed, 26 Feb 2014 07:17:35 +0000 (+0000) Subject: Add test for various ways of 'add' operations and ghost entries X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=536333453d08313cb35fa8ce11082adac8675f70;p=dealii-svn.git Add test for various ways of 'add' operations and ghost entries git-svn-id: https://svn.dealii.org/trunk@32555 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/mpi/parallel_vector_14.cc b/tests/mpi/parallel_vector_14.cc new file mode 100644 index 0000000000..7ec1440295 --- /dev/null +++ b/tests/mpi/parallel_vector_14.cc @@ -0,0 +1,145 @@ +// --------------------------------------------------------------------- +// $Id$ +// +// Copyright (C) 2011 - 2014 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE at +// the top level of the deal.II distribution. +// +// --------------------------------------------------------------------- + + +// check that handling of ghost elements in parallel distributed vectors works +// appropriately when assigning from ghosted to non-ghosted vectors + +#include "../tests.h" +#include +#include +#include +#include +#include +#include + + +void test () +{ + unsigned int myid = Utilities::MPI::this_mpi_process (MPI_COMM_WORLD); + unsigned int numproc = Utilities::MPI::n_mpi_processes (MPI_COMM_WORLD); + + if (myid==0) deallog << "numproc=" << numproc << std::endl; + + + // processor 0 and 1 own 2 indices each, higher processors nothing, all are + // ghosting global elements 1 and 3 + IndexSet local_owned(std::min(numproc*2, 4U)); + if (myid < 2) + local_owned.add_range(myid*2,myid*2+2); + IndexSet local_relevant(local_owned.size()); + local_relevant = local_owned; + local_relevant.add_range(1,2); + if (numproc > 1) + local_relevant.add_range(3,4); + + // run this twice, once where the vectors have called update_ghost_values + // and once where they have not + for (unsigned int run = 0; run < 2; ++run) + { + parallel::distributed::Vector v(local_owned, local_relevant, MPI_COMM_WORLD); + + // set local values + if (myid < 2) + { + v(myid*2)=myid*2.0; + v(myid*2+1)=myid*2.0+1.0; + } + + v.compress(VectorOperation::insert); + + parallel::distributed::Vector w(v), u(v); + u = 0; + + v*=2.0; + v.add(1.0); + + if (run == 1) + { + v.update_ghost_values(); + w.update_ghost_values(); + u.update_ghost_values(); + } + + if (myid < 2) + { + Assert(v(myid*2) == myid*4.0+1, ExcInternalError()); + Assert(v(myid*2+1) == myid*4.0+3.0, ExcInternalError()); + } + + // copy vector content to non-ghosted vectors, manually created. + parallel::distributed::Vector v_dist(local_owned, MPI_COMM_WORLD), + w_dist(v_dist), u_dist(v_dist); + + v_dist = v; + w_dist = w; + u_dist = u; + + u_dist.add(1.0, v_dist, -1.0, w); + + // copy back to a ghosted vector and update ghost values there + u = u_dist; + u.update_ghost_values(); + + if (myid < 2) + { + Assert(u_dist(myid*2) == myid*2.0+1, ExcInternalError()); + Assert(u_dist(myid*2+1) == myid*2.0+2.0, ExcInternalError()); + Assert(u(myid*2) == myid*2.0+1, ExcInternalError()); + Assert(u(myid*2+1) == myid*2.0+2.0, ExcInternalError()); + } + + Assert(u(1) == 2., ExcInternalError()); + if (numproc > 1) + { + if (run == 1) + Assert(v(3) == 7., ExcInternalError()); + Assert(u(3) == 4., ExcInternalError()); + } + + // check l2 norm + const double l2_norm = u.l2_norm(); + if (myid == 0 && run == 1) + deallog << "L2 norm: " << l2_norm << std::endl; + } + + if (myid == 0) + deallog << "OK" << std::endl; +} + + + +int main (int argc, char **argv) +{ + Utilities::System::MPI_InitFinalize mpi_initialization(argc, argv); + + unsigned int myid = Utilities::MPI::this_mpi_process (MPI_COMM_WORLD); + deallog.push(Utilities::int_to_string(myid)); + + if (myid == 0) + { + std::ofstream logfile("output"); + deallog.attach(logfile); + deallog << std::setprecision(4); + deallog.depth_console(0); + deallog.threshold_double(1.e-10); + + test(); + } + else + test(); + +} diff --git a/tests/mpi/parallel_vector_14.mpirun=1.output b/tests/mpi/parallel_vector_14.mpirun=1.output new file mode 100644 index 0000000000..ec655fb783 --- /dev/null +++ b/tests/mpi/parallel_vector_14.mpirun=1.output @@ -0,0 +1,4 @@ + +DEAL:0::numproc=1 +DEAL:0::L2 norm: 2.236 +DEAL:0::OK diff --git a/tests/mpi/parallel_vector_14.mpirun=3.output b/tests/mpi/parallel_vector_14.mpirun=3.output new file mode 100644 index 0000000000..eff429eec6 --- /dev/null +++ b/tests/mpi/parallel_vector_14.mpirun=3.output @@ -0,0 +1,4 @@ + +DEAL:0::numproc=3 +DEAL:0::L2 norm: 5.477 +DEAL:0::OK