module npolydiv2 { use npolydiv real, dimension (nd), allocatable :: tt #% _init (nd, aa) integer, intent (in) :: nd type( nfilter), intent (in) :: aa call npolydiv_init (nd, aa) #% _lop (xx, yy) integer :: stat1, stat2 if (adj) { stat1 = npolydiv_lop (.false.,.false.,yy,tt) stat2 = npolydiv_lop (.true.,.true.,xx,tt) } else { stat1 = npolydiv_lop (.false.,.false.,xx,tt) stat2 = npolydiv_lop (.true.,.true.,yy,tt) } #% _close call npolydiv_close () }