Dear Developers,
I am working with EEG and MRI.
I read here that there are two possible reasons for this error:
the source space or the localization of your dipoles are outside of the GM or the Brain.
there are holes in the mesh.
I tried to force a grey layer when generating the head model, or to remerge the mesh. None of them work.
Interestingly, when selecting the "brain" option, everything is working.
Do you have any suggestions?
I can send you the T1 file if it is necessary.
Thanks.
Bests,
Hi, thanks for reaching us,
Are you using this option to define the source space?
When you select the "head", the sources are defined all over the head, and thus you will have sources outside of the brain.
Hi!
Yes.
Do you have any idea to solve this issue?
Thanks!
Ádám
Dear tmedani,
I attach some images to clarify my issue.https://neuroimage.usc.edu/brainstorm/Tutorials/TutVolSource .
This is why I need to use the "Head" option: to get a grid for it, but the computation fails.
Best regards,
Can you please copy-paste here the exact error message you get?
Dear Fraincois,
Thanks for the reply!
Here is everything what I get when I opt for "head":
>> brainstorm
BST> Starting Brainstorm:
BST> =================================
BST> Version: 16-Feb-2023
BST> Compiling main interface files...
BST> Deleting old process reports...
BST> Loading configuration file...
BST> Checking internet connectivity... ok
BST> Update available online: 04-Mar-2023
BST> Initializing user interface...
BST> Starting OpenGL engine... hardware
BST> Reading process folder...
BST> Loading current protocol...
BST> =================================
BST> Checking latest online version for duneuro...
BST> Adding plugin duneuro to path: C:\Users\Berki Ádám József\.brainstorm\plugins\duneuro\bst_duneuro
BST> Adding plugin duneuro to path: C:\Users\Berki Ádám József\.brainstorm\plugins\duneuro\bst_duneuro\bin
BST> Emptying temporary directory...
DUNEURO> Writing temporary files to: D:\BRAINSTORM\tmp
DUNEURO> System call: "C:\Users\Berki Ádám József\.brainstorm\plugins\duneuro\bst_duneuro\bin\bst_duneuro_meeg_win64.exe" "D:\BRAINSTORM\tmp\duneuro_minifile.mini"
Reading 3d Gmsh grid...
version 2.2 Gmsh file detected
file contains 516498 nodes
file contains 2862522 elements
number of real vertices = 516498
number of boundary elements = 0
number of elements = 2862522
volume_conductor.time_reading_gmsh = 259.836
tensors = 5
time = 0.0019782
volume_conductor.time_reading_tensors = 0.0019782
volume_conductor.time_reordering_indices = 0.11899
volume_conductor.time = 259.957
degree = 1
element_type = tetrahedron
WARNING: You are using AMG without SuperLU! Please consider installing SuperLU, or set the usesuperlu flag to false to suppress this warning.
time = 0.0009978
time = 1.58603
This is a new model. Compute the transfer Matrix.
solver.electrode_1.time_rhs_assembly = 0
=== matrix setup (max) 4.33056 s
solver.electrode_1.linear_system_solver.time_matrix_setup = 4.33056
solver.electrode_1.linear_system_solver.time_matrix_assembly = 14.631
solver.electrode_1.linear_system_solver.iterations = 98
solver.electrode_1.linear_system_solver.reduction = 9.56699e-011
solver.electrode_1.linear_system_solver.conv_rate = 0.790247
solver.electrode_1.linear_system_solver.time_solution = 56.6679
solver.electrode_1.linear_system_solver.time = 75.6295
solver.electrode_1.linear_system_solver.time = 75.6705
solver.electrode_1.time_solution = 75.6705
solver.electrode_1.time = 75.6705
solver.electrode_2.time_rhs_assembly = 0
solver.electrode_2.linear_system_solver.iterations = 170
solver.electrode_2.linear_system_solver.reduction = 9.45973e-011
solver.electrode_2.linear_system_solver.conv_rate = 0.873041
solver.electrode_2.linear_system_solver.time_solution = 93.989
solver.electrode_2.linear_system_solver.time = 93.99
solver.electrode_2.linear_system_solver.time = 94.03
solver.electrode_2.time_solution = 94.03
solver.electrode_2.time = 94.03
solver.electrode_3.time_rhs_assembly = 0.0010071
solver.electrode_3.linear_system_solver.iterations = 144
solver.electrode_3.linear_system_solver.reduction = 9.46251e-011
solver.electrode_3.linear_system_solver.conv_rate = 0.851901
solver.electrode_3.linear_system_solver.time_solution = 78.06
solver.electrode_3.linear_system_solver.time = 78.061
solver.electrode_3.linear_system_solver.time = 78.101
solver.electrode_3.time_solution = 78.101
solver.electrode_3.time = 78.102
solver.electrode_4.time_rhs_assembly = 0
solver.electrode_4.linear_system_solver.iterations = 148
solver.electrode_4.linear_system_solver.reduction = 9.09263e-011
solver.electrode_4.linear_system_solver.conv_rate = 0.855369
solver.electrode_4.linear_system_solver.time_solution = 79.298
solver.electrode_4.linear_system_solver.time = 79.299
solver.electrode_4.linear_system_solver.time = 79.341
solver.electrode_4.time_solution = 79.341
solver.electrode_4.time = 79.341
solver.electrode_5.time_rhs_assembly = 0
solver.electrode_5.linear_system_solver.iterations = 116
solver.electrode_5.linear_system_solver.reduction = 9.9285e-011
solver.electrode_5.linear_system_solver.conv_rate = 0.81991
solver.electrode_5.linear_system_solver.time_solution = 62.025
solver.electrode_5.linear_system_solver.time = 62.027
solver.electrode_5.linear_system_solver.time = 62.068
solver.electrode_5.time_solution = 62.068
solver.electrode_5.time = 62.068
solver.electrode_6.time_rhs_assembly = 0
solver.electrode_6.linear_system_solver.iterations = 144
solver.electrode_6.linear_system_solver.reduction = 8.13614e-011
solver.electrode_6.linear_system_solver.conv_rate = 0.851008
solver.electrode_6.linear_system_solver.time_solution = 89.8794
solver.electrode_6.linear_system_solver.time = 89.8804
solver.electrode_6.linear_system_solver.time = 89.9224
solver.electrode_6.time_solution = 89.9224
solver.electrode_6.time = 89.9224
solver.electrode_7.time_rhs_assembly = 0
solver.electrode_7.linear_system_solver.iterations = 436
solver.electrode_7.linear_system_solver.reduction = 9.88693e-011
solver.electrode_7.linear_system_solver.conv_rate = 0.948534
solver.electrode_7.linear_system_solver.time_solution = 278.925
solver.electrode_7.linear_system_solver.time = 278.926
solver.electrode_7.linear_system_solver.time = 278.971
solver.electrode_7.time_solution = 278.971
solver.electrode_7.time = 278.971
solver.electrode_8.time_rhs_assembly = 0
solver.electrode_8.linear_system_solver.iterations = 119
solver.electrode_8.linear_system_solver.reduction = 9.83492e-011
solver.electrode_8.linear_system_solver.conv_rate = 0.823959
solver.electrode_8.linear_system_solver.time_solution = 77.8673
solver.electrode_8.linear_system_solver.time = 77.8683
solver.electrode_8.linear_system_solver.time = 77.9153
solver.electrode_8.time_solution = 77.9153
solver.electrode_8.time = 77.9153
solver.electrode_9.time_rhs_assembly = 0.0010026
solver.electrode_9.linear_system_solver.iterations = 123
solver.electrode_9.linear_system_solver.reduction = 9.8517e-011
solver.electrode_9.linear_system_solver.conv_rate = 0.829175
solver.electrode_9.linear_system_solver.time_solution = 79.6459
solver.electrode_9.linear_system_solver.time = 79.6469
solver.electrode_9.linear_system_solver.time = 79.7548
solver.electrode_9.time_solution = 79.7548
solver.electrode_9.time = 79.7558
solver.electrode_10.time_rhs_assembly = 0.0010052
solver.electrode_10.linear_system_solver.iterations = 129
solver.electrode_10.linear_system_solver.reduction = 9.00882e-011
solver.electrode_10.linear_system_solver.conv_rate = 0.835852
solver.electrode_10.linear_system_solver.time_solution = 83.0362
solver.electrode_10.linear_system_solver.time = 83.0372
solver.electrode_10.linear_system_solver.time = 83.0832
solver.electrode_10.time_solution = 83.0832
solver.electrode_10.time = 83.0842
solver.electrode_11.time_rhs_assembly = 0.0009979
solver.electrode_11.linear_system_solver.iterations = 115
solver.electrode_11.linear_system_solver.reduction = 9.69254e-011
solver.electrode_11.linear_system_solver.conv_rate = 0.818324
solver.electrode_11.linear_system_solver.time_solution = 71.5779
solver.electrode_11.linear_system_solver.time = 71.5789
solver.electrode_11.linear_system_solver.time = 71.6229
solver.electrode_11.time_solution = 71.6229
solver.electrode_11.time = 71.6239
solver.electrode_12.time_rhs_assembly = 0
solver.electrode_12.linear_system_solver.iterations = 131
solver.electrode_12.linear_system_solver.reduction = 9.57576e-011
solver.electrode_12.linear_system_solver.conv_rate = 0.838533
solver.electrode_12.linear_system_solver.time_solution = 80.4437
solver.electrode_12.linear_system_solver.time = 80.4447
solver.electrode_12.linear_system_solver.time = 80.4877
solver.electrode_12.time_solution = 80.4877
solver.electrode_12.time = 80.4877
solver.electrode_13.time_rhs_assembly = 0.0010023
solver.electrode_13.linear_system_solver.iterations = 109
solver.electrode_13.linear_system_solver.reduction = 9.60083e-011
solver.electrode_13.linear_system_solver.conv_rate = 0.809272
solver.electrode_13.linear_system_solver.time_solution = 64.7548
solver.electrode_13.linear_system_solver.time = 64.7558
solver.electrode_13.linear_system_solver.time = 64.7968
solver.electrode_13.time_solution = 64.7978
solver.electrode_13.time = 64.7988
solver.electrode_14.time_rhs_assembly = 0
solver.electrode_14.linear_system_solver.iterations = 308
solver.electrode_14.linear_system_solver.reduction = 9.88016e-011
solver.electrode_14.linear_system_solver.conv_rate = 0.927931
solver.electrode_14.linear_system_solver.time_solution = 179.843
solver.electrode_14.linear_system_solver.time = 179.844
solver.electrode_14.linear_system_solver.time = 179.888
solver.electrode_14.time_solution = 179.888
solver.electrode_14.time = 179.888
solver.electrode_15.time_rhs_assembly = 0.0010116
solver.electrode_15.linear_system_solver.iterations = 160
solver.electrode_15.linear_system_solver.reduction = 8.94715e-011
solver.electrode_15.linear_system_solver.conv_rate = 0.865362
solver.electrode_15.linear_system_solver.time_solution = 93.4663
solver.electrode_15.linear_system_solver.time = 93.4673
solver.electrode_15.linear_system_solver.time = 93.5104
solver.electrode_15.time_solution = 93.5104
solver.electrode_15.time = 93.5114
solver.electrode_16.time_rhs_assembly = 0
solver.electrode_16.linear_system_solver.iterations = 154
solver.electrode_16.linear_system_solver.reduction = 9.56359e-011
solver.electrode_16.linear_system_solver.conv_rate = 0.860873
solver.electrode_16.linear_system_solver.time_solution = 90.1958
solver.electrode_16.linear_system_solver.time = 90.1968
solver.electrode_16.linear_system_solver.time = 90.2378
solver.electrode_16.time_solution = 90.2378
solver.electrode_16.time = 90.2378
solver.electrode_17.time_rhs_assembly = 0.0009962
solver.electrode_17.linear_system_solver.iterations = 150
solver.electrode_17.linear_system_solver.reduction = 8.95288e-011
solver.electrode_17.linear_system_solver.conv_rate = 0.857064
solver.electrode_17.linear_system_solver.time_solution = 87.2974
solver.electrode_17.linear_system_solver.time = 87.2984
solver.electrode_17.linear_system_solver.time = 87.3404
solver.electrode_17.time_solution = 87.3404
solver.electrode_17.time = 87.3414
solver.electrode_18.time_rhs_assembly = 0
solver.electrode_18.linear_system_solver.iterations = 156
solver.electrode_18.linear_system_solver.reduction = 9.62794e-011
solver.electrode_18.linear_system_solver.conv_rate = 0.862565
solver.electrode_18.linear_system_solver.time_solution = 90.2676
solver.electrode_18.linear_system_solver.time = 90.2686
solver.electrode_18.linear_system_solver.time = 90.3096
solver.electrode_18.time_solution = 90.3096
solver.electrode_18.time = 90.3096
solver.electrode_19.time_rhs_assembly = 0.0010208
solver.electrode_19.linear_system_solver.iterations = 862
solver.electrode_19.linear_system_solver.reduction = 9.94561e-011
solver.electrode_19.linear_system_solver.conv_rate = 0.973635
solver.electrode_19.linear_system_solver.time_solution = 506.098
solver.electrode_19.linear_system_solver.time = 506.099
solver.electrode_19.linear_system_solver.time = 506.145
solver.electrode_19.time_solution = 506.145
solver.electrode_19.time = 506.146
solver.electrode_20.time_rhs_assembly = 0
solver.electrode_20.linear_system_solver.iterations = 248
solver.electrode_20.linear_system_solver.reduction = 9.81305e-011
solver.electrode_20.linear_system_solver.conv_rate = 0.911264
solver.electrode_20.linear_system_solver.time_solution = 146.281
solver.electrode_20.linear_system_solver.time = 146.283
solver.electrode_20.linear_system_solver.time = 146.324
solver.electrode_20.time_solution = 146.324
solver.electrode_20.time = 146.324
solver.electrode_21.time_rhs_assembly = 0
solver.electrode_21.linear_system_solver.iterations = 149
solver.electrode_21.linear_system_solver.reduction = 9.72239e-011
solver.electrode_21.linear_system_solver.conv_rate = 0.856651
solver.electrode_21.linear_system_solver.time_solution = 86.631
solver.electrode_21.linear_system_solver.time = 86.633
solver.electrode_21.linear_system_solver.time = 86.674
solver.electrode_21.time_solution = 86.674
solver.electrode_21.time = 86.674
solver.electrode_22.time_rhs_assembly = 0.0010049
solver.electrode_22.linear_system_solver.iterations = 134
solver.electrode_22.linear_system_solver.reduction = 9.72702e-011
solver.electrode_22.linear_system_solver.conv_rate = 0.841944
solver.electrode_22.linear_system_solver.time_solution = 78.2856
solver.electrode_22.linear_system_solver.time = 78.2866
solver.electrode_22.linear_system_solver.time = 78.3266
solver.electrode_22.time_solution = 78.3266
solver.electrode_22.time = 78.3276
solver.electrode_23.time_rhs_assembly = 0.0010029
solver.electrode_23.linear_system_solver.iterations = 144
solver.electrode_23.linear_system_solver.reduction = 8.86608e-011
solver.electrode_23.linear_system_solver.conv_rate = 0.851516
solver.electrode_23.linear_system_solver.time_solution = 85.2373
solver.electrode_23.linear_system_solver.time = 85.2393
solver.electrode_23.linear_system_solver.time = 85.2763
solver.electrode_23.time_solution = 85.2763
solver.electrode_23.time = 85.2773
solver.electrode_24.time_rhs_assembly = 0
solver.electrode_24.linear_system_solver.iterations = 197
solver.electrode_24.linear_system_solver.reduction = 9.8341e-011
solver.electrode_24.linear_system_solver.conv_rate = 0.889614
solver.electrode_24.linear_system_solver.time_solution = 117.625
solver.electrode_24.linear_system_solver.time = 117.626
solver.electrode_24.linear_system_solver.time = 117.669
solver.electrode_24.time_solution = 117.669
solver.electrode_24.time = 117.669
solver.electrode_25.time_rhs_assembly = 0.0010015
solver.electrode_25.linear_system_solver.iterations = 2013
solver.electrode_25.linear_system_solver.reduction = 9.98711e-011
solver.electrode_25.linear_system_solver.conv_rate = 0.988626
solver.electrode_25.linear_system_solver.time_solution = 1206.2
solver.electrode_25.linear_system_solver.time = 1206.2
solver.electrode_25.linear_system_solver.time = 1206.24
solver.electrode_25.time_solution = 1206.24
solver.electrode_25.time = 1206.24
solver.electrode_26.time_rhs_assembly = 0
solver.electrode_26.linear_system_solver.iterations = 231
solver.electrode_26.linear_system_solver.reduction = 9.14676e-011
solver.electrode_26.linear_system_solver.conv_rate = 0.904779
solver.electrode_26.linear_system_solver.time_solution = 126.389
solver.electrode_26.linear_system_solver.time = 126.389
solver.electrode_26.linear_system_solver.time = 126.447
solver.electrode_26.time_solution = 126.447
solver.electrode_26.time = 126.447
solver.electrode_27.time_rhs_assembly = 0
solver.electrode_27.linear_system_solver.iterations = 172
solver.electrode_27.linear_system_solver.reduction = 9.76066e-011
solver.electrode_27.linear_system_solver.conv_rate = 0.87458
solver.electrode_27.linear_system_solver.time_solution = 86.3932
solver.electrode_27.linear_system_solver.time = 86.3932
solver.electrode_27.linear_system_solver.time = 86.4244
solver.electrode_27.time_solution = 86.4244
solver.electrode_27.time = 86.4244
solver.electrode_28.time_rhs_assembly = 0
solver.electrode_28.linear_system_solver.iterations = 250
solver.electrode_28.linear_system_solver.reduction = 9.85128e-011
solver.electrode_28.linear_system_solver.conv_rate = 0.911956
solver.electrode_28.linear_system_solver.time_solution = 126.197
solver.electrode_28.linear_system_solver.time = 126.197
solver.electrode_28.linear_system_solver.time = 126.228
solver.electrode_28.time_solution = 126.228
solver.electrode_28.time = 126.228
solver.electrode_29.time_rhs_assembly = 0
solver.electrode_29.linear_system_solver.iterations = 98
solver.electrode_29.linear_system_solver.reduction = 9.56483e-011
solver.electrode_29.linear_system_solver.conv_rate = 0.790245
solver.electrode_29.linear_system_solver.time_solution = 50.9266
solver.electrode_29.linear_system_solver.time = 50.9266
solver.electrode_29.linear_system_solver.time = 50.9735
solver.electrode_29.time_solution = 50.9735
solver.electrode_29.time = 50.9735
solver.electrode_30.time_rhs_assembly = 0
solver.electrode_30.linear_system_solver.iterations = 174
solver.electrode_30.linear_system_solver.reduction = 9.24171e-011
solver.electrode_30.linear_system_solver.conv_rate = 0.875653
solver.electrode_30.linear_system_solver.time_solution = 87.4864
solver.electrode_30.linear_system_solver.time = 87.4864
solver.electrode_30.linear_system_solver.time = 87.5176
solver.electrode_30.time_solution = 87.5176
solver.electrode_30.time = 87.5176
solver.electrode_31.time_rhs_assembly = 0
solver.electrode_31.linear_system_solver.iterations = 167
solver.electrode_31.linear_system_solver.reduction = 9.48411e-011
solver.electrode_31.linear_system_solver.conv_rate = 0.870928
solver.electrode_31.linear_system_solver.time_solution = 83.8043
solver.electrode_31.linear_system_solver.time = 83.8043
solver.electrode_31.linear_system_solver.time = 83.8511
solver.electrode_31.time_solution = 83.8511
solver.electrode_31.time = 83.8511
solver.electrode_32.time_rhs_assembly = 0
solver.electrode_32.linear_system_solver.iterations = 400
solver.electrode_32.linear_system_solver.reduction = 9.94037e-011
solver.electrode_32.linear_system_solver.conv_rate = 0.944047
solver.electrode_32.linear_system_solver.time_solution = 200.528
solver.electrode_32.linear_system_solver.time = 200.528
solver.electrode_32.linear_system_solver.time = 200.56
solver.electrode_32.time_solution = 200.56
solver.electrode_32.time = 200.56
solver.electrode_33.time_rhs_assembly = 0
solver.electrode_33.linear_system_solver.iterations = 91
solver.electrode_33.linear_system_solver.reduction = 9.31328e-011
solver.electrode_33.linear_system_solver.conv_rate = 0.775837
solver.electrode_33.linear_system_solver.time_solution = 45.7983
solver.electrode_33.linear_system_solver.time = 45.7983
solver.electrode_33.linear_system_solver.time = 45.8452
solver.electrode_33.time_solution = 45.8452
solver.electrode_33.time = 45.8452
solver.electrode_34.time_rhs_assembly = 0
solver.electrode_34.linear_system_solver.iterations = 122
solver.electrode_34.linear_system_solver.reduction = 9.90643e-011
solver.electrode_34.linear_system_solver.conv_rate = 0.827941
solver.electrode_34.linear_system_solver.time_solution = 61.0073
solver.electrode_34.linear_system_solver.time = 61.0073
solver.electrode_34.linear_system_solver.time = 61.0386
solver.electrode_34.time_solution = 61.0698
solver.electrode_34.time = 61.0698
solver.electrode_35.time_rhs_assembly = 0
solver.electrode_35.linear_system_solver.iterations = 145
solver.electrode_35.linear_system_solver.reduction = 8.98049e-011
solver.electrode_35.linear_system_solver.conv_rate = 0.852535
solver.electrode_35.linear_system_solver.time_solution = 72.4273
solver.electrode_35.linear_system_solver.time = 72.4273
solver.electrode_35.linear_system_solver.time = 72.4586
solver.electrode_35.time_solution = 72.4586
solver.electrode_35.time = 72.4586
solver.electrode_36.time_rhs_assembly = 0
solver.electrode_36.linear_system_solver.iterations = 168
solver.electrode_36.linear_system_solver.reduction = 8.91199e-011
solver.electrode_36.linear_system_solver.conv_rate = 0.871321
solver.electrode_36.linear_system_solver.time_solution = 84.2893
solver.electrode_36.linear_system_solver.time = 84.2893
solver.electrode_36.linear_system_solver.time = 84.3362
solver.electrode_36.time_solution = 84.3362
solver.electrode_36.time = 84.3362
solver.electrode_37.time_rhs_assembly = 0
solver.electrode_37.linear_system_solver.iterations = 127
solver.electrode_37.linear_system_solver.reduction = 9.4032e-011
solver.electrode_37.linear_system_solver.conv_rate = 0.833776
solver.electrode_37.linear_system_solver.time_solution = 63.4277
solver.electrode_37.linear_system_solver.time = 63.4277
solver.electrode_37.linear_system_solver.time = 63.459
solver.electrode_37.time_solution = 63.459
solver.electrode_37.time = 63.459
solver.electrode_38.time_rhs_assembly = 0
solver.electrode_38.linear_system_solver.iterations = 232
solver.electrode_38.linear_system_solver.reduction = 9.46929e-011
solver.electrode_38.linear_system_solver.conv_rate = 0.905304
solver.electrode_38.linear_system_solver.time_solution = 115.336
solver.electrode_38.linear_system_solver.time = 115.336
solver.electrode_38.linear_system_solver.time = 115.383
solver.electrode_38.time_solution = 115.383
solver.electrode_38.time = 115.383
solver.electrode_39.time_rhs_assembly = 0
solver.electrode_39.linear_system_solver.iterations = 203
solver.electrode_39.linear_system_solver.reduction = 9.49017e-011
solver.electrode_39.linear_system_solver.conv_rate = 0.892539
solver.electrode_39.linear_system_solver.time_solution = 102.668
solver.electrode_39.linear_system_solver.time = 102.668
solver.electrode_39.linear_system_solver.time = 102.699
solver.electrode_39.time_solution = 102.699
solver.electrode_39.time = 102.699
solver.electrode_40.time_rhs_assembly = 0
solver.electrode_40.linear_system_solver.iterations = 125
solver.electrode_40.linear_system_solver.reduction = 9.84051e-011
solver.electrode_40.linear_system_solver.conv_rate = 0.831657
solver.electrode_40.linear_system_solver.time_solution = 62.349
solver.electrode_40.linear_system_solver.time = 62.349
solver.electrode_40.linear_system_solver.time = 62.4115
solver.electrode_40.time_solution = 62.4115
solver.electrode_40.time = 62.4115
solver.electrode_41.time_rhs_assembly = 0
solver.electrode_41.linear_system_solver.iterations = 157
solver.electrode_41.linear_system_solver.reduction = 9.59601e-011
solver.electrode_41.linear_system_solver.conv_rate = 0.863359
solver.electrode_41.linear_system_solver.time_solution = 79.2867
solver.electrode_41.linear_system_solver.time = 79.2867
solver.electrode_41.linear_system_solver.time = 79.3179
solver.electrode_41.time_solution = 79.3179
solver.electrode_41.time = 79.3179
solver.electrode_42.time_rhs_assembly = 0
solver.electrode_42.linear_system_solver.iterations = 369
solver.electrode_42.linear_system_solver.reduction = 9.81079e-011
solver.electrode_42.linear_system_solver.conv_rate = 0.939458
solver.electrode_42.linear_system_solver.time_solution = 185.721
solver.electrode_42.linear_system_solver.time = 185.736
solver.electrode_42.linear_system_solver.time = 185.767
solver.electrode_42.time_solution = 185.767
solver.electrode_42.time = 185.767
solver.electrode_43.time_rhs_assembly = 0
solver.electrode_43.linear_system_solver.iterations = 264
solver.electrode_43.linear_system_solver.reduction = 9.45635e-011
solver.electrode_43.linear_system_solver.conv_rate = 0.916282
solver.electrode_43.linear_system_solver.time_solution = 132.369
solver.electrode_43.linear_system_solver.time = 132.369
solver.electrode_43.linear_system_solver.time = 132.432
solver.electrode_43.time_solution = 132.432
solver.electrode_43.time = 132.432
solver.electrode_44.time_rhs_assembly = 0
solver.electrode_44.linear_system_solver.iterations = 220
solver.electrode_44.linear_system_solver.reduction = 9.08953e-011
solver.electrode_44.linear_system_solver.conv_rate = 0.900237
solver.electrode_44.linear_system_solver.time_solution = 110.945
solver.electrode_44.linear_system_solver.time = 110.945
solver.electrode_44.linear_system_solver.time = 110.977
solver.electrode_44.time_solution = 110.977
solver.electrode_44.time = 110.977
solver.electrode_45.time_rhs_assembly = 0
solver.electrode_45.linear_system_solver.iterations = 138
solver.electrode_45.linear_system_solver.reduction = 9.82522e-011
solver.electrode_45.linear_system_solver.conv_rate = 0.846215
solver.electrode_45.linear_system_solver.time_solution = 69.0374
solver.electrode_45.linear_system_solver.time = 69.0374
solver.electrode_45.linear_system_solver.time = 69.0686
solver.electrode_45.time_solution = 69.0686
solver.electrode_45.time = 69.0686
solver.electrode_46.time_rhs_assembly = 0
solver.electrode_46.linear_system_solver.iterations = 128
solver.electrode_46.linear_system_solver.reduction = 9.43067e-011
solver.electrode_46.linear_system_solver.conv_rate = 0.83498
solver.electrode_46.linear_system_solver.time_solution = 64.584
solver.electrode_46.linear_system_solver.time = 64.584
solver.electrode_46.linear_system_solver.time = 64.6309
solver.electrode_46.time_solution = 64.6309
solver.electrode_46.time = 64.6309
solver.electrode_47.time_rhs_assembly = 0
solver.electrode_47.linear_system_solver.iterations = 240
solver.electrode_47.linear_system_solver.reduction = 9.45085e-011
solver.electrode_47.linear_system_solver.conv_rate = 0.908304
solver.electrode_47.linear_system_solver.time_solution = 120.961
solver.electrode_47.linear_system_solver.time = 120.961
solver.electrode_47.linear_system_solver.time = 120.993
solver.electrode_47.time_solution = 120.993
solver.electrode_47.time = 120.993
solver.electrode_48.time_rhs_assembly = 0
solver.electrode_48.linear_system_solver.iterations = 162
solver.electrode_48.linear_system_solver.reduction = 9.30808e-011
solver.electrode_48.linear_system_solver.conv_rate = 0.86712
solver.electrode_48.linear_system_solver.time_solution = 81.7908
solver.electrode_48.linear_system_solver.time = 81.7908
solver.electrode_48.linear_system_solver.time = 81.8377
solver.electrode_48.time_solution = 81.8377
solver.electrode_48.time = 81.8377
solver.electrode_49.time_rhs_assembly = 0
solver.electrode_49.linear_system_solver.iterations = 201
solver.electrode_49.linear_system_solver.reduction = 9.70781e-011
solver.electrode_49.linear_system_solver.conv_rate = 0.89163
solver.electrode_49.linear_system_solver.time_solution = 100.29
solver.electrode_49.linear_system_solver.time = 100.29
solver.electrode_49.linear_system_solver.time = 100.321
solver.electrode_49.time_solution = 100.321
solver.electrode_49.time = 100.321
solver.electrode_50.time_rhs_assembly = 0
solver.electrode_50.linear_system_solver.iterations = 106
solver.electrode_50.linear_system_solver.reduction = 9.09071e-011
solver.electrode_50.linear_system_solver.conv_rate = 0.804025
solver.electrode_50.linear_system_solver.time_solution = 53.1781
solver.electrode_50.linear_system_solver.time = 53.1781
solver.electrode_50.linear_system_solver.time = 53.2094
solver.electrode_50.time_solution = 53.2094
solver.electrode_50.time = 53.2094
solver.electrode_51.time_rhs_assembly = 0
solver.electrode_51.linear_system_solver.iterations = 128
solver.electrode_51.linear_system_solver.reduction = 9.31959e-011
solver.electrode_51.linear_system_solver.conv_rate = 0.834903
solver.electrode_51.linear_system_solver.time_solution = 64.1153
solver.electrode_51.linear_system_solver.time = 64.1153
solver.electrode_51.linear_system_solver.time = 64.2247
solver.electrode_51.time_solution = 64.2247
solver.electrode_51.time = 64.2247
solver.electrode_52.time_rhs_assembly = 0
solver.electrode_52.linear_system_solver.iterations = 882
solver.electrode_52.linear_system_solver.reduction = 9.94576e-011
solver.electrode_52.linear_system_solver.conv_rate = 0.974225
solver.electrode_52.linear_system_solver.time_solution = 444.261
solver.electrode_52.linear_system_solver.time = 444.261
solver.electrode_52.linear_system_solver.time = 444.308
solver.electrode_52.time_solution = 444.308
solver.electrode_52.time = 444.308
solver.electrode_53.time_rhs_assembly = 0
solver.electrode_53.linear_system_solver.iterations = 258
solver.electrode_53.linear_system_solver.reduction = 9.60455e-011
solver.electrode_53.linear_system_solver.conv_rate = 0.914476
solver.electrode_53.linear_system_solver.time_solution = 130.139
solver.electrode_53.linear_system_solver.time = 130.139
solver.electrode_53.linear_system_solver.time = 130.17
solver.electrode_53.time_solution = 130.17
solver.electrode_53.time = 130.17
solver.electrode_54.time_rhs_assembly = 0
solver.electrode_54.linear_system_solver.iterations = 119
solver.electrode_54.linear_system_solver.reduction = 9.69775e-011
solver.electrode_54.linear_system_solver.conv_rate = 0.823862
solver.electrode_54.linear_system_solver.time_solution = 60.2414
solver.electrode_54.linear_system_solver.time = 60.2414
solver.electrode_54.linear_system_solver.time = 60.2727
solver.electrode_54.time_solution = 60.2727
solver.electrode_54.time = 60.2727
solver.electrode_55.time_rhs_assembly = 0
solver.electrode_55.linear_system_solver.iterations = 213
solver.electrode_55.linear_system_solver.reduction = 9.79489e-011
solver.electrode_55.linear_system_solver.conv_rate = 0.897448
solver.electrode_55.linear_system_solver.time_solution = 107.619
solver.electrode_55.linear_system_solver.time = 107.619
solver.electrode_55.linear_system_solver.time = 107.65
solver.electrode_55.time_solution = 107.65
solver.electrode_55.time = 107.65
solver.electrode_56.time_rhs_assembly = 0
solver.electrode_56.linear_system_solver.iterations = 138
solver.electrode_56.linear_system_solver.reduction = 9.64349e-011
solver.electrode_56.linear_system_solver.conv_rate = 0.846101
solver.electrode_56.linear_system_solver.time_solution = 69.3825
solver.electrode_56.linear_system_solver.time = 69.3825
solver.electrode_56.linear_system_solver.time = 69.4294
solver.electrode_56.time_solution = 69.4294
solver.electrode_56.time = 69.4294
solver.electrode_57.time_rhs_assembly = 0
solver.electrode_57.linear_system_solver.iterations = 128
solver.electrode_57.linear_system_solver.reduction = 9.71299e-011
solver.electrode_57.linear_system_solver.conv_rate = 0.835173
solver.electrode_57.linear_system_solver.time_solution = 65.5077
solver.electrode_57.linear_system_solver.time = 65.5077
solver.electrode_57.linear_system_solver.time = 65.5545
solver.electrode_57.time_solution = 65.5545
solver.electrode_57.time = 65.5545
solver.electrode_58.time_rhs_assembly = 0
solver.electrode_58.linear_system_solver.iterations = 190
solver.electrode_58.linear_system_solver.reduction = 7.50631e-011
solver.electrode_58.linear_system_solver.conv_rate = 0.88453
solver.electrode_58.linear_system_solver.time_solution = 96.7641
solver.electrode_58.linear_system_solver.time = 96.7641
solver.electrode_58.linear_system_solver.time = 96.811
solver.electrode_58.time_solution = 96.811
solver.electrode_58.time = 96.811
solver.electrode_59.time_rhs_assembly = 0
solver.electrode_59.linear_system_solver.iterations = 102
solver.electrode_59.linear_system_solver.reduction = 9.61584e-011
solver.electrode_59.linear_system_solver.conv_rate = 0.797616
solver.electrode_59.linear_system_solver.time_solution = 51.5965
solver.electrode_59.linear_system_solver.time = 51.5965
solver.electrode_59.linear_system_solver.time = 51.6278
solver.electrode_59.time_solution = 51.6278
solver.electrode_59.time = 51.6278
solver.electrode_60.time_rhs_assembly = 0
solver.electrode_60.linear_system_solver.iterations = 99
solver.electrode_60.linear_system_solver.reduction = 9.25959e-011
solver.electrode_60.linear_system_solver.conv_rate = 0.791867
solver.electrode_60.linear_system_solver.time_solution = 49.9874
solver.electrode_60.linear_system_solver.time = 49.9874
solver.electrode_60.linear_system_solver.time = 50.0968
solver.electrode_60.time_solution = 50.0968
solver.electrode_60.time = 50.0968
solver.electrode_61.time_rhs_assembly = 0
solver.electrode_61.linear_system_solver.iterations = 131
solver.electrode_61.linear_system_solver.reduction = 9.07485e-011
solver.electrode_61.linear_system_solver.conv_rate = 0.83819
solver.electrode_61.linear_system_solver.time_solution = 67.5086
solver.electrode_61.linear_system_solver.time = 67.5086
solver.electrode_61.linear_system_solver.time = 67.5555
solver.electrode_61.time_solution = 67.5555
solver.electrode_61.time = 67.5555
solver.electrode_62.time_rhs_assembly = 0
solver.electrode_62.linear_system_solver.iterations = 130
solver.electrode_62.linear_system_solver.reduction = 9.71338e-011
solver.electrode_62.linear_system_solver.conv_rate = 0.83749
solver.electrode_62.linear_system_solver.time_solution = 65.0533
solver.electrode_62.linear_system_solver.time = 65.0533
solver.electrode_62.linear_system_solver.time = 65.0846
solver.electrode_62.time_solution = 65.0846
solver.electrode_62.time = 65.0846
solver.electrode_63.time_rhs_assembly = 0
solver.electrode_63.linear_system_solver.iterations = 119
solver.electrode_63.linear_system_solver.reduction = 9.6961e-011
solver.electrode_63.linear_system_solver.conv_rate = 0.823861
solver.electrode_63.linear_system_solver.time_solution = 60.9288
solver.electrode_63.linear_system_solver.time = 60.9445
solver.electrode_63.linear_system_solver.time = 60.9757
solver.electrode_63.time_solution = 60.9757
solver.electrode_63.time = 60.9757
dipole_0.density = sparse
dipole_0.time = 0
dipole_1.density = sparse
dipole_1.time = 0
dipole_2.density = sparse
dipole_2.time = 0
dipole_3.density = sparse
dipole_3.time = 0
dipole_4.density = sparse
dipole_4.time = 0
dipole_5.density = sparse
dipole_5.time = 0
dipole_6.density = sparse
dipole_6.time = 0
dipole_7.density = sparse
dipole_7.time = 0
dipole_8.density = sparse
dipole_8.time = 0
dipole_9.density = sparse
dipole_9.time = 0
dipole_10.density = sparse
dipole_10.time = 0
dipole_11.density = sparse
dipole_11.time = 0
dipole_12.density = sparse
dipole_12.time = 0
dipole_13.density = sparse
dipole_13.time = 0
dipole_14.density = sparse
dipole_14.time = 0
dipole_15.density = sparse
dipole_15.time = 0
dipole_16.density = sparse
dipole_16.time = 0
dipole_17.density = sparse
dipole_17.time = 0
dipole_18.density = sparse
dipole_18.time = 0
dipole_19.density = sparse
dipole_19.time = 0
dipole_20.density = sparse
dipole_20.time = 0
dipole_21.density = sparse
dipole_21.time = 0
dipole_22.density = sparse
dipole_22.time = 0
dipole_23.density = sparse
dipole_23.time = 0
dipole_24.density = sparse
dipole_24.time = 0
dipole_25.density = sparse
dipole_25.time = 0
dipole_26.density = sparse
dipole_26.time = 0
dipole_27.density = sparse
dipole_27.time = 0
dipole_28.density = sparse
dipole_28.time = 0
dipole_29.density = sparse
dipole_29.time = 0
dipole_30.density = sparse
dipole_30.time = 0
dipole_31.density = sparse
dipole_31.time = 0
dipole_32.density = sparse
dipole_32.time = 0
dipole_33.density = sparse
dipole_33.time = 0
dipole_34.density = sparse
dipole_34.time = 0
dipole_35.density = sparse
dipole_35.time = 0
dipole_36.density = sparse
dipole_36.time = 0
dipole_37.density = sparse
dipole_37.time = 0
dipole_38.density = sparse
dipole_38.time = 0
dipole_39.density = sparse
dipole_39.time = 0
dipole_40.density = sparse
dipole_40.time = 0
dipole_41.density = sparse
dipole_41.time = 0
dipole_42.density = sparse
dipole_42.time = 0
dipole_43.density = sparse
dipole_43.time = 0
dipole_44.density = sparse
dipole_44.time = 0
dipole_45.density = sparse
dipole_45.time = 0
dipole_46.density = sparse
dipole_46.time = 0
dipole_47.density = sparse
dipole_47.time = 0
dipole_48.density = sparse
dipole_48.time = 0
dipole_49.density = sparse
dipole_49.time = 0.0468702
dipole_50.density = sparse
dipole_50.time = 0
dipole_51.density = sparse
dipole_51.time = 0
dipole_52.density = sparse
dipole_52.time = 0
dipole_53.density = sparse
dipole_53.time = 0
dipole_54.density = sparse
dipole_54.time = 0
dipole_55.density = sparse
dipole_55.time = 0
dipole_56.density = sparse
dipole_56.time = 0
dipole_57.density = sparse
dipole_57.time = 0
dipole_58.density = sparse
dipole_58.time = 0
dipole_59.density = sparse
dipole_59.time = 0
dipole_60.density = sparse
dipole_60.time = 0
dipole_61.density = sparse
dipole_61.time = 0
dipole_62.density = sparse
dipole_62.time = 0
dipole_63.density = sparse
dipole_63.time = 0
dipole_64.density = sparse
dipole_64.time = 0
dipole_65.density = sparse
dipole_65.time = 0
dipole_66.density = sparse
dipole_66.time = 0
dipole_67.density = sparse
dipole_67.time = 0
dipole_68.density = sparse
dipole_68.time = 0
dipole_69.density = sparse
dipole_69.time = 0
dipole_70.density = sparse
dipole_70.time = 0
dipole_71.density = sparse
dipole_71.time = 0
dipole_72.density = sparse
dipole_72.time = 0
dipole_73.density = sparse
dipole_73.time = 0
dipole_74.density = sparse
dipole_74.time = 0
dipole_75.density = sparse
dipole_75.time = 0
dipole_76.density = sparse
dipole_76.time = 0
dipole_77.density = sparse
dipole_77.time = 0
dipole_78.density = sparse
dipole_78.time = 0
dipole_79.density = sparse
dipole_79.time = 0
dipole_80.density = sparse
dipole_80.time = 0
dipole_81.density = sparse
dipole_81.time = 0
dipole_82.density = sparse
dipole_82.time = 0
dipole_83.density = sparse
dipole_83.time = 0
dipole_84.density = sparse
dipole_84.time = 0
dipole_85.density = sparse
dipole_85.time = 0
dipole_86.density = sparse
dipole_86.time = 0
dipole_87.density = sparse
dipole_87.time = 0
dipole_88.density = sparse
dipole_88.time = 0
dipole_89.density = sparse
dipole_89.time = 0
dipole_90.density = sparse
dipole_90.time = 0
dipole_91.density = sparse
dipole_91.time = 0
dipole_92.density = sparse
dipole_92.time = 0
dipole_93.density = sparse
dipole_93.time = 0
dipole_94.density = sparse
dipole_94.time = 0
dipole_95.density = sparse
dipole_95.time = 0
dipole_96.density = sparse
dipole_96.time = 0
dipole_97.density = sparse
dipole_97.time = 0
dipole_98.density = sparse
dipole_98.time = 0
dipole_99.density = sparse
dipole_99.time = 0
dipole_100.density = sparse
dipole_100.time = 0
dipole_101.density = sparse
dipole_101.time = 0
dipole_102.density = sparse
dipole_102.time = 0
dipole_103.density = sparse
dipole_103.time = 0.015643
dipole_104.density = sparse
dipole_104.time = 0
dipole_105.density = sparse
dipole_105.time = 0
dipole_106.density = sparse
dipole_106.time = 0
dipole_107.density = sparse
dipole_107.time = 0
dipole_108.density = sparse
dipole_108.time = 0
dipole_109.density = sparse
dipole_109.time = 0
dipole_110.density = sparse
dipole_110.time = 0
dipole_111.density = sparse
dipole_111.time = 0
dipole_112.density = sparse
dipole_112.time = 0
dipole_113.density = sparse
dipole_113.time = 0
dipole_114.density = sparse
dipole_114.time = 0
dipole_115.density = sparse
dipole_115.time = 0
dipole_116.density = sparse
dipole_116.time = 0
dipole_117.density = sparse
dipole_117.time = 0
dipole_118.density = sparse
dipole_118.time = 0
dipole_119.density = sparse
dipole_119.time = 0
dipole_120.density = sparse
dipole_120.time = 0
dipole_121.density = sparse
dipole_121.time = 0
dipole_122.density = sparse
dipole_122.time = 0
dipole_123.density = sparse
dipole_123.time = 0
dipole_124.density = sparse
dipole_124.time = 0
dipole_125.density = sparse
dipole_125.time = 0
dipole_126.density = sparse
dipole_126.time = 0
dipole_127.density = sparse
dipole_127.time = 0
dipole_128.density = sparse
dipole_128.time = 0
dipole_129.density = sparse
dipole_129.time = 0
dipole_130.density = sparse
dipole_130.time = 0
dipole_131.density = sparse
dipole_131.time = 0
dipole_132.density = sparse
dipole_132.time = 0
dipole_133.density = sparse
dipole_133.time = 0
dipole_134.density = sparse
dipole_134.time = 0
dipole_135.density = sparse
dipole_135.time = 0
dipole_136.density = sparse
dipole_136.time = 0.0937473
dipole_137.density = sparse
dipole_137.time = 0
dipole_138.density = sparse
dipole_138.time = 0
dipole_139.density = sparse
dipole_139.time = 0
dipole_140.density = sparse
dipole_140.time = 0
dipole_141.density = sparse
dipole_141.time = 0
dipole_142.density = sparse
dipole_142.time = 0
dipole_143.density = sparse
dipole_143.time = 0
dipole_144.density = sparse
dipole_144.time = 0
dipole_145.density = sparse
dipole_145.time = 0
dipole_146.density = sparse
dipole_146.time = 0
dipole_147.density = sparse
dipole_147.time = 0
dipole_148.density = sparse
dipole_148.time = 0
dipole_149.density = sparse
dipole_149.time = 0
dipole_150.density = sparse
dipole_150.time = 0
dipole_151.density = sparse
dipole_151.time = 0
dipole_152.density = sparse
dipole_152.time = 0
dipole_153.density = sparse
dipole_153.time = 0
dipole_154.density = sparse
dipole_154.time = 0
dipole_155.density = sparse
dipole_155.time = 0
dipole_156.density = sparse
dipole_156.time = 0
dipole_157.density = sparse
dipole_157.time = 0
dipole_158.density = sparse
dipole_158.time = 0
dipole_159.density = sparse
dipole_159.time = 0
dipole_160.density = sparse
dipole_160.time = 0
dipole_161.density = sparse
dipole_161.time = 0
dipole_162.density = sparse
dipole_162.time = 0
dipole_163.density = sparse
dipole_163.time = 0
dipole_164.density = sparse
dipole_164.time = 0
dipole_165.density = sparse
dipole_165.time = 0
dipole_166.density = sparse
dipole_166.time = 0
dipole_167.density = sparse
dipole_167.time = 0
dipole_168.density = sparse
dipole_168.time = 0
dipole_169.density = sparse
dipole_169.time = 0
dipole_170.density = sparse
dipole_170.time = 0
dipole_171.density = sparse
dipole_171.time = 0
dipole_172.density = sparse
dipole_172.time = 0
dipole_173.density = sparse
dipole_173.time = 0
dipole_174.density = sparse
dipole_174.time = 0
dipole_175.density = sparse
dipole_175.time = 0
dipole_176.density = sparse
dipole_176.time = 0
dipole_177.density = sparse
dipole_177.time = 0
dipole_178.density = sparse
dipole_178.time = 0
dipole_179.density = sparse
dipole_179.time = 0
dipole_180.density = sparse
dipole_180.time = 0.0156118
dipole_181.density = sparse
dipole_181.time = 0
dipole_182.density = sparse
dipole_182.time = 0
dipole_183.density = sparse
dipole_183.time = 0
dipole_184.density = sparse
dipole_184.time = 0
dipole_185.density = sparse
dipole_185.time = 0
dipole_186.density = sparse
dipole_186.time = 0
dipole_187.density = sparse
dipole_187.time = 0
dipole_188.density = sparse
dipole_188.time = 0
dipole_189.density = sparse
dipole_189.time = 0
dipole_190.density = sparse
dipole_190.time = 0
dipole_191.density = sparse
dipole_191.time = 0
dipole_192.density = sparse
dipole_192.time = 0
dipole_193.density = sparse
dipole_193.time = 0
dipole_194.density = sparse
dipole_194.time = 0
dipole_195.density = sparse
dipole_195.time = 0
dipole_196.density = sparse
dipole_196.time = 0
dipole_197.density = sparse
dipole_197.time = 0
dipole_198.density = sparse
dipole_198.time = 0
dipole_199.density = sparse
dipole_199.time = 0
dipole_200.density = sparse
dipole_200.time = 0
dipole_201.density = sparse
dipole_201.time = 0
dipole_202.density = sparse
dipole_202.time = 0
dipole_203.density = sparse
dipole_203.time = 0
dipole_204.density = sparse
dipole_204.time = 0
dipole_205.density = sparse
dipole_205.time = 0
dipole_206.density = sparse
dipole_206.time = 0
dipole_207.density = sparse
dipole_207.time = 0
dipole_208.density = sparse
dipole_208.time = 0
dipole_209.density = sparse
dipole_209.time = 0
dipole_210.density = sparse
dipole_210.time = 0
dipole_211.density = sparse
dipole_211.time = 0
dipole_212.density = sparse
dipole_212.time = 0
dipole_213.density = sparse
dipole_213.time = 0
dipole_214.density = sparse
dipole_214.time = 0
dipole_215.density = sparse
dipole_215.time = 0
dipole_216.density = sparse
dipole_216.time = 0
dipole_217.density = sparse
dipole_217.time = 0
dipole_218.density = sparse
dipole_218.time = 0
dipole_219.density = sparse
dipole_219.time = 0
dipole_220.density = sparse
dipole_220.time = 0
dipole_221.density = sparse
dipole_221.time = 0.0937618
dipole_222.density = sparse
dipole_222.time = 0
dipole_223.density = sparse
dipole_223.time = 0
dipole_224.density = sparse
dipole_224.time = 0
dipole_225.density = sparse
dipole_225.time = 0
dipole_226.density = sparse
dipole_226.time = 0
dipole_227.density = sparse
dipole_227.time = 0
dipole_228.density = sparse
dipole_228.time = 0
dipole_229.density = sparse
dipole_229.time = 0
dipole_230.density = sparse
dipole_230.time = 0
dipole_231.density = sparse
dipole_231.time = 0
dipole_232.density = sparse
dipole_232.time = 0
dipole_233.density = sparse
dipole_233.time = 0
dipole_234.density = sparse
dipole_234.time = 0
dipole_235.density = sparse
dipole_235.time = 0
dipole_236.density = sparse
dipole_236.time = 0
dipole_237.density = sparse
dipole_237.time = 0
dipole_238.density = sparse
dipole_238.time = 0
dipole_239.density = sparse
dipole_239.time = 0
dipole_240.density = sparse
dipole_240.time = 0
dipole_241.density = sparse
dipole_241.time = 0
dipole_242.density = sparse
dipole_242.time = 0
dipole_243.density = sparse
dipole_243.time = 0
dipole_244.density = sparse
dipole_244.time = 0
dipole_245.density = sparse
dipole_245.time = 0
dipole_246.density = sparse
dipole_246.time = 0
dipole_247.density = sparse
dipole_247.time = 0
dipole_248.density = sparse
dipole_248.time = 0
dipole_249.density = sparse
dipole_249.time = 0
dipole_250.density = sparse
dipole_250.time = 0
dipole_251.density = sparse
dipole_251.time = 0
dipole_252.density = sparse
dipole_252.time = 0
dipole_253.density = sparse
dipole_253.time = 0
dipole_254.density = sparse
dipole_254.time = 0
dipole_255.density = sparse
dipole_255.time = 0
dipole_256.density = sparse
dipole_256.time = 0
dipole_257.density = sparse
dipole_257.time = 0
dipole_258.density = sparse
dipole_258.time = 0
dipole_259.density = sparse
dipole_259.time = 0
dipole_260.density = sparse
dipole_260.time = 0
dipole_261.density = sparse
dipole_261.time = 0
dipole_262.density = sparse
dipole_262.time = 0
dipole_263.density = sparse
dipole_263.time = 0
dipole_264.density = sparse
dipole_264.time = 0
dipole_265.density = sparse
dipole_265.time = 0
dipole_266.density = sparse
dipole_266.time = 0
dipole_267.density = sparse
dipole_267.time = 0
dipole_268.density = sparse
dipole_268.time = 0
dipole_269.density = sparse
dipole_269.time = 0
dipole_270.density = sparse
dipole_270.time = 0
dipole_271.density = sparse
dipole_271.time = 0
dipole_272.density = sparse
dipole_272.time = 0
dipole_273.density = sparse
dipole_273.time = 0
dipole_274.density = sparse
dipole_274.time = 0
dipole_275.density = sparse
dipole_275.time = 0
dipole_276.density = sparse
dipole_276.time = 0
dipole_277.density = sparse
dipole_277.time = 0
dipole_278.density = sparse
dipole_278.time = 0
dipole_279.density = sparse
dipole_279.time = 0
dipole_280.density = sparse
dipole_280.time = 0
dipole_281.density = sparse
dipole_281.time = 0
dipole_282.density = sparse
dipole_282.time = 0
dipole_283.density = sparse
dipole_283.time = 0
dipole_284.density = sparse
dipole_284.time = 0
dipole_285.density = sparse
dipole_285.time = 0
dipole_286.density = sparse
dipole_286.time = 0
dipole_287.density = sparse
dipole_287.time = 0
dipole_288.density = sparse
dipole_288.time = 0
dipole_289.density = sparse
dipole_289.time = 0
dipole_290.density = sparse
dipole_290.time = 0
dipole_291.density = sparse
dipole_291.time = 0
dipole_292.density = sparse
dipole_292.time = 0
dipole_293.density = sparse
dipole_293.time = 0
dipole_294.density = sparse
dipole_294.time = 0
dipole_295.density = sparse
dipole_295.time = 0
dipole_296.density = sparse
dipole_296.time = 0
dipole_297.density = sparse
dipole_297.time = 0
dipole_298.density = sparse
dipole_298.time = 0
dipole_299.density = sparse
dipole_299.time = 0
dipole_300.density = sparse
dipole_300.time = 0
dipole_301.density = sparse
dipole_301.time = 0
dipole_302.density = sparse
dipole_302.time = 0
dipole_303.density = sparse
dipole_303.time = 0
dipole_304.density = sparse
dipole_304.time = 0
dipole_305.density = sparse
dipole_305.time = 0
dipole_306.density = sparse
dipole_306.time = 0.0937606
dipole_307.density = sparse
dipole_307.time = 0
dipole_308.density = sparse
dipole_308.time = 0
dipole_309.density = sparse
dipole_309.time = 0
dipole_310.density = sparse
dipole_310.time = 0
dipole_311.density = sparse
dipole_311.time = 0
dipole_312.density = sparse
dipole_312.time = 0
dipole_313.density = sparse
dipole_313.time = 0
dipole_314.density = sparse
dipole_314.time = 0
dipole_315.density = sparse
dipole_315.time = 0
dipole_316.density = sparse
dipole_316.time = 0
dipole_317.density = sparse
dipole_317.time = 0
dipole_318.density = sparse
dipole_318.time = 0
dipole_319.density = sparse
dipole_319.time = 0
dipole_320.density = sparse
dipole_320.time = 0
dipole_321.density = sparse
dipole_321.time = 0
dipole_322.density = sparse
dipole_322.time = 0
dipole_323.density = sparse
dipole_323.time = 0
dipole_324.density = sparse
dipole_324.time = 0
dipole_325.density = sparse
dipole_325.time = 0
dipole_326.density = sparse
dipole_326.time = 0
dipole_327.density = sparse
dipole_327.time = 0
dipole_328.density = sparse
dipole_328.time = 0
dipole_329.density = sparse
dipole_329.time = 0
dipole_330.density = sparse
dipole_330.time = 0
dipole_331.density = sparse
dipole_331.time = 0
dipole_332.density = sparse
dipole_332.time = 0
dipole_333.density = sparse
dipole_333.time = 0
dipole_334.density = sparse
dipole_334.time = 0
dipole_335.density = sparse
dipole_335.time = 0
dipole_336.density = sparse
dipole_336.time = 0
dipole_337.density = sparse
dipole_337.time = 0
dipole_338.density = sparse
dipole_338.time = 0
dipole_339.density = sparse
dipole_339.time = 0
dipole_340.density = sparse
dipole_340.time = 0
dipole_341.density = sparse
dipole_341.time = 0
dipole_342.density = sparse
dipole_342.time = 0
dipole_343.density = sparse
dipole_343.time = 0
dipole_344.density = sparse
dipole_344.time = 0
dipole_345.density = sparse
dipole_345.time = 0
dipole_346.density = sparse
dipole_346.time = 0
dipole_347.density = sparse
dipole_347.time = 0
dipole_348.density = sparse
dipole_348.time = 0
dipole_349.density = sparse
dipole_349.time = 0
dipole_350.density = sparse
dipole_350.time = 0
dipole_351.density = sparse
dipole_351.time = 0
dipole_352.density = sparse
dipole_352.time = 0
dipole_353.density = sparse
dipole_353.time = 0
dipole_354.density = sparse
dipole_354.time = 0
dipole_355.density = sparse
dipole_355.time = 0
dipole_356.density = sparse
dipole_356.time = 0
Dune reported error: Dune::Exception [findEntityImpl:/home/juan/bst-duneuro/src/duneuro/duneuro/common/edgehopping.hh:84]: coordinate is outside of the grid, or grid is not convex
status =
-1
***************************************************************************
** Error: Error during the DUNEuro computation, see logs in the command window.
***************************************************************************
I would not need the "head" option if the "brain" option would include the cerebellum.
I understand, there is a missing in-between solution.
committed 12:37PM - 13 Mar 23 UTC
It does not work, the same error message appears. I also tried with forcing a grey layer, it did not help.
https://wetransfer.com/downloads/08ba2124c33f38e377440ebcd5d5903e20230313204445/745e4dd133fe7fad5019313644ee132e20230313204457/313328
This is an MRI file.
Thank you!
It does not work, the same error message appears.
Which error message?
This is an MRI file.
If order for us to reproduce and try to fix your error, we would need a full imported dataset to reproduce the error, not only the original MRI. Could you please do the following:
Right-click on the subject folder > File > Duplicate subject
In the Anatomy tab, for the duplicated subject: delete all the unused files, and keep only the reference MRI (in green), the head surface, the low-resolution cortex surface, the tissues segmentation, and the FEM mesh generated by SimNIBS
In the Functional data tab: keep only one imported folder (no "RAW" folder), containing only the channel file. Delete everything
Right-click on the channel file > Compute head model => Make sure you can still reproduce exactly the same error with this minimal dataset.
Right-click on the subject > File > Export subject
Upload the .zip file somewhere and post the download link here.
Thanks
Thank you for the example data.
Let's see what @tmedani 's recommendations are.
1 Like
I really hope that you will find a solution to include the cerebellum.
Thank you!
Dear Francois,
A colleague of mine solved the issue: extracting the surfaces from FEM and then setting the grey layer as cortex will lead to a source grid which includes the cerebellum as well.
Thanks for the effort!
Best regards,
1 Like
Hi all,
Indeed, you need to use a surface that includes both the cortex and the cerebellum (let's label it S0) and then use it to generate the volume source space.
Before running the forward computation, you can check the placement of the sources, and you can overlay the FEM mesh on it and check if all the dipoles are well positioned.
If you find some dipoles outside the brain area, you can duplicate the surface S0 and then shrink the new surface by a few millimeters, then use this new surface S0' to generate the volume source.
PS: Please also note that the FEM may have some instabilities on the results when the dipoles are located on tissue interfaces with high conductivity gradients. This issue is still under deep investigation in the FEM modeling field.
Hi everyone, I am also encountering the same issue:
Dune reported error: Dune::Exception [findEntityImpl:/home/juan/bst-duneuro/src/duneuro/duneuro/common/edgehopping.hh:84]: coordinate is outside of the grid, or grid is not convex
status =
-1
***************************************************************************
** Error: Error during the DUNEuro computation, see logs in the command window.
***************************************************************************
I made different attempts, e.g. limiting the sources within the inner skull, limiting to the gray, selecting gray as cortex, neck resection, surface and volume models, but none of it works. By visually inspecting the preview of the dipoles, they seem to correctly be placed from the selected surface inwards. I thought there could be a problem with the tessellation performed by SimNIBS, hence I tried to use the iso2mesh function, but doesn't work in any way (2021 version, old version, fast or slow method).
Any clue?
Thank you,
Hi @romboide2004
Can you send me a reduced version of your protocol, including the head model, the source space you use, and the channel position?
Thanks
