-- pde_process_ucd.ads   format pde_read_ucd for solution
--
-- nx, ny, nz, nt  number of unique boundary coordinates
-- xg(1..nx), yg(1..ny), zg(1..nz), tg(1..nt) coordinates
--                                            in increasing order
--
with Ada.Text_Io;
use  Ada.Text_Io;
with Pde_Read_Ucd;
use  Pde_Read_Ucd;
with Real_Arrays;
use  Real_Arrays;

package Pde_Process_Ucd is
   Nx, Ny, Nz, Nt, Ndirchilet, Nneumann : Integer;
   type Real_Vectora is access Real_Vector;
   Xg, Yg, Zg, Tg, Ubd, Ubn  : Real_Vectora;

   function Ub2(X:Real; Y:Real) return Real; -- 2D boundary function
   function Ub3(X:Real; Y:Real; Z:Real) return Real; -- 3D boundary function
   function Ub4(X:Real; Y:Real; Z:Real; T:Real) return Real;  -- 4D
   function Dirchilet(I : Integer) return Real;
   function Neumann(I : Integer) return Real;
   function U(X:Real; Y:Real) return Real; -- optional solution for checking
   procedure Pde_Process_Ucd_File(File_Name : String);
end Pde_Process_Ucd;