thanks, Barry.
I have params like coordinate, particle position, density, viscosity etc, which definitely prefer not to be PetscScalar type. I try to localize PetscScalar. it works with double precision. quadruple precision really gives me running problem. i haven't configured why.
