Vishal, I don't think anyone has ever tried to use periodic boundary conditions with GPU. The way we apply constraints on the GPU is very different than what is done on the CPU. So, I am not surprised that it doesn't work. I'll add that to my todo list but I have no idea when I will be able to look at this. If you need this capability any time soon, you will probably need to implement it yourself. If you choose to work on this, I can help you.
Best, Bruno On Saturday, December 12, 2020 at 11:35:56 PM UTC-5 [email protected] wrote: > Hello, > > I am facing an error while I am using matrix-free in GPU with periodic > boundary conditions. I am attaching a minimal example that illustrates the > issue I am facing. I am using deal.II -9.3.0-pre. > > The minimal example is derived from step-64 of the tutorials. In this > code, > 1) Create a single element using the hypercube function. > 2) Create HEX27 finite element based dof_handler and also create the > constraint matrices > 3) Create matrix-free objects on the host and GPU. > 4) Create a host input vector compatible with the constraints ( I set the > values at the unconstrained nodes to be its global Id). > 5) Send the input vector from the host to the GPU > 6) Perform a single vmult operation with the Laplace operator on the host > and the GPU. > 7) Send the output from the GPU to the host > 8) Compare the two outputs > > When I ran the code in debug mode on a single MPI task and compared the > two outputs, the values at the unconstrained nodes do not seem to match. > To ensure there are no bugs in my minimal example, I have a periodicBC > flag. When The periodicBC is set to true, periodic + homogeneous Dirichlet > boundary condition are imposed. If it is set to false, an homogeneous > Dirichlet BC is imposed at the interior node. In this case, the output > values do match. This flag affects how the constraint matrix is created and > nothing else. > > I would be very grateful if someone can tell me what mistake I am making. > Any help is greatly appreciated. > > > thanks and regards, > Vishal Subramanian > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/57368055-0def-4f37-8fcc-423b74e84afen%40googlegroups.com.
