epicf/ef

example conductind spehre potential

dumbman opened this issue · 1 comments

dumbman commented

Потенциал сферы стал хуже считаться. Это из-за перехода к другому способу расчета поля??

noooway commented

Да, похоже на то.

В Field_solver'е есть параметр max_Jacobi_iterations = 150;

void Field_solver::solve_poisson_eqn_Jacobi( Spatial_mesh &spat_mesh,
                                             Inner_regions_manager &inner_regions )
{
    max_Jacobi_iterations = 150;
    int iter;
    .....

Если его поставить по-больше, например max_Jacobi_iterations = 2500; , результат будет несколько лучше.

Вцелом, нужно наконец-то поменять метод Якоби на что-то более приличное.
Несколько улучшить текущую ситуацию можно, если заменить Jacobi на SOR -- для сходимости будет требоваться меньше итераций.