  Pseudodisks are dense structures formed perpendicular to the direction of the
magnetic field during the gravitational collapse of a molecular cloud core.
Numerical simulations of the formation of pseudodisks are usually
computationally expensive with conventional CPU codes. To demonstrate the
proof-of-concept of a fast computing method for this numerically costly
problem, we explore the GPU-powered MHD code Astaroth, a 6th-order finite
difference method with low adjustable finite resistivity implemented with sink
particles. The formation of pseudodisks is physically and numerically robust
and can be achieved with a simple and clean setup for this newly adopted
numerical approach for science verification. The method's potential is
illustrated by evidencing the dependence on the initial magnetic field strength
of specific physical features accompanying the formation of pseudodisks, e.g.
the occurrence of infall shocks and the variable behavior of the mass and
magnetic flux accreted on the central object. As a performance test, we measure
both weak and strong scaling of our implementation to find most efficient way
to use the code on a multi-GPU system. Once suitable physics and
problem-specific implementations are realized, the GPU-accelerated code is an
efficient option for 3-D magnetized collapse problems.
