We do not compute M or C explicitly. Rather, we compute only matrix-vector products, e.g. M • p, using a rapid spectrally-accurate Ewald method discussed elsewhere and solve for the particle moments numerically using an iterative linear solver, GMRES. By solving for the dipoles induced by three orthogonal incident field polarizations E 0 = êx , êy , and êz , we construct the entire C tensor.