function [ u, v, p ] = uvp_gms ( nu, rho, n, x, y, t )

%*****************************************************************************80
%
%% uvp_gms returns velocity and pressure for the GMS flow.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    19 August 2020
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Jean-Luc Guermand, Peter Minev, Jie Shen,
%    An overview of projection methods for incompressible flows,
%    Computer methods in applied mechanics and engineering, 
%    Volume 105, pages 6011-6045, 2006.
%
%  Input:
%
%    real NU, the kinematic viscosity.
%
%    real RHO, the fluid density.
%
%    integer N, the number of nodes.
%
%    real X(N), Y(N), the coordinates of nodes.
%
%    real T, the current time.
%
%  Output:
%
%    real U(N), V(N), the X and Y velocity.
%
%    real P(N), the pressure.
%
  [ u, dudt, dudx, dudxx, dudy, dudyy, us, ...
    v, dvdt, dvdx, dvdxx, dvdy, dvdyy, vs, ...
    p, dpdt, dpdx, dpdxx, dpdy, dpdyy, ps ] = all_gms ( nu, rho, n, x, y, t );

  return
end

