Matlab crashes when trying to calculate ODFs

[erica-gabriele]

Hello!
I'm new to the MTEX ​​community, but with the help of the documentation and some user scripts, I was able to perform some analyses using the MTEX ​​toolbox, but when I try to calculate the ODF maps (even using the GUI 3pt5 version), Matlab closes with the following log: "Illegal instruction detected". I'm using Matlab version R2023a and MTEX ​​version 6.0.0. I'm trying to perform analysis to replace the use of software and use it in my master's dissertation... What could I be doing wrong?

Here's my script:

`%% ================================
%% 1. Import and configure data
%% ================================

% Load EBSD data
% Use import_wizard

% Define crystal symmetry
CS = {... 
  'notIndexed',...
  crystalSymmetry('m-3m', [3.6 3.6 3.6], 'mineral', 'Ni-superalloy', 'color', [0.53 0.81 0.98])};

% plotting convention
pC = plottingConvention(-zvector,xvector);

% Verify that the data was loaded correctly
if isempty(ebsd)
    error('No EBSD data was loaded. Check the file and path.');
end

%% ==============================================
%% 2. Restrict variables
%% ==============================================

% Filter only the Ni-superalloy phase
ebsd = ebsd('Ni-superalloy');

% Check guidelines
disp(ebsd.orientations);

%% ========================
%% 3. Plot IPF map
%% ========================

% Plot IPF map
plot(ebsd, ebsd.orientations, 'micronbar', 'off');

%% ==================================
%% 4. Grains and subgrains
%% ==================================

% Calculate and correct grains
[grains, ebsd.grainId] = calcGrains(ebsd, 'threshold', 10*degree, 'minPixel', 5);
grains = smooth(grains, 5);

% Plot grains with mean orientation
figure;
plot(grains, 'FaceAlpha', 0.5, 'micronbar', 'off');
hold on;
plot(grains, grains.meanOrientation);
legend off;
hold off;
title('Grãos e Orientações Médias');

% Aspect Ratio
figure;
plot(grains, grains.aspectRatio, 'micronbar', 'off');
setColorRange([1 5]);

% Grain orientation spread (GOS)
GOS = grains.GOS ./ degree;
figure;
plot(grains, GOS, 'micronbar', 'off');
mtexColorbar('title', 'GOS in degree');

%% ========================================
%% 5. Grain boundaries
%% ========================================

% Plot grain contours over IPF
figure;
plot(ebsd, ebsd.orientations, 'micronbar', 'off');
hold on;
plot(grains.boundary, 'lineWidth', 2, 'micronbar', 'off');
hold off;
title('Contorno de Grãos sobre IPF');

% Highlight low and high angle contours
[grains, ebsd.grainId] = calcGrains(ebsd, 'threshold', [1*degree, 10*degree], 'minPixel', 5);

%% ==============================
%% 6. Calculation and ODF visualization
%% ==============================

%calculate and plot the ODF maps with different axis directions
psi = calcKernel(grains.meanOrientation);
odf = calcDensity(ebsdS.orientations,'kernel',psi);

h = Miller(1,1,0,odf.CS);
plotPDF(odf,h,'antipodal','silent');
caxis([0,12]);

close all
h = Miller(1,1,1,odf.CS);
plotPDF(odf,h,'antipodal','silent');
caxis([0,12]);

%% ============================
%% 7. Error rating: MAD
%% ============================

% Plot MAD directly on the map
figure;
plot(ebsd, ebsd.mad, 'micronbar', 'off');
mtexColorbar;
title('Distribuição de MAD');

% Filter measurements with MAD < 0.8 (best quality)
ebsd = ebsd(ebsd.mad < 0.8);

% Calculate and correct grains
[grains, ebsd.grainId] = calcGrains(ebsd, 'threshold', 10*degree, 'minPixel', 5);
grains = smooth(grains, 5);

% Plot filtered EBSD
figure;
plot(ebsd, ebsd.orientations, 'micronbar', 'off');
hold on;
plot(grains.boundary, 'lineWidth', 2, 'micronbar', 'off');
hold off;
title('EBSD - Apenas medições com MAD < 0.8');

%% =========================================
%% 8. CSL (Coincidence Site Lattice) contour analysis
%% =========================================

% Grain calculation and smoothing
[grains, ebsd.grainId] = calcGrains(ebsd);
grains = smooth(grains, 5);

% Plot IQ + average grain orientation (ERROR: Unrecognized method, property, or field 'imagequality' for class 'EBSD'.)
figure;
plot(ebsd, log(ebsd.imagequality), 'figSize', 'large');
mtexColorMap black2white;
setColorRange([.5, 5]);
hold on;
plot(grains, grains.meanOrientation, 'FaceAlpha', 0.5, 'linewidth', 3);
hold off;
title('IQ + Grãos');

% Contours Ni-superalloy x Ni-superalloy
gB = grains.boundary('Ni-superalloy', 'Ni-superalloy');

% Select and plot CSL type Σ3 contours
cs = grains.CS; % Certifique-se que 'cs' está definido
gB3 = gB(angle(gB.misorientation, CSL(3, cs)) < 3*degree);
hold on;
plot(gB3, 'lineColor', 'gold', 'linewidth', 3, 'DisplayName', 'CSL 3');
hold off;

% Identify triple points with at least two CSL(3) contours
isCSL3 = angle(gB.misorientation, CSL(3, cs)) < 3*degree;
tPid = sum(isCSL3(grains.triplePoints.boundaryId), 2) >= 2;

% Plot these triple points
hold on;
plot(grains.triplePoints(tPid), 'color', 'red', 'linewidth', 2, 'MarkerSize', 8);
hold off;
title('Pontos Triplos com múltiplos contornos CSL(3)');

% Join grains with CSL boundaries(3)
[mergedGrains, parentIds] = merge(grains, gB3);

% Plotting blended grain contours
hold on;
plot(mergedGrains.boundary, 'linecolor', 'w', 'linewidth', 3);
hold off;
title('Grãos mesclados com fronteiras CSL(3)');

%% =========================================
%% 9. Twinning
%% =========================================
% Plot histogram of disorientation
figure;
plot(gB.misorientation.angle ./ degree, 'linewidth', 2);
mtexColorbar;
close all;
histogram(gB.misorientation.angle ./ degree, 40);
xlabel('misorientation angle (degree)');

% Check which contours are close to a twinning (< 5°)
isTwinning = angle(gB.misorientation, twinning) < 5*degree;
twinBoundary = gB(isTwinning);

% Plots twinning contours
plot(grains, grains.meanOrientation);
hold on;
plot(twinBoundary, 'linecolor', 'w', 'linewidth', 4, 'displayName', 'twin boundary');
hold off;

%% =========================================
%% 10. Slip Systems and Schmid Factor
%% =========================================

% Get the sliding systems
sS = slipSystem.fcc(cs);

% Sets the direction
r = vector3d.Z;

% View the crystal and sliding system
cS = crystalShape.cube(cs);
plot(cS, 'faceAlpha', 0.5);
hold on;
plot(cS, sS, 'facecolor', 'blue', 'label', 'b');
arrow3d(-0.8 * sS.n, 'faceColor', 'black', 'linewidth', 2, 'label', 'n');
arrow3d(0.4 * r, 'faceColor', 'red', 'linewidth', 2, 'label', 'r');
plottingConvention.default3D().setView;
hold off;

% Calculate Schmid factor 
SF = sS.SchmidFactor(r);

% Visualize SF as a function of stress direction
plot(SF);
[SFMax, pos] = max(SF);
annotate(pos);`

[kilir]

Hi Erica,

which part exactly causes the crash?

odf = calcDensity(ebsdS.orientations,'kernel',psi);

?

What OS are you using? Can you try to use the latest Mtex version? If that also fails, can you try to see if there´s any log/information related to the Matlab crash?

Cheers,
Rüdiger