BUG: surface charges should be negative

gyral_structure/scripts/gs_fit

@@ -101,7 +101,9 @@ else:
     radial = white_radial
     density = white_density
 
-no_dyads = args.wdensity * density + args.wradial * radial + args.wL2 * L2 + args.wsmoothness * smoothness
+no_dyads = args.wdensity * density + args.wradial * radial + args.wL2 * L2
+if args.wsmoothness != 0:
+    no_dyads = no_dyads + args.wsmoothness * smoothness
 with_dyads = no_dyads + args.wdyad * aligned
 
 
@@ -117,7 +119,7 @@ mm_max = ref_img.affine[:3, :3].dot(voxel_max) + ref_img.affine[:3, -1]
 masked_pial = CorticalMesh.read(args.pial)[surf_mask]
 points_pial = masked_pial.vertices[:, masked_pial.faces].mean(1).T
 
-charge_pial = target_density
+charge_pial = -target_density
 
 charge_max = -np.atleast_1d(charge_pial.sum())