import mir.math.common: log;
import mir.test: shouldApprox;

1.0.paretoLPDF(1, 3).shouldApprox == log(paretoPDF(1.0, 1, 3));
2.0.paretoLPDF(1, 3).shouldApprox == log(paretoPDF(2.0, 1, 3));
3.0.paretoLPDF(2, 4).shouldApprox == log(paretoPDF(3.0, 2, 4));