MLPG

I have few general queries regarding MLPG meshfree method

1. Can the trial and test functions be taken same for solving a problem

2. In the problem if there are zones with different properties, for example different Coeffiencient of Permeability in the aquifers (Plan view-wise), how to handle the interface between the zones and assemble the global matrix

I am currently implementing some meshfree methods in Matlab. I have already implemented the EFG method. Now I'm trying to implement the MLPG mixed collocation method, and here I run into a problem.