File tree Expand file tree Collapse file tree 1 file changed +12
-4
lines changed Expand file tree Collapse file tree 1 file changed +12
-4
lines changed Original file line number Diff line number Diff line change @@ -1965,10 +1965,18 @@ namespace proteus
19651965 ck.gradTrialFromRef (&p_grad_trial_ref.data ()[k*nDOF_trial_element*nSpace],jacInv,p_grad_trial);
19661966 ck_v.gradTrialFromRef (&vel_grad_trial_ref.data ()[k*nDOF_v_trial_element*nSpace],jacInv,vel_grad_trial);
19671967
1968- std::copy_n (p_trial_ref.data () + k*nDOF_trial_element, nDOF_trial_element, p_trial);
1969- std::copy (std::begin (p_grad_trial), std::end (p_grad_trial), std::begin (p_grad_trial_ib));
1970- std::copy_n (vel_trial_ref.data () + k*nDOF_v_trial_element, nDOF_v_trial_element, p_trial);
1971- std::copy (std::begin (vel_grad_trial), std::end (vel_grad_trial), std::begin (vel_grad_trial_ib));
1968+ for (int i=0 ; i < nDOF_trial_element; i++)
1969+ {
1970+ p_trial[i] = p_trial_ref.data ()[k*nDOF_trial_element + i];
1971+ p_grad_trial_ib[i*nSpace + 0 ] = p_grad_trial[i*nSpace+0 ];
1972+ p_grad_trial_ib[i*nSpace + 1 ] = p_grad_trial[i*nSpace+1 ];
1973+ }
1974+ for (int i=0 ; i < nDOF_v_trial_element; i++)
1975+ {
1976+ vel_trial[i] = vel_trial_ref.data ()[k*nDOF_v_trial_element + i];
1977+ vel_grad_trial_ib[i*nSpace + 0 ] = vel_grad_trial[i*nSpace+0 ];
1978+ vel_grad_trial_ib[i*nSpace + 1 ] = vel_grad_trial[i*nSpace+1 ];
1979+ }
19721980
19731981 if (icase == 0 )
19741982 {
You can’t perform that action at this time.
0 commit comments