Commit 52cb49fa authored by Andreas Zilian's avatar Andreas Zilian
Browse files

Simplify application of ICs.

parent f3a12ab9
......@@ -37,11 +37,10 @@ ana.analyse();
cla; clf;
fig = figure(1); subplot('Position',[0.05 0.05 0.90 0.90]); hold on;
% apply initial conditions
modeshape = real(ana.X(:,1:2:end)); % get some selected (or all) distinct modes
u0 = 0.20 * rand(fep.ndofs,1);
u0 = 0.20 * rand(fep.ndofs,1); % put some (small) random values
v0 = 0.05 * rand(fep.ndofs,1);
% u0 = 0.2 * modeshape(:,1) + 0.10 * modeshape(:,2);
% v0 = 0.0 * modeshape(:,1) + 0.00 * modeshape(:,2);
% u0 = 0.2 * ana.X(:,1) + 0.10 * ana.X(:,3); % put a combination of modes
% v0 = 0.0 * ana.X(:,1) + 0.00 * ana.X(:,3); % 1 (=1st distinct mode) and 3 (=2nd distinct mode)
ana.applyInitialConditions( u0, v0 );
% obtain and plot the amplitude spectrum
omega = abs(imag(ana.L(1:2:end)));
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment