Error with Custom Spatially Varying Elastic Moduli Material

[josebastiase]

Hi experts,

I'm encountering an error while running an input file for my custom material. The material computes a spatially varying elastic tensor. Not sure what is going on...

Please find the material files, a minimal example, and log files in debug mode.

Files

[GiudGiud]

Hello

my guess without having run your code is that the residualSetup routine for the derived class needs to call the residualSetup for the parent class first.

If it s not that, please run your input in a debugger and get us the exact line for the "out of bounds" error that you receive currently.

Guillaume

[josebastiase]

Im sorry, should be updated now. Could it MooseArray<double>::operator[] the issue?

[GiudGiud]

That operator is indexing into the array of the material properties at the index of the quadrature point. And for some reason, that material property has not been initialized, the size is too small

[GiudGiud]

ok nevermind.

it looks to me like the seg fault is happening when accessing the variable values.

[GiudGiud]

so with that in mind, is residualSetup the right place for what you are doing?

Residual setup is not in a quadrature point loop maybe?
I think you may be wanting to override another material property routine.
Is _cijk supposed to depend on space? Isnt it a constant in the base class?

[josebastiase]

Hi @GiudGiud

Thank you! Yeah the issue was that I didnt need to override residualSetup. I removed it and it worked ;)